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ABSTRACT 


A  pipeline  technique  by  Loomis  and  Sinha  has  been  applied  to  the  design  of  recur¬ 
sive  digital  filters.  Recursive  digital  filters  operating  at  hitherto  impossibly  high  rates  can 
be  designed  by  this  technique. 

An  alternate  technique  by  R.  Gnanasekaran  allows  high  speed  implementation  using 
the  state-space  structure  directly.  High  throughput  is  also  achieved  by  use  of  pipelined 
multiply-add  modules.  The  actual  hardware  complexity  will  depend  upon  the  number 
of  pipeline  stages. 

These  techniques  are  used  for  the  design  of  the  IIR  notch  filter  and  finally,  a  com¬ 
parison  of  the  performance  and  complexity  of  these  two  techniques  is  presented. 
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I.  INTRODUCTION 


A.  DIRECT  FORM  STRUCTURE  REALIZATION 

A  scheme  for  high  speed  recursive  digital  filter  realization  was  suggested  by  H.  H. 
Loomis  and  B.  Sinha  [Ref.  1].  The  main  aspect  of  the  design  is  in  modifying  the  filter 
difference  equation  so  that  the  recursive  pan  requires  only  the  Lth  past  output  instead 
of  the  immediate  past  output.  This  removes  the  tight  recursive  requirement  thus  enabling 
pipelining  of  multipliers  and  adders  in  the  direct  form  implementation  resulting  in  high 
throughput.  One  of  the  difficulties  with  this  method  is  that  the  modilication  introduces 
extra  poles  (not  necessarily  stable)  and  zeros  in  the  filter  transfer  function  which  cancel 
each  other  in  the  infinite  arithmetic  implementation.  In  practice  due  to  the  finite  word 
length  effect  this  cancellation  may  not  occur.  This  may  result  in  an  unstable  realization 
even  if  the  original  filter  were  a  stable  one.  However,  if  L  is  chosen  to  be  large  enough, 
then  it  has  been  proven  that  the  new  poles  introduced  will  also  be  within  the  unit  circle 
thus  guaranteeing  stability  even  if  the  pole-zero  cancellation  does  not  occur.  Tliis 
method  is  called  the  direct  form  structure,  or  Use  scalar  structure. 

IL  STATE  SPACE  REALIZATION 

Another  method  by  R.  Gaanasekaran  is  an  alternate  technique  which  relaxes  the 
recursive  nature  of  the  state  equation  (state  space  realization)  (Ref.  2).  That  is,  the 
present  state  vector  computation  does  not  require  the  immediate  past  state  vector.  The 
proposed  technique  allows  high  speed  implementation  using  the  state-spate  structure 
directly.  High  throughput  is  achieved  by  use  of  pipelined  multiply  -add  modules.  The 
actual  hardware  complexity  will  depend  'ipon  the  method  of  pipelining,  though  very*  high 
speed  is  possible  by  increasing  the  number  of  states  in  the  pipelining.  This  method  im¬ 
plements  A1  instead  of  A.  the  sy  stem  matrix.  The  statespace  realization  will  be  stable 
if  the  original  filter  is  stable.  Furthermore,  the  A*  filter  has  improved  tolerance  to  coef¬ 
ficient  error. 

C  OVERVIEW  OF  THESIS 

Chapter  11  describes  the  UR  notch  filter.  The  transfer  function  and  the  coefficients 
that  are  appropriate  for  the  It R  notch  filter  are  developed. 

Chapter  lit  presents  the  realization  of  the  high  speed  UR  notch  filter  by  using  the 
pipeline  technique  of  Loomis  and  Sinha.  Stability  testing,  frequency  response  and  unit 
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sample  response  are  described.  System  designs  for  pipelined  multioly-add  modules  are 
included. 

Chapter  IV  presents  Gnanasckaran's  state  space  realization  which  includes  a  modi¬ 
fied  state  space  equation.  System  designs  for  pipelined  multiply-add  modules  are  de¬ 
scribed. 

Chapter  V  contains  a  comparison  of  the  results  of  the  Loorais/Sinha  pipeline  re¬ 
cursive  design  and  the  Gnanasekaran  state  space  realization  design. 

Appendix  A  contains  the  program  and  data  for  frequency  response  analysis  of  the 
UR  notch  filters  in  Chapter  III. 

Appendix  B  contains  the  program  and  data  for  unit  sample  response  analysis  of  the 
1IR  notch  filters  in  Chapter  III. 


II.  IIR  NOTCH  FILTER  DESIGN 


From  Rowe,  [Ref,  3:  p.  72]  the  IIR  notch  filter  has  the  transfer  function 


m 


b0z 2  +  b\Z  +  b0 

2  * 
OqZ  4 -b^z  +  a^ 


(2.1) 


Dividing  the  numerator  and  denominator  of  (2.1)  by  a 0,  we  have 


H(z) - 


bp  2  Jh  _  b0 
«o  1  +  %  2  +  % 


j 


«b  *  +  °o 


b  b  ci 

Let  ”  =  Bq.  =  B{,  and  =  A0.  Therefore  the  transfer  function  is 


B0z2  +  Bxz  4-  B0 
z“  4-  B\Z  4-  /1q 


(2.2) 


(23) 


From  (2.1)  and  (2.3)  we  can  see  that 

1.  B{  in  the  numerator,  must  be  equal  to  Bx  in  the  denominator;  and, 

2.  the  first  B0  in  the  numerator  must  be  equal  to  the  last  Z?0  in  the  numerator. 

Rewriting  (2.3)  as 


tf(z)  = 


2 

M  +  ~d~  2  + 1 ) 

_ £o _ 

[z2  4-  B}z  4-  A0) 


(2.4) 


the  transfer  function  N(z)  is  of  the  form 


H(z)- 


B0(:  -  t?i)(z  -  n7) 
(z-d,){z-d2)  • 


(2.5) 


Let  the  zeros  be  placed  at  le~\  and  let  the  poles  be  placed  at  Therefore 

the  transfer  function  is 


(2.6) 


Hh) .  B0(z  -  \e-J6)(z  -  \e+je) 

(z  -  JA0  e~je){z  -  JA0  e+J0) 


H{z)  = 


B0jz2  -\e+J6z-  \e~jez +12 
[z2  -  JA^ e+Jdz  -  Ja^ e~jez  +  V^o"] 


2?0[V-2(  —  )z+l] 

~  +/g  T  .-/a 

[z2-2v%  . -)g  +  4J 


r„  ,  5o[z-2cos02+l‘j 

//(,)  „  — - - — - - 

[z*-2N.  .l0  cos  Oz  +  /i0] 


fljiZ  —  250  cos  +  St 


Vw  T* 

-  '  ‘A^  cos  0z  +  A0 


(2.10) 


According  to  (2.3)  \vc  can  see  that 


-2Bq  cos  Oz  m  -2^A0  cos  Oz 


(2.11) 


which  means 


D0  q  v  <1()  • 


(2.12) 


For  our  example  of  the  HR  notch  filter  we  want  the  i.eros  to  be  placed  at 

lc±J0M61lWMr:  ^ 


(2.13) 


and  we  want  the  poles  to  be  placed  at 

0  1.99^^?  t  tom*  . 


(2.14) 


We  can  see  that 


0  =  0.4167U0339k 


/.'o  •  0.9*101 


(2.15) 

(2.16) 


Bq  “  v  *^o  “0.99. 
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(2.17) 


This  realization  is  shown  in  Figure  1  on  page  6. 

From  equation  (2.1*  *) 

,  B0z2  -  2B0  cos  Bz  4-  B0 

H(Z)  m  — - — - . 

z  -  2 *JA0  cos  $z  +  A0 
Substituting  6,  B0  and  A0,  we  then  have 

N  0.99z2  -  2(0.99)  co$(0 .44671 10339n)z  +  0.99 

■—  *  11  1 

z  -  2(0.99)  cos(0.44671 10339^  +  0.9801 

H(y)  =  °"*2  -  °-32"3;  ±  0-99 

z2  —  0.32993;  +0.9801 

0-99  -  0.32993;-'  4-  0.99;"2 

*i\z  J  .  ,  .  ' 

1  —  0.32993;“  4-0.9S0i;"2 

Because  we  have  the  transfer  function  (Ref.  1:  p.  273] 
i-V -hz~2 

from  (2.21)  and  (2.22)  we  can  sec  that  ail  of  coefficients  arc 


«o  =  0.99 

а,  »  -0.32993 
a2  -  0.99 

б,  -  0.32993 
«  -0.9S01 . 


M 

A’U) 


(2.18) 


(2.19) 

(2.20) 

(2.21) 


(2.22) 


(2.23) 


(2.24) 


V(r)  o,99  -  Q.32993^!  4-  Q.99;"2 

•%)  m  1  —  0.32993;“'  4-  O.9S01:"2  ^ 

>'(*)[  1  -  0.32993;"'  4-  0.980 ir2]  =  A'(.*)[0.99  -  0.329932"'  +  0.99;'2]  (2.26) 

F(r)  “ 0.329932"' >'(;) 4  0.9S0l2"2F(;)  - 0.99A'(;)  -  0.329932"'A'(r)  4-  0.99;  "2A'(;)  (2.27) 
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Figure  I .  Foies  and  Zeros  Location 


Y(z)  -  0.99. \‘(z)  -  0.32993i'-1AX?)  +  0.9 Pz  ~JA'(z)  +  0  32993z”1  >'(z) 
-0.9S01z~3F(z) 


(2.28) 


M  *•  0.99.v(«)  -  0.32993.v(w  -  1)  +  0.99^\fl  -  2)  +  0.32993><«  -  1) 
-0.9$0!>{«-2). 


(229) 


This  is  the  basic  second-order  rccu.sivc  filter  (p^-O).  The  number  of  p  will  be  ex¬ 
plained  later  in  Chapter  111  and  the  rc-’zation  is  shown  in  Figure  2  on  page  7. 

We  can  design  this  HR  notch  filter  'o  be  a  high  speed  HR  notch  filter  ’■>  using 
pipeline  techniques  so  that  it  can  be  operated  at  hithero  impossibly  high  rates.  In 
Chapter  III  we  will  apply  Loomis'Sinha  method  and  in  Chapter  IV  Cnanasekaran's 
state  space  realization.  A  comparison  of  tfc;  results  of  these  two  techniques  will  be 
shown  in  Chapter  V. 
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— 0—  - 

Unit  Delay 

Pipeline  Multiplier 

H 0—  - 

Pipeline  Add  Unit 

x(n) 

input 

y(n) 

Output 

Figure  2.  Second  Order  UR  Notch  Filter:  p  -  0 
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III.  REALIZATION  OF  HR  NOTCH  FILTER  BY  USING  THE 
LOOMIS/SINHA  METHOD 

A.  PIPELINE  RECURSIVE  DESIGN 

From  Chapter  II  we  have  the  second  order  HR  notch  filter  that  is  described  by  the 
transfer  function 


tti  — K  +  '  +  *72  2) 

n(2  )  = - : - x- 

(l  -  V  1  - 


(3.1) 


and  the  difference  equation 

y{n) »  b^n  -1)4-  b$'(n  -  2)  4-  ^(h)  4-  a^n  -1)4-  a2x(a  -  2).  (3.2) 

The  computation  of><»)  requires  the  immediate  past  output  y(n  -  1) .  This  filter  can 
be  developed  into  high  speed  HR  notch  filter  by  using  the  pipeline  technique  of  Loomis 
and  Sinha  (Ref.  1).  The  main  aspect  of  the  design  is  in  modifying  this  filter  difference 
equation  so  that  the  recursive  part  requires  only  from  the  Lth  past  output  instead  of  the 
immediate  past  output.  This  means  the  speed  limits  can  be  overcome  by  modifying  this 
difference  equation  so  thatyfn)  depends  not  on  the  immediate  past  output,  ><«  -  1),  but 
on  yin  -  L)  for  some  L  » 1.  This  can  be  accomplished  as  follows: 

For  L  <=»  1  the  difference  equation  is 

y{ n) »  b^in  -1)4-  -  2)  4-  a^n)  4-  4*  «,*(«  -1)4-  a2x(n  -  2).  (3.3) 

We  can  see  that  this  equation  is  a  2*order  difference  equation. 

From  (Ref.  1:  p.  274)  the  difference  equation  has  (p  4-  m)*»order  where  m  is  the  order 
of  the  difference  equation  and  p  is  the  pipeline  delay  used  to  overcome  the  speed  limit. 
Therefore  equation  (3.3)  has  the  number  of  p=0.  We  can  call  this  equation  the  second 
order  UR  notch  filter  with  p  =  0. 

Delaying  equation  (3.3)  by  1  yields 

y(«  -  1)  s*  b^n  -  2)  4-  b^yin  -  3)  4-  -1)4-  u,jc(/»  -  2)  4-  a2jk(«  -  3).  (3.4) 

For  L  ®  2,  putting  (3.4)  into  (3.3)  results  in  the  difference  equation 


$ 


(3.5) 


>•(«)  =  b^b^n  -  2)  +  b&in  -  3)  4-  %*(«  -  1)  +  axx(n  -  2)  +  a2x(rt  -  3)] 

+  bjin  -  2)  +  a0.r(«)  +  axx{n  -  1)  +  a2x(n  -  2) 

y{n)  -  b\y{n  -  2)  +  bxb^\n  -  3)  +  a^b^xin  -  1)  +  axbxx{n  -  2)  +  a2bxx{n  -  3) 
+  btfin  -  2)  +  Oqx(h)  +  axx(n  -  1)  +  a2x{n  -  2) 

j<w)  =  {b2  +  b^in  -  2)  +  b^in  -  3)  +  aox(n)  +  M ,  +  a^JcC/i  -  1) 

+  (axbx  +  a2)x(fi  -  2)  +  a2bxx(n  -  3). 


We  see  that  the  computation  of  y(»)  requires  only  the  second  past  output  y{n  -  2) 
instead  of  immediate  past  output  y(n  - 1),  and  we  can  use  this  technique  for  the 
higher- order  of  L  . 

We  see  that  equation  (3.7)  is  a  3-order  difference  equation  and  has  p«  1.  Thus  we 
can  call  this  filter  the  augmented  HR  notch  filter  for  p=  1. 

The  z  transform  of  equation  (3.7)  is 


}'=  (bx  +  b2)z~2  Y  +  bxb2z~3  Y  +  a0X  +  M»  +  ax)z~]X 
+  (axbx  +  a2)z~2X  +  a2bxz~3X 


y  a9  +  M)  +  +  (axbx  4  a2)z~ 2  +  a2bxz~3 

^  1  -  [b2  4  b):'1  -  bxbiz~3 

-|k  :  2)(l+/>,r  ’) 

//  (X  636 - Zi - Is - 7~  • 

(l -V  )U  +  M  ) 


(3.8) 


(3.9) 


(3.10) 


That  means  the  transfer  function  of  the  augmented  HR  notch  filter,  p«  1,  is  equal 
to  the  transfer  function  of  the  second  order  HR  notch  filter,  p“0,  multiplied  and  divided 
by  the  same  polynomials. 

We  can  design  the  higher-order  difference  equation  by  starting  from  a  second-order 
difference  equation.  That  means  we  can  design  the  higher  order,  p-augmented  UR  notch 
filter  by  starting  from  our  second  order  HR  notch  filter.  p*0,  as  follows: 

From  Chapter  ti  the  second  order  HR  notch  filter  transfer  fimetion  is  the  original 
transfer  function 


11(2-') 


+  ti)-.-1  +  q,i~3) 

(I  -  ' 


(3.11) 
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Multiply  the  numerator  and  denominator  of  the  original  transfer  function  by  the 
same  polynomial  equation.  This  can  be  thought  of  as  equation  (3.12)  As  developed  in 
(Ref.  1:  p.  281],  the  transfer  function  of  the  p-augmented  filter  is  as  shown  in  equation 
(3.13)  We  know  from  other  considerations  that  the  transfer  function  of  the  higher  order 
realization  must  in  fact  be  the  same  as  the  original  and  therefore  must  have  the  same 
poles  and  zeroes  in  the  z  -  plane. 

The  coefficients  of  the  a  polynomial  depend  on  the  original  by  and  b2,  and  p  is  de¬ 
termined  by  the  number  of  stages  in  the  pipeline  multiply  and  add  units. 


+  atz-1  +  Oj2~2)(1  +  a^2~]  H - )-  ap2~p) 

(i  -  v'1  -  v’xi + «“»■' + - + 4'^) 


(3.12) 


«V') 


jV  +  -  + 


(3.13) 


The  values  of  bf\  a f\ and  af'  are  given  in  Tables  1,  2  and  3  (Ref.  1:  pp.  278-279]. 


Table  1.  b  ‘s  for  />  ”  0.1 . 5 


Coeflicient  of  />-' 

Coelficient  of  />** 

r 

AJS. 

bfli 

0 

l>\ 

b, 

i 

/>f  +  l>; 

bA 

2 

bl  + 

bA  +  b] 

b*  +  )bfl>:  +  bj 

b\b}  4-  2b A 

b\  +  #(/»,  +  U\b] 

bibj  +  3 />(/']  r  h) 

5 

l<\  +  5  by>:  +  ()b-b}  +  b\ 

b>b}  +  4  bA  +  3 
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Table  2.  ot  s  for  p  -  0,1,— ,5 


EH 

ES 

af 

otf> 

a? 

o#> 

0 

1 

0 

0 

0 

0 

0 

1 

I 

b\ 

0 

0 

0 

0 

2 

1 

b\ 

6?4  62 

0 

0 

0 

3 

1 

6, 

6?4  63 

6?  +  26, 62 

0 

:eh 

1 

EM 

6?  4  62 

6?  +  26,6, 

6;  4-  36J6,  4-  if 

0 

5 

1 

bx 

6?4-62 

6J  +  26,6, 

6f4-36?6,4-6? 

6?  4-46?6,4-36,6j 

Table  3.  a  ‘s  for  p  —  0,1,..., 3 


n 

HH 

«?> 

<#> 

a 

E9 

BHH 

0 

0 

0 

D 

El 

+  o, 

«,6,  4-  a, 

a.6, 

0 

0 

2 

«o 

«g6,  + 

fl#(6?  +  6.0 

4-  +  fl, 

«i(6?  4-  6,) 

4-  <r26, 

<6(6?  4-  60 

0 

H 

4> 

«o(6f  4-  60 

4-  £?,/»,  +  a. 

«c(6?  4-  26,60 

4*  «,(6(  4*  60 

4  flj6, 

4-  fl,(6,  4-  26,60 

4-  t»:(6f  4-  6,) 

4-  <jj(6J  4  26,60 

(3.14) 

(3.15) 

(3.16) 


IL  SECOND  ORDER  UR  NOTCH  FILTER;  P  •=  0 

From  equation  (3.13)  in  section  A,  we  let  p  =>  0.  This  equation  becomes: 

jru-*)  -ia±a£l±aO. 


F(;)  _  (<tp  4-  -r  a?;~2) 

%)  "  (1  - 6,2“'  -  ^r“J) 


Manipulating  (3.16)  will  result  in  the  difference  equation,  (3.20). 
>W  -  -  6Ir“13  “  A'(2)Coo  4  «,r*  +  «v"23 


(3.17) 


F(z)  -  bxz~l  Y(z)  -  b2z~2Y(z)  =  a^X{z)  +  axz~xX{z)  +  a2z~2X(z) 
Y(z)  =  %T(z)  +  axz~xX{z)  +  a2z~2X(z)  +  bxz~x  Y(z)  +  b2z~ 2  Y(z) 
y{n)  =  aox(rt)  +  axx(n  -  1)  +  a2x(n  -  2)  +  b^y{n  -  1)  -f  b^n  -  2). 


(3.18) 


(3.19) 

(3.20) 


This  realization  is  shown  in  Figure  3  on  page  13. 

From  chapter  II  we  know  that  for  our  example  notch  filter,  that  the  coefficient 
values  are 


Oq  =  0.99 
a,  -  -0.32993 
a2 «  0.99 
bx  =0.32993 
&2  =  -0.9801 . 


The  transfer  function  is 

(1  —  bxt~x  —  hz~') 
The  difference  equation  is 


(3-21) 


(3.22) 


>(«) «  %v(m)  +  axx(n  -  1)  +  a2x(n  -  2)  +  b^{n  -  1)  +  b^in  -  2). 


(3,23) 


C.  AUGMENTED  HR  NOTCH  FILTER:  P  *  1 

From  equations  (3.12)  and  (3.13)  in  section  A.  we  let  p  « 

come: 

1.  Those  equations  be* 

f/*(r')» 

(ttp  +  rqr"1  4*  <ijr“J)(l  +  c^V) 

(3.24) 

(l-M“,-M"aMi  +  «!V‘) 

tf'(r')  - 

(3.25) 

i-*SV2-*<V3 

//'(2“l)  - 

Hr) 

A*U) 

(3.26) 

Figure  3.  Second  Order  HR  Notch  Filter:  p  »  0 


m)  _  oS'^alVUaiVHaSV3 
AW  “  i  _  -  4JV* 


Manipulating  (3.27)  will  result  in  the  difference  equation,  (3.31). 

\\:)l\  -  <4V2  -  ^jV3]  -  AV’JCoJ0  +  ajV‘  -f  a^V2  +  uJV3] 

Y(s)  -  bi'^Yiz)  -  &5VJ>'(r) «  eg'#*)  +  ajV’A'W  +  a}VJAW 

+  o},)r,A'(i) 

Y(s) -  a5u.V(r)  +  «jV*JIW  «f  «iVJAW  4-  ajV3A'(;)  +  6j,>x“al W 

+  *SVJ>W 


><«)  -  Oo'^rt)  +  a{3>jc<ft  -1)  +  flj'W  -  2)  +  aj’lxf/i  -  3)  +  “  2) 

+  *?!><* -3). 


(3.27) 

(3.2S) 

(3.29) 

(3.30) 

(3.31) 


This  realization  is  shown  in  Figure  4  on  page  14. 

From  Tables  1,  2  and  3  (Ref.  1:  pp.  278-279],  we  have 


Figure  4.  The  Augmented  HR  Notch  Filter:  p  *  1 


Bo'*  -  1 

af 1 »  h,  «  0.32993 
^  -  0.99 

a}’1  -  ooh,  -f  a,  «  -0.0032993 
aj13  ■»  +  Oj  m  0.881 1462 

**0.3266307 

6j‘>-h?  + ^--0.8712462 
hj»  -0.32336439. 

The  transfer  function  is 


mr') 


4,,  +  B!VUaiV3  +  njV3 
i-hSVJ-65V3 


(3.32) 


(3.33) 


The  difference  equation  is 


(3.34) 


y{n)  =  af  ^(w)  +  a?]x{n  -  1)  +  af*(«  -  2)  +  a\l)x(n  -  3)  +  b[])y{n  -  2) 

+  lg%  -  3). 

D.  AUGMENTED  IIR  NOTCH  FILTER:  P  =  2 

From  equations  (3.12)  and  (3.13)  in  section  A,  we  let  p  =  2.  Those  equations  be¬ 
come: 


ff.  -k  _  (H  +  <y~'  +  +  njV1  +  a fi~a) 

(1  -  b,z~:  -  I  +  af’j"1  +  a®*-1) 


rnr'i 


af  +  ai V^+agr’  +  af.--1  +  a«V 


I  -  b?r3  -  if.-' 


W(*  '-1W 


)W  a®  -t-  +  «.-» +  af +  „f '  ■-* 

*<-)"  I -if'.--* -if..-1 

Manipulating  (3.3S)  will  result  in  the  difference  equation*  (3,42). 

>\z)i\  -  -  if**-4]  -  a*u)[4,) + ajV‘ + af r2 + af  r3  +  ^p*,-4] 

>V)  -  *>jV3  >(.-)  -  -  afA(x)  +  afr’A'u)  +  af  )rJA(s) 

+  afrJ.V(r)-faf,5“4A'U) 

>V)  -  «f  *Y(r)  +  a}V*A'{;)  +  af  x“2A(x)  4*  ajV'Afc)  4-  «f x^AV) 

+  bf):-)\X:)  +  b?s'*ns) 


(3.35) 

(3.36) 

(3.3?) 

(MS) 

(3.39) 

(3.40) 

(3.41) 


yin)  m  of 4»)  +  af a<h  -  1)  +  afjtfa  -  2)  +  af 4*  -  3)  +  qfjr(i*  -  4) 
+  ^i.-3)  +  6ia,>(«-4). 

This  realization  is  shown  in  Figure  5  on  page  16. 

From  Tables  1,  2  and  3  (Ref.  1:  pp.278-279),  wc  have 


15 


16 


af  =  bx  =0.32993 

af  =  b\  +  i2  =  -0.8712462 

of  =  *0  =  0.99 

*f  -  Oq^  +  a,  =  -0.0032993 

«f  —  +  b 2)  +  +  a2  =  0.01861246 

af  =  aj(if  +  b2)  +  a2ij  =  0.61408096 

a?  -  a2(6f  +  i2)  =  -0.86253373 

if  =  ^  +  2^  =  -0.61081465 

if  =  b\b2  4-  i2  =  0.8539084  . 


The  transfer  function  is 


J2)  ,  „<2),-l  ,  <2)  — 2  ,  J2)„ -3  (2)  -4 

^0  +  ^1  -  +  q2  2  +  a3  *  +  a4  2 


1  -  MV3  -  bpz-4 


The  difference  equation  is 

y{n)  -  ctQ’x{n)  +  af\x{n  -  1)  +  «f  x(w  -  2)  +  af\x{n  -  3) 
+  «f *(«  -  4)  +  bf]y(n  -  3)  +  b^y(n  -  4). 


(3.44) 


(3.45) 


E.  AUGMENTED  UR  NOTCH  FILTER:  P  =  3 

From  equations  (3.12)  and  (3.13)  in  section  A,  we  let  p  a  3,  Those  equations  be¬ 
come: 


_k  (flp  +  fltz  1  +  (^  *)(l+ofz  1  +  «f;  2  4  ot^r  3) 
U  j~  (l-V^-^Kl+alV^^  +  ^V3) 


**  \2  ) .  m  -4  m  -5 


-5 


i-ifr4-#*- 


1W 


m 


Y(z)  _  af  +  a<V>  +  «<V2  +  a<V3  +  *f c"4  +  ojV5 
1  -  if-’-4  -  if r5 


(3.46) 


(3.47) 


(3.4S) 


(3.49) 
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Manipulating  (3.49)  will  result  in  the  difference  equation,  (3.53). 


Y(z)[l  -  if z-4  -  ifz~5']  - 

*(z)[af  +  *SVl  +  a<  V2  +  af z~3  +  af z"4  +  af  z“5] 

V(z)  -  if z~*y(z)  -  bf}z~5Y(?)  =•  af *(z)  +  aj3>z_1;f(z)  4-  afz“2*(z) 

+  af)z~3^(z)  +  a?)z"4Ar(z) 

4-  af^Xlz) 

Y{z)  =  afjF(z)  +  afz'^z)  +  af z~2J(z)  +  af >z‘3*(z)  +  af z^z) 

+  af}z~5X(z)  +  if z"4  Y(z)  +  if  z'5  Y(z) 

y(rt)  =  af *(<*»)  4-  a[3lr(n  —  1)  +  afv(«  -  2)  +  af *(«  —  3)  +  af  jc(«  -  4) 
4-  af\x{n  -  5)  +  if><w  -  4)  4-  bf]y(n  -  5). 


(3.50) 

(3.51) 

(3.52) 

(3.53) 


This  realization  is  shown  in  Figure  6  on  page  19. 

From  Tables  1,  2  and  3  (Ref.  1:  pp.  278-279],  we  have 

«f=l 

of  =  i,  =  0.32993 

af  -if  +  i,  =  ~0.87 12462 

of  o  if  +  25,  i2 »  -0.61081465 

of  » %  “  0.99 

af  =>  Oyi,  +  rt,  =  — 0.0032993 

af  “  a^i2  4-  i2)  4-  a,5,  4*  «2  »  0.01861246 

af  «  Oo(bl  4-  25, ij)  4-  0,(6?  4-  ij)  4-  a,5,  *  0.00937445 

of  =*  a,(53  4-  25,5^  4-  a2(i2  4*  5,) *»  —0.66100766 

af  *  a2(if  4-  25, i2) «  -0.6047065 

if  -  5?  4-  35, ^ii  4-  5{  -  0.65238232 

if  -  ifo  4-  25,  i3 «  0.59865944 . 


The  transfer  function  is 

of  4-  a}V  4-  a<V2  4-  afz“3  4-  tfV4  4-  a<V5 


WV1)- 


1  -  if z"4  -  if  z“5 


The  difference  equation  is 


(3.54) 


(3.55) 


18 


19 


(3.56) 


y(n)  =  ^ 3lx(«)  4-  *'('i  -  1)  +  -  2) 

+  af\x(n  -  3)  +  af]x{n  -  4) 

+  af]x(n  -  5)  +  bf]y(n  -  4) 

+  bf'An  -  5). 

F.  STABILITY  TESTING 

From  section  A  if  we  start  with  a  stable  filter,  we  would  Tike  to  be  assured  that  the 
poles  of  the  augmented-order  filter  all  be  outside  the  unit  circle  in  the  z~‘  plane.  We  de¬ 
sired  this  because,  even  though  the  added  poles  are  cancelled  by  the  zeroes  of 

(1  +  awz-1  -I - f  awz-'),  realization  imperfections  will  prevent  exact  cancellation  and  the 

augmented  filter  would  be  unstable.  In  order  to  examine  the  stability  question,  we  will 
use  Jury's  stability  test  as  applied  to  successively  higher-degree  denominator  polynomials 
in  z  as  p  increases  [Ref.  1:  p.  281].  The  results  from  this  paper  are  summarized  below: 

1.  Condition  for  Stability  of  the  Original  IIR  Notch  Filter:  p  “  0 

The  original  IIR  filter  is  described  by  equation  (3.14)  and  to  test  for  stability, 
we  analyze  the  following  characteristic  polynomial  resulting  from  it; 

F(z)  =  z2-iV-&2»  0  (3.57) 

The  conditions  for  stability  arc: 

1  -  -  &2  >  0 ...  © 

1  +  -  !>j  >0 ...  ©  (3.58) 

\b2\  <1...© 

These  conditions  are  shown  graphically  in  Figure  7  on  page  21. 

Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 
(3.58) 

Checking  for  stability  of  our  example  IIR  notch  filter;  we  have  from  ©  of(3.5S) 

1  -  6,  -  bi  >  0 

1  —  (0.32993)  -  ( -0.9S0I)  >  0 
1.65017  >0 


from  ©  of  (3.58) 
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Figure  7.  Condition  for  Stability  of  the  Original  Filter  [Ref.  1:  p.  283) 


1  +  6,  -  b2  >  0 
1  +  (0.32993)  -( -0.9801)  >0 
2.31003  >0 


from  ©  of  (3.58) 

I  i  <  l 

[0.9S01 1  <1 
0.9801  <  1 

Thus  the  conditions  are  all  satisfied  and  the  original  filter  is  stable. 

2.  Condition  for  Stability  of  the  Augmented  HR  Notch  Filter:  p- 1 

The  augmented  UR  notch  filter,  p ■“  1,  is  described  by  equation  (3.33)  and  to  test 
for  stability,  we  analyze  the  following  characteristic  polynomial  resulting  from  it: 

FU)-2  +  b,-0  (3.59) 

The  conditions  for  stability  are: 
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1+  b}  >  0 ...  © 
1  -  bx  >  0 ...  © 


(3.60) 


These  conditions  are  shown  graphically  in  Figure  8  on  page  23. 

Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 
(3.60) 

Checking  for  stability  of  our  example  HR  notch  filter;  we  have  from  ©  of  (3.60) 

l  +  6t>0 

1  + 0.32993  >0 
1.32993  >0 

from  ©  of  (3.60) 

l-6t>0 
0.67007  >0 

Thus  the  conditions  are  all  satisfied  and  the  augmented  IIR  notch  filter,  p  =  1,  is  stable. 
3.  Condition  for  Stability  of  the  Augmented  HR  Notch  Filter:  p  =  2 

The  augmented  HR  notch  filter,  p=  2,  is  described  by  equation  (3.44)  and  to  test 
for  stability,  we  analyze  the  following  characteristic  polynomial  resulting  from  it: 

F(r)<=:2  +  bxz  +  (!>?  +  bj  (3.61) 

The  conditions  for  stability  are: 

1  +  b\  +  (l>f  +  b±)  >  0 ...  © 

l-6,  +  (/>,2  +  &2)>0...©  (3.62) 

U>,2  +  *2|  <!...© 


These  conditions  arc  shown  graphically  in  Figure  9  on  page  24. 

Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 
(3.62) 

Checking  for  stability  of  our  example  HR  notch  filter;  we  have  from  ©  of  (3.62) 

i  +  *,  +  (*J  4-  b})  >  0 

1  +  (0.32993)  +  [(0.32993)1  +  ( -0.9S01)]  >  0 

0.4586S3S  >  0 
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Figure  8.  Condition  for  Stability  of  the  Augmented  Filter:  p  =  1  [Ref.  1:  p.  2S4] 
from  ©  of  (3.62) 

l  -  6,  +  {b}  +  52)  >  0 

1  -  (0.32993)  4-  [(0.32993)2  +  ( -0.9801)3  <  0 

-0.2011762  <0 

therefore  the  filler  is  unstable,  from  ©  of  (3.62) 

I  b*  +  I  <1 

l(0.32993)J  +  (  -0.9801)1  <1 
0.S7 12462  <  l 

Note  that  ©  is  violated  and  here  the  realization  will  be  unstable. 

4.  Condition  for  Stability  of  the  Augmented  UR  Notch  Filter:  p  «  3 

The  augmented  UR  notch  filter,  p»  3,  is  described  by  equation  (3.55)  and  to  test 
for  stability,  wc  analyze  the  following  characteristic  polynomial  resulting  from  it: 

ffc)  -i3  +  6,zJ  +  (&;  +  bfc  +  (bf  4-  25, 6j)  (3.63) 
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Figure  9.  Condition  for  Stability  of  the  Augmented  Filter:  p  «  2  (Ref.  i:  p.  285) 


The  conditions  for  stability  arc: 

1  +  6,+(6?  +  6,)  +  (6?  +  26,6,)>0...© 

1  -  h,  +  (6?  +  6,)  -  (6?  +  2 6,6,)  >  0 ...  © 

|6?  + 26,6,1  <1  ...© 

I  (6, 3  +  2 6,6J2  -  1 1  >  1 6,(6?  +  26,6,)  -  (6?  +  6,)  I  ...  © 

These  conditions  arc  shown  graphically  in  Figure  10  on  page  25. 

Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 
(3.64) 

Checking  for  stability  of  our  example  HR  notch  filter;  we  have  from  ©  of  (3.64) 

1  +  6,  +  (6?  +  6,)  4*  (6?  +  26,6,)  >  0 

1  +0.32993  +  [(0.32993)2  +  ( -0.9301)3  +  [(0.32993)3  +  2(0.32993)(  -0.9801)]  >  0 

-0. 1 521 3085  <0 
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Figure  10.  Condition  for  Stability  of  the  Augmented  Filter:  p  -  3  (Ref.  1:  p.  2S6] 
therefore  the  filter  is  unstable,  from  ©  of  (3.64) 

1  —  +  (if  +  hj)  —  (if  +  2i]ij)  >  0 

1  -  0.32993  +  [(0.32993)2  +  ( -0.9801)]  -  [(0.32993)3  +  2(0.32993)(  -0.9801)]  >  0 

0.40963S46  >  0 


from  ©  of  (3.64) 


I  if  +  2i)bjjl  <  1 

I (0.32993)3  +  2(0.32993)(  -0.9S01)!  <  1 

0.61081465  <  l 

from  ©  of  (3.64) 


i(if  +  2i,i2)J-l{  >  |it(i?  -f  2i,ij)  -  (if  +  ij)| 

I  {(0.32993)3  +  2(0.32993)(  -0.9801)}1  -  1 1  < 

I  (0.32993){(0.32993)5  +  2(0.32993)(  -0.9801))  -  {{0. 32993)*  +  ( -0.9801)}  | 

0,62690546 <  0.66972012 
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therefore  the  filter  is  unstable. 

Note  that  (D  and  ®  are  violated  and  therefore  the  realization  will  be  unstable. 

G,  FREQUENCY  RESPONSE 

From  the  transfer  functions  (3.22),  (3.33),  (3.44),  and  (3.55),  we  use  the  program  in 
Appendix  A  to  calculate  the  magnitude  and  phase  of  the  frequency  response  of  the  se¬ 
veral  realizations.  We  plot  the  magnitude  of  these  responses  for  the  original  filter  and  for 
the  pipeline  augmented  that  all  of  the  frequency  responses  are  the  same  from  those  fig¬ 
ures. 

1.  Second  Order  IIR  Notch  Filter:  p  =  0 

The  magnitude  response  is  shown  in  Figure  11  on  page  27. 

2.  The  Augmented  IIR  Notch  Filter:  p  =  1 

The  magnitude  response  is  shown  in  Figure  12  on  page  28. 

3.  The  Augmented  IIR  Notch  Filter:  p  =  2 

The  magnitude  response  is  shown  in  Figure  13  on  page  29. 

4.  The  Augmented  IIR  Notch  Filter:  p  =  3 

The  magnitude  response  is  shown  in  Figure  14  on  page  30. 

We  know  from  section  F  that  the  example  augmented  IIR  notch  filter  for  p  » 
2,  and  3  are  unstable  but  we  cannot  see  that  from  these  magnitude  responses.  We  have 
to  find  another  way  to  show  which  one  is  stable  or  unstable.  This  will  be  shown  in  sec¬ 
tion  H. 

H.  UNIT  SAMPLE  RESPONSE 

From  the  difference  equations  (3.23),  (3.34),  (3.45),  and  (3.56),  wc  use  the  program 
in  Appendix  B  to  calculate  the  output  for  a  unit  sample  input  to  each  of  the  pipeline 
realizations.  We  plot  the  outputs  in  Figures  15,  16,  17  and  18.  We  can  sec  that  the  out¬ 
puts  for  Figure  15  (p»0)  and  Figure  16  (p®  1)  are  stable,  however  the  outputs  for  Fig¬ 
ure  17  (p«2)  and  Figure  18  (p“3)  arc  unstable  as  predicted  by  the  stability  testing  in 
section  F. 

L  Second  Order  IIR  Notch  Filter:  p  -  0 

The  output  is  shown  in  Figure  15  on  page  31. 

2.  The  Augmented  IIR  Notch  Filter:  p  °  1 

The  output  is  shown  in  Figure  16  on  page  32. 

3.  The  Augmented  HR  Notch  Filter:  p  =  2 

The  output  is  shown  in  Figure  17  on  page  33. 
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Figure  12.  Magnitude  Response  of  the  Augmented  HR  Notch  Filter:  p  «•  1 
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Figure  13.  Magnitude  Response  of  the  Augmented  HR  Notch  Filter:  p  -  2 
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Figure  14.  Magnitude  Respouse  of  the  Augmented  UR  Notch  Filter,  p  ■  3 


X  10** 


SAMPLE  «  (ns) 


4.  The  Augmented  HR  Notch  Filter:  p  =  3 
The  output  is  shown  in  Figure  18  on  page  35. 

I.  PIPELINE  IMPLEMENTATION 

All  of  the  IIR  filters  in  sections  B,  C,  D,  and  E  can  be  developed  into  high-speed  HR 
filters  by  applying  pipeline  techniques.  We  use  Ka-stage  pipeline  add  units  and  Km-stage 
pipeline  multipliers  to  increase  the  speed  of  the  filters.  The  Km  stages  are  involved  in  the 
coefficient  multiplies  and  the  Ka  stages  are  involved  in  the  adds.  The  increase  in  order 
of  the  filter  is  P  =  2Ka  +  Km  -  2  (Ref.  1:  p.  280]. 

1.  Second  Order  IIR  Notch  Filter:  p  =  0 

To  achieve  a  p  =  0  realization,  we  cannot  use  pipeline  modules  with  longer  de¬ 
lays  than  Ka  -  1  and  Km  =  0. 


P=2Ka+  Km-2 
—  2(1)  +  0-2 
=  0 

We  have  a  realization  of  the  original  IIR  notch  filter  but  with  an  additional  3 
delays.  The  equation  is 

r“3y  =  r3[<vY+  «,*■'*  +  a2:~2X+  b^lY+  ^z"2y].  (3.65) 

This  realization  is  shown  in  Figure  19  on  page  36. 

2.  The  Augmented  1IR  Notch  Filter:  p  =  i 

To  achieve  a  p=  1  realization,  we  cannot  use  pipeline  modules  with  longer  de¬ 
lays  than  Ka  «=*  l  and  Km  *=  1. 


P  «  2Ka  +  Km  -  2 
=  2(1)  -4-1  —  2 

ea  ] 


We  have  a  realization  of  the  augmented  IIR  notch  filter,  p  ■=  1,  but  with  an 
additional  4  delays.  The  equation  is 

Y  -  ^[ooA'  +  axz~'X  +  a^X  +  a3z~}X  +  Y  +  T].  (3.66) 


This  realization  is  shown  m  Figure  20  on  page  37. 
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Figure  19.  Second  Order  I1R  Notch  Filter,  p  «  0  with  Ka=  1,  Km  =  0 

3.  The  Augmented  IIR  Notch  Filter:  p  *»  2 

To  achieve  a  p  =  2  realization,  we  cannot  use  pipeline  modules  with  longer  de¬ 
lays  than  Ka  *=  1  and  Km  =  2. 


P=2Ka+  Am -2 
“  2(1)  +  2  -  2 

“  * 

Wc  have  a  realization  of  the  augmented  HR  notch  filter,  p  «  2,  but  with  an 
additional  5  delays.  The  equation  is 

Y -  +  a2r"lv  +  a}:~3X  +  Ck^x  +  V'J  Y  +  K].  (3.67) 

This  realization  is  shown  in  Figure  21  on  page  38. 

4.  The  Augmented  HR  Notch  Filter  p  -  3 

To  achieve  a  p  -  3  realization,  we  cannot  use  pipeline  modules  with  longer  de¬ 
lays  than  Ka  •*  1  and  Km  «  3. 
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P  —  2Ka  +  Km  —  2 
=  2(1)  +  3-2 
=  3 


We  have  a  realization  of  the  augmented  IIR  notch  filter,  p  =  3,  but  with  an 
additional  6  delays.  The  equation  is 


OqX  ■+■  fljZ  X  4-  Q-jZ  X  +  Q$Z  X 
+  af*X  +  a$z~5X  +  V"4  Y  +  bsz~s  Y 


(3.68) 


This  realization  is  shown  in  Figure  22  on  page  40. 


Table  4.  LOQMIS/SINHA  PIPELINE  RECURSIVE  DESIGN 


p-order 

Ka-add  units 

Km-multipliers 

D  (output  delayed) 

Figure 

0 

1 

0 

3  • 

19 

1 

1 

1 

4 

20 

■> 

1 

2 

5 

21 

3 

1 

3 

6 
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J.  SUMMARY 

This  chapter  investigated  a  method  for  applying  pipeline  techniques  to  the  design 
of  a  high-speed  UR  no!  *h  (liter.  Using  the  Loomis  Sinha  method,  HR  notch  filters  op¬ 
erating  at  rates  hitherto  impossible  can  be  designed.  The  general  structure  of  the  HR 
notch  filter  and  the  method  of  calculating  the  multiplier  coefficients  have  been  presented. 
The  stability  of  the  resulting  realization  has  been  investigated  and  a  technique  for  as¬ 
certaining  the  stability  of  the  realization  has  been  presented.  The  delay  and  complexity 
of  each  realization  are  summarized  in  Table  4. 
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IV.  GNANASEKARAN'S  STATE  SPACE  REALIZATIONS 


IIR  filter  realizations  are  possible  and  most  of  them  can  be  presented  within  the 
general  time  domain  representation  of  the  system  under  consideration.  In  continuous 
time  it  was  convenient  to  represent  a  linear  constant  coefficient  differential  equation  by 
a  system  of  first-order  differential  equations  for  analysis  purposes.  It  is,  therefore,  not 
surprising  that  it  is  useful  to  represent  a  linear  constant  coefficient  difference  equation 
by  a  system  of  first-order  linear  constant  coefficient  difference  equations.  Such  a  for¬ 
mulation  results  in  a  special  form  of  representation  of  the  difference  equation  called  a 
state  variable  representation, 

A.  STATE  SPACE  REPRESENTATION  OF  A  NTH  ORDER  RECURSIVE 
DIGITAL  FILTER 

We  are  used  to  working  with  a  difference  equation  of  the  form 
y(n)  =>  b^n  -  1)  +  -  +  b^yin  -N)  +  a^x^n)  +  Gjx(/t  -  I)  +  •••  +  aLx(n  -  \ f)  (4.1) 
or  in  the  notation 


S  M 

=  +  (4.2) 

km  l  A«0 

that  is  called  the  Nth-order  form.  An  alternative  is  to  describe  an  IT1  system  with  N 
first-order  difference  equations  and  an  output  equation.  This  alternative  form  for  de¬ 
scribing  a  linear  time-invariant  system  is  known  as  the  state-space  description  and  can 
be  written  us  a  matrix  equation: 

v(«)  e  Av(/t  -  1)  H*  Bx(n)  (state  equation)  (4.3) 

with  the  system  output  given  by 

yin)  =  Crv(«)  +  dx(n)  ( output  equation).  (4.4) 

In  Eqs.  (4.3)  and  (4.4),  v(«)  is  a  vector  with  N  rows  and  l  column  (an  <V  x  I  vector  or 
column  matrix)  of  newly  defined  variables  called  the  states,  x(n)  is  the  system  input,  and 
yin)  its  output.  In  the  matrix  state  equation  of  Eq.  (4.3)  the  matrix  A  must  be  an 
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NxN matrix  and  the  B  matrix  an  N x  1  vector.  The  scalar  output  >•(«)  is  defined  by  Eq. 
(4.4)  where  C  is  a  N  x  1  column  matrix  and  d  is  simply  a  real  number  that  multiplies  (or 
scales)  the  input  x(n).  The  elements  in  A,  B,  C,  and  d  are  constants  that  depend  only 
upon  the  system  parameters. 

B.  WRITING  THE  STATE  SPACE  EQUATIONS 

To  illustrate  this  approach  to  modeling  of  LTI  systems  we  return  to  a  version  of  the 
Nth-order  difference  equation  described  as  a  second  order  HR  filter,  namely 

y(n)  =*  b^n  —  I)  +  b-y{n  —  2)  4-  UqJ :(n)  +  axx(n  —  1)  +  c2x(n  —  2)  f4.5) 


Y**bxz~xY  +  b2z~2  Y+OQX+  axz~xX  +  a2z~2X. 

(4.6) 

We  define  the  state  \arible  as:  Vx,  V2.  Then,  we  can  establish  the  following  re- 
lationships: 

Vx  *  J'=»  bxz~x  Y  +  b2z~2  Y+OqX+  axz~xX  +  a2z~2X 

(4.7) 

Therefore  the  state  equation  is: 

t',  -  b,;~‘  +  4;-_1  Kj  +  oqX  +  o,r'A  +  a,z'\f 

(4.3) 

We  express  it  in  the  matrix  state  equation: 

v(/») «  A  v(n  -  1)  4-  B.v(m) 

(4.9) 

V  «  Az“lV  +  BA' 

f4.ro) 

[3-ir 

.  f4.n> 

That  means: 


(4.13) 


B 


Qq  4*  Q\Z 

<> 


The  output  equation  will  be: 

><«)  -  Crv(/j)  +  dx(n) 

rv  1 

v-cc.,  .ill  1t*r 


y=c,K, +C2  V2  4-  dX. 

C.  MODIFIED  STATE  SPACE  EQUATION 

From  the  state  equation 

v(«)  =  Av(«  -1)4-  B.v(n), 


(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 


the  computation  of  v(n)  requires  the  immediate  past  state  vector  v(n*l).  This  recursive 
nature  severely  limits  the  speed  that  can  be  achieved  by  direct  implementation  of  the 
state  equation.  The  speed  limits  can  be  overcome  by  modifying  this  state  equation  so 
that  v(n)  depends  not  on  the  immediate  past  vector,  v{n  -  1),  but  on  vpt  -  L)  for  some 
L  >  1.  This  can  be  accomplished  as  follows: 

For  L  «  1  the  state  equation  is 


v(«)  ■»  Av(«  -1)4*  Bx{h) 

(4.19) 

V  «  A^V  +  BX. 

(4.20) 

Delaying  equation  (4.19)  by  1,  2,  and  3  yields 

\{n  -  1)  =  Av(h  -  2)  -f  B4«  -  l) 

(4.21) 

\\n  -  2)  =  Av(«  -  3)  4-  lk(n  —  2) 

(4.22) 

v(«  -  3)  =>  Av(«  -  4)  4*  Hf(«  -  3). 

(4.23) 

For  1  =  2  putting  (4.21)  into  (4.22),  the  state  equation,  results  is 
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v(«)  =  A[Av(m  -  2)  +  Bjt(«  -  1)]  +  Bx(>i)  (-4-24) 

v(/j)  =  A2v(/t  -  2)  +  AEx(n  -  1)  +  B*(/0  (4.25) 

V  =  A2z~2V  +  ABr'1^  +  EX.  (4.26) 

For  L  =  3  putting  (4.22)  into  (4.25),  the  state  equation,  results  is 

y(n)  =  A2[Av(rt  -3)  +  Bx(«  -  2)]  +  ABx(it  -  1)  +  Bx(n)  (4.27) 

v(n) »  A3v(«  -  3)  +  A2Baf(w  -  2)  +  AB*(*  - 1)  +  B x(n)  (4.28) 

V  =  A3z“3V  +  A2Bz~2,V  +  ABz"1  A  +  EX.  (4.29) 

For  L  =»  4  putting  (4.23)  into  (4.28),  the  state  equation,  results  is 


v(«)  -  A3[  Av(«  -  4)  +  B x(n  -  3)3  +  A2B*(«  -2)  +  ABx(«  -  1)  +  B*(it)  (4.30) 
v(/i) «  A4v(«  -  4)  +  A3Ba(/i  ~  3)  +  A2Ba(m  -  2)  +  AB.r(/i  -  1)  +  B.v(/i)  (4.31) 

V  .  AJz~Sf  4-  A3Bz~3A'  4-  A!Bz“2A'  +  ABz“'A'  4-  EX.  (4.32) 

Equation  (4.32)  can  be  written  in  the  following  form  for  any  L>  1: 


A  L:~LV  4- 


/.-i 

I 

M) 


AfB s~*X. 


(4.33) 


One  of  the  benefits  of  the  state  space  realization  is  that  if  the  original  filter  (L*  1) 
is  stable,  ali  higher  order  realizations  are  al;o  stable,  and  furthermore,  A1  has  improved 
error  properties  over  A.  This  advantage  of  the  state  space  method  is  discussed  in  the 
paper  by  Lu,  Lee,  and  Mcsscrschmitt  (Ref.  4). 

D.  SECOND  ORDER  1IR  FILTER  STATE  SPACE  REALIZATIONS 

From  the  state  space  equation  of  the  second  order  UR  filter  in  section  B  and  the 
modified  stale  space  equation  in  section  C,  we  can  develop  techniques  for  high  speed 
implementation  of  second  order  UR  filters  for  several  values  of  L.  litis  section  shows 
how  to  apply  these  concepts  for  L  m  1,2,  3,  and  4: 
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1.  L  =  1 


Beginning  with  equation  (4.20), 

V  =  Az"1V  +  BX 


and  substituting  for  V,  A,  and  B  from  section  B  yields: 


(4.34) 


(4.35) 


Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  23, 
we  note  that  the  outpiu  of  the  multiplier  d  is  dz~7X.  This  results  is  a  revised  output 
equation  of 

Y  =  c,  Vx  +  c2  V2  +  dz'2X.  (4.36) 

Substituting  the  state  equation  into  the  output  equation,  results  in 

r«(Vi+c2)rlF, 

+  b2cxz  V2 

+  (4.37) 

A 

+  (a2c,  +  d)z~2  X 

Y~(lY\  +  c2)z-'Y 
+  b2cxz~2Y 

+  floM'  (4.3S) 

+  axc{2~'X 
+  (fitjCj  +  d)z~2X. 

To  demonstrate  equality  with  the  Y  of  the  original  MR  filter  as  expressed  in  equation 
(3.23)  for  p  “  0,  let  c,  =  1,  Cj  =  0,  and  d=  0  .  Then  the  output  equation  is: 

Y  =  fc,z~l  Y  +  ^r2  Y  +  OqX  +  axz~'X  +  a2z~2X.  (4.39) 


By  substituting: 
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(4.40) 


62  =  bx  2 
%  —  an 
a\  ~  an 

a2  ~  a13 

we  can  write  this  equation  as: 

Y  -  6,  j2~ '  r  4-  i>i22"2  r  +  anJf  +  a,22-1*  +  a132~2X 

This  realization  is  shown  in  Figure  23  on  page  47. 

2.  L  =  2 

Beginning  with  the  state  equatic  \  (4.26), 

V  =  A22"2V  +  AB2',A'+BJ 
and  substituting  for  V,  A,  and  B  from  section  B  yields: 


- j 

_ i 

"*i 

’i  V 

,~2 

L<  J 

J 

o  JLj 

o 

c 

1 _ 

‘6,  Vf 

.1  OJ 


i  i  -n 

a0  4-  a,2  4-  a22 


A' 


<4)  4-  rt,2~'  4-  fl22~2 


A' 


£3 

4?  +  4, 

V>2 

A 

A 

/l 

*3  . 

v2 

<%*,  +  1  +  1 

_2  2  A 

4-  a,s  4-  o2?  j 

T 


'0,  +  Oij"1  +os.."! 

+  o  |A 


Oq  4*  (<2o^i  4*  rtO2”1  +  (<*A  +  Oj);-2  4- 


-3 


Uqa  4*  4“  $2* 


.V. 


(4.41) 


(4.42) 


(4.43) 


(4.44) 
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Figure  23.  Second  Order  HR  Filter.  L  -  1 
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Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  24, 
we  note  that  the  output  of  the  multiplier  d  is  dz~lX.  This  results  is  a  revised  output 
equation  of 

Y  =  Cx  Vx  +  c2  V2  +  dz~3X.  (4.46) 


Substituting  the  state  equation  into  the  output  equation,  results  in 

4"  c2q  C\X 

+  C(«0*1  +  a0ci  +  aoC2]2-1  J 
+  (axbxcx  +  a2cx  4-  axc2)z~2X 
+  (a2bxcx  +  a2c2  +  d)z~3X 

Y=[(bt  +  b2)cx+bxc2y2Y 

+  4*  Y 

4*  OoC,/V 

+  CK^I  +  <*\)C\  4-  00^]^“’^ 

4-  (<j,^,c,  4-  ajCj  4-  axc2)z~2X 
4-  (fl26jC1  4-  fljCj  4*  d)z  3X. 


(4.47) 


(4.4S) 


To  demonstrate  equality  with  the  Y  of  the  augmented  UR  filter  as  expressed  in  equation 
(3.34)  for  p  “  1,  let  c,  »  1,  c}  «=  0,  and  </<=  0.  Then  the  output  equation  is: 

>•«  (i>f + y + w1  y + «oX + m + 

(Hj6;  +  «j);  3 A  +  ;  3 A"- 


By  substituting: 

* 

+  bj  =  bx  i 

bih  “  6,j 

Oobx  4-  fl)  “  0)2 
axbx  4-  a2  =  a(j 
a26,  «■  a,  4 


wc  can  writs  this  equation  as: 


(4.50) 
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(4.51) 


7=  bnz~2Y  +  bnz~3Y+  auX+  anz~lX 
+  ^132  X  4*  {*142  X. 

This  is  the  same  as  the  equation  (3.34)  for  p=  1  found  on  page  15.  This  realization  is 
shown  in  Figure  24  on  page  50. 

3.  L  -  3 

Beginning  with  the  state  equation  (4.29), 

V  =  A32~3V  +  A2B2~2*  +  ABz~{X  +  BX  (4.52) 

and  substituting  for  V,  A,  and  B  from  section  B  yields: 


Vlfa  +  fii*"1  +<V'21  -iv 

.  oJL  0  y  x 

«o  4*  a,z"‘  4-  «2z~2] 

°  J 


(4.53) 


+ 2^)  _,[>,] 

i.,fc  J*  W 

■+•  bj)  +  +  fe,)z  1  +  +  bj)i  1 

4*  i 

a^i  4- a|i>|2  4- 036,2 

f  a^,  +  0,6,2"'  4-  a26,2_2l 
4-  ]  _i  2  A 

[  Oq  4*  a,*  4-  a3: 

L  0  f 


(4.54) 
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Figure  24.  The  Augmented  HR  Fitter:  L  *=  2 


SO 


bJht  +  Tbd 
b]  +  b2 


b2(bt  +  b2) 
b\b2 


-3 


lA 


+ 


Oq  +  (a^by  +  +  (a0by  +  +  a\b\  +  a2)2-2 

+  (a^f  +  a]i2  +  o2bx)z~i  +  (a2bx  +  a2b2)z~* 

%z~x  +  (oq^i  +  a\)z~2  +  (a\bx  +  a2)z~3  +  a^z-4 


X. 


(4.55) 


Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  25, 
we  note  that  output  of  the  multiplier  d  is  dz~*X.  This  result  is  a  revised  output  equation 
of 


Y  -  Ci  K,  +  c2  V2  +  dz~*X.  (4.56) 

Substituting  the  state  equation  into  the  output  equation,  results  in 

Y  «  Of  +  Ibybjcy  +  (b}  +  62)c2]2~3F, 

+  Ofo  +  b2)Cy  -r  6t62c2]2_3  V2 
+  OoC,A' 

+  l(%bx  +  a,)c,  +  «oC2J:_iA’  (4.57) 

+  [(oq/jj2  +  +  «\b\  +  +  foby  +  a,)^];-2^ 

+  [(a,6f  +  +  a2by)cy  +  (ayby  +  U2)C2>“3A' 

+  +  «2^lc2  + 


v-0,34-2Mih  +  (/h3  +  ^h>‘3^ 

+  [(fcf£»2  +  622)c,  +  bybfalz^Y 

+  Upt'jA 

+  C(^|  +  (4.5S) 

+  [(tyfrf  +  +  axby  +  a2)c,  +  (a^,  +  UjJcjJj  *JA( 

+  +  aib\)t\  + 

+  (](a2i>f  4-  a2bi)cy  +  o2byC2  +  dAr. 

To  demonstrate  equality  with  the  Y  of  the  augmented  UR  filter  as  expressed  in  equation 
(3.45)  for  p  =  2,  let  c,  =  1,  c3  =  0,  and  d  =  0.  Then  the  output  equation  is: 
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y«  (bl  4-  2  blb2)z~*Y+  {b\b2  +  b})z~AY  +  acX+  {%bx  +  ax)z~xX 

+  £ <^[b\  4-  i>2)  +  &\b}  4*  o-^z  2Af  4-  4-  bj)  +  d2b^ X  (4.59) 

+  a2{b]  4-  b2)z~*X. 


By  substituting: 

b\  +  2b^b2  ~  b i  j 

b\b2~Y  b%  —  b\ 2 
ao  =  «u 
a§b\  +  Qj  =  d\2 

Oq  ,(&?  +  b2)  4-  axbx  4-  a2  =  n13 

a\{b\  +  b2)  +  <*2^1  ~  a\i 
a2{b*  +  b2)  m  a15 

we  can  write  this  equation  as: 

Y  m  b,  ,z~3  Y  +  &j2z~*  y+anA'+  a, 2z~lX 
+  X  4"  X  +  <J;s2  Af* 

Thi*  is  the  same  as  the  equation  (3.45)  for  ptB2  found  on  page  17.  This  realization  is 
shown  in  Figure  25  on  page  53. 

4.  L  ™  4 

Beginning  with  equation  (4.32), 

V  -  A4z“*V  +  A3BT3A'  ■<*  AJBT2A'  4-  ABz~'X  4-  BA'  (4.62) 

and  substituting  for  V,  A,  and  B  from  section  B  yields: 


(4.60) 


(4.61) 


Figure  23.  The  Augmented  UR  Filter  L  *  3 
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Vi 
V 2 


bx{b2x+2b2)  b2(b2  +  b2) 


bx  +b2 


b,b 


\v2 


A  b2 
1  0 


Vi 


4- 


b\(b\  4-  2b2)  b2{b\  +  52)"lr<2Q  4"  a\z  1  +  a2z  2 


b{  4-  b2 


b\b2 


JL 


0 


z~2X 


K 

[:■ 

[ 


+  b2 

V 

0 


bxb2 

b2 


—1  — 2T 

Oq  4-  axz  4-  a2z 


-i 


%  +  axz  4-  a2z 
0 


0 

-2 


h 


z'2* 


Oq  +  axz~ 1  4-  a2z~ 2_l 


(4.63) 


bx  4-  3 b\b2  +  b l  bxb2{b\  4-  2 b2) 


iA 


lAJ  [b:(bt  +  2b2)  b2{b{  +  b2) 

aabx(bx  4-  2 b2)  4-  axbx(bx  4-  2 b2)z~l  4*  a2bx(bf  +  2 b2)z' 
flo(6,  4-  b 2)  4"  0\{bx  4-  i2)z  4"  ct2{bx  4"  &2)z 


Oq {bx  4"  b2)  4"  4"  b2)z  4-  o2(bx  4"  b2)z 

Oo^t  4-  «,6,z*'1  4-  a2bxz~ 2 


r2* 


I  «o£>t  4-  axbxz~'  4-  a2&,r2 1 
L  «o  4-  4-  flj2“2  !' 

[flo  4- a,r'  4- a2z’21 

L  «  JA 


z-3* 


(4.64) 


bx  4-  4-  hi  +  2/>2)l  JV 

6,(6?  4*  26j)  ^4-6,)  J2  |/2 


+ 


4o  4*  (og^t  +  01)2“'  4*  (^C  4*  a^b}  4-  4-  a^z""2 

4*  (tfg62  4*  2^6)62  4-  4-  4*  3 

4*  (ajfrj  4*  2oxbxb2  4-  o2b2)z  4"  (flitj  4" 

ao2“'  4-  (flof>,  4-  U|)2-2  4-  (a^,2  4-  Oq^  4-  a, 6,  4-  a2)*“3 

4*  (<lji>j  4“  4"  fl2&j)z  4“  (djiij  4*  <^^>2)2 


A'. 


(4.65) 
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Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  26, 
we  note  that  the  output  of  the  multiplier  d  is  dz~sX.  This  results  is  a  revised  output 
equation  of 


Y—  cxVx  +  c2V2  +  dz~sX. 
Substituting  the  state  equation  into  the  output  equation,  results  in 


(4.66) 


Y  =  |>f  +  3  b\b2  +  b\)cx  +  {b\  +  2^2)c2]z‘V, 

+  Of*2  +  2^2  )C,  +  {b\b2  +  bl)c1\z~AV1 
+  %cxX 

+  [><A  +  a\)c\  +  ^C2]z-1  X 
+  [(a0fef  +  %b2  4-  axbx  +  «2)cj  4-  4-  <Jj)c2]z“2A: 

+  {_(%bx  4*  2  %bxb2  +  axb]  4-  at62  4-  a2bx)cx  +  (oq^2  4-  %b2  4-  axbx  4-  a2)c2]z-3^ 

+  Oi&f  +  2*,^  4-  a2b*  4-  a2b2)cx  +  (atif  +  axb2  4*  a261)c2]z“4A' 

+  £(<j2^j  4“  ‘2a2bxb2)cx  4*  {o2bx  +  fl2^>2)c2  +  d Jz  A 

Y=  Of  +  3  b}b2  +  &2)c,  4-  {bx  +  2bxb2)c2']z~AY 
4-  [(^>f/*2  +  2  bxbl)cx  +  (6fo  +  i>?)c2]r“5y 

*4“  (fyCjA 

+  [(V’l  +  «iki  +  X 

+  Oo/>f  +  «o^2  +  Q\b\+  02)c,  4-  (a^,  +  a,)c2]z“2A'  ^4'68^ 

+  Oo^f  +  2oQbxb2  +  axb f  4-  axb2  4-  a2bx)cx  4-  (<?0i>,2  4-  4-  axbx  4-  «j)c2]z'3,V 

4-  [(a,(>3  4-  2axbxb2  4-  a2bx  4-  a2b2)cx  4-  {axbx  4-  axb2  + 

4*  [(tf2^3  4*  2a2fc,&2)C|  4-  {a2bx  4-  a2b2)c2  4-  <f]z“sA*. 

To  demonstrate  equality  with  the  Y  of  the  augmented  UR  filter  as  expressed  in  equation 
(3.56)  for  p  «  3,  let  c, »  1,  c2 »  0,  and  d<*  0.  Then  the  output  equation  is: 


r«(6f  4-  3 bfa  4-  6jJ)z_4K4*  (frfo  4-  26,fr2V5>'4-  q>A'  +  («<A  +  ^“‘AT 
4-  0(&f  4*  (>2)  4-  axbx  4-  a2]z”2A’' 

4“  f^^3  4*  2bxbj)  +  <i|(62  4*  bj)  4*  Jz  ^ 

4*  £ fl|(63  4*  2 bxbj)  4*  o2(bx  4*  4*  o2{b\  4-  2bxb2)z  *X. 


(4.69) 


Bv  substituting: 
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(4.70) 


bt  +  ^b2  +  bl  =  bu 
b\b2  +  2b^b2  =  b]2 

flo  =  an 
^ji?i  +  tfj  =  d\2 

%{b\  +  b2)  +  Q\b\  +  a2  =  a13 
ijg  (b}  +  2b\b2j  +  +  ij)  +  a2b\  —  Gq 4 

•+■  2  b^b^)  +  q2(&]  +  b2)  —  Gjj 
<*2(^1  +  2b\b2)  =  Qjg 


we  can  write  this  equation  as: 

y  =  &HZ-4  y  +  *>122_S  Y+anX  +  a122-1A: 

+  a,32"JA'+al42"3A'  (4.71) 

+  a^z^X  +  a^z~sX. 


This  is  the  same  as  the  equation  (3.56)  for  p  =  3  on  page  20.  This  realization  is  shown 
in  Figure  26  on  page  57. 

From  equation  (4.33)  we  can  see  that  the  feedback  portion  of  the  realization 
must  have  a  delay  of  L  stages  in  any  loop.  The  general  realization  of  that  portion  is 
shown  in  Figure  27  on  page  58. 

E.  PIPELINED  IMPLEMENTATION 

All  of  the  UR  filters  in  section  D  can  be  developed  into  high-speed  UR  filters  by 
applying  pipeline  techniques.  We  use  Ka-stage  pipeline  add  units  and  Km-stage  pipeline 
multipliers  to  increase  the  speed  of  the  filters.  The  Km  stages  are  involved  in  the  cocfli- 
cicnt  multiplies  and  the  Ka  stages  are  involved  in  the  adds.  The  net  result  is  a  (Km  *f 
Ka)-stagc  pipeline  multiply-add  unit.  The  equations  for  the  filter  outputs  for  different 
L's  arc  given  below: 

I.  L  -  I 

The  equation  is 


-<SA‘a+3£«)  v  .-tSA’a+lAm) 
2  ■  4=3  * 


■J  bnz~«*0+*m)Y+  bl2:~ilk't'km+')Y 
•+•  rijjA  +  flj2z  X  +  a)3i  A 


(4.72) 


This  rcalizn*ic;ns  shown  in  Figure  28  on  page  59. 
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Figure  26.  The  Augmented  HR  Filter.  L  -  4 
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Figure  27.  Heart  of  the  Pipeline  Recursive  Filter  Realization 


?  In? 

«*  L  w  a* 

The  equation  is 

a  /  **  . 


h,  ,2"«*a+A‘w)  y  4-  b^-®**'***"  Y 


4  a,  j.V  4  Ujjs  *,V  4  o,j2  W  4  0(45  *A 


-2, 


-1 


This  realization  is  shown  in  Figure  20  on  page  60. 
3.  L  -  3 

The  equation  is 


bui-V*°+km)Y+  u^^^Y 
4  <tj  ,AT  4  0,j£  A  4  0,j2  A 

4  0,42^4  0,52^ 


This  realization  is  shown  in  Figure  30  on  page  61- 
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(473) 


(4.74) 


L 
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Figure  29,  The  Augmented  HR  Filter  L  «  2  with  Ka,  Km 
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Figure  30.  The  Augmented  UR  Filter.  L  ■  3  with  Ka,  Km 
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4.  L  =  4 

The  equation  is 

(4.75) 

This  realization  is  shown  in  Figure  31  on  page  63. 

From  analysis  of  the  feed  back  loop  portion  of  the  state  space  realization,  as 
shown  in  Figure  27  and  equation  (4.33),  we  have  seen  that  the  delay  around  any  loop 
must  be  L.  A  general  realization  for  the  feedback  portion  of  the  pipeline  case  is  shown 
in  Figure  32  on  page  64. 

Therefore  we  can  see  that 


-{1  \Ka+2Km)y  _  z~{\  lKa+2Km) 


-tanA-tanz  a  -t-  al3z  a 
+  Q\^Z  X  +  flu  2  X  +  jZ 


L  =  2Ka  +  Km.  (4.76) 

The  following  examples  show  the  analysis  of  the  cases  for  L  =»  1, 2,  3  and  4. 

5.  Example  for  L  =  1: 

a.  Ka  »  0,  Km  *=  1 

L  »  2  Ka  +  Km 
—  2(0)  +  1 
*=  1 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation  (4.76)  are 
Ka  «  0.  Km  <=*  !.  This  means  we  cannot  use  pipeline  modules  with  longer  delays  than 
Ka  «  0,  Km  «  1  . 

This  realizes  the  original  UR  filter  for  L  «  1,  but  with  an  additional  2  de¬ 
lays  and  the  equation  is 

f'Y  M  r2[f»lt2“*  y+  bn:~2Y  +  «,,*  +  +  (4.77) 

This  realization  is  shown  in  Figure  33  on  page  65. 

6.  Example  for  L  =  2: 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation 
(4.76)  arc  Ka  =  0,  Km  =  2  or  Ka  <*  1 ,  Km  =  0. 
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Figure  31.  The  Augmented  UR  Filter:  L  c  4  with  K»,  Km 
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Figure  32,  Heart  of  More  General  Pipeline  Filter  Realization 


a.  Ka  ®  0,  Km  =  2 


L  =  2A'a  +  Am 
-  2(0)  +  2 
CB  2 

For  these  values,  \vc  use  pipeline  modules  with  delays  Ka  =  0,  Am  «=  2  , 

This  realizes  the  augmented  HR  filter  for  L  ■»  2,  but  with  an  additional  4 
delays.  The  equation  is 

r4  Y  m  Y  +  A,  j2~J  r  +  a„  A  +  al22~’,v  +  a132“lA  +  auz‘JA'].  (4.78) 

This  realization  is  shown  in  Figure  34  on  page  66, 
b.  Ka  “  /,  Km  «  0 


L  -  2Aa  +  Am 
-2(l)  +  0 
«  2 

For  these  values,  we  us**  pipeline  modules  with  delays  Ka  «  1,  Am  *»  0 , 
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Figure  33.  Second  Order  UR  Filter;  L  -  1  with  Ki=0,Km=  1 
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Figure  34.  The  Augmented  !IR  Filter.  L  “  2  Hith  Ka**Q,  Km*2 
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This  realizes  the  augmented  IIR  filter  for  L  =  2,  but  with  an  additional  9 
delays.  The  equation  is 

z~9Y-z~9[buz~2Y+  bnz~zY+auX+  anz~xX+  auz~2X  +  auz~3J].  (4.79) 

This  realization  is  shown  in  Figure  35  on  page  68. 

7.  Example  for  L  =  3: 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation 
(4.76)  are  Ka  -  0,  Km  =  3  or  Ka  -  1,  Km  -  1. 
a.  Ka  -  0,  Km  -  3 


# 


L  =  2Ka  +  Km 
=  2(0)  +  3 
*=  3 


For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  0,  Km  =»  3  . 

This  realizes  the  augmented  IIR  filter  for  L  »  3,  but  with  an  additional  6 
delays.  The  equation  is 


b\\Z  2Y+b\2Z  "Y  +  +  fljjZ  'A' 

+  Gj 'X  *f  a14r  3  A  +  au?”4 A' 


(4.S0) 


This  realization  is  shown  in  Figure  36  on  page  69. 
b.  Ka  ■  I,  Km  m  / 


L  o  2Ka  +  Km 
«  2(1) +  J 

es  J 


For  these  values,  vve  use  pipeline  modutes  with  delays  Ka  =*  I,  Km  «  t  . 

This  realizes  the  augmented  UR  filter  for  L  =  3.  but  with  an  additional  12 
delays.  The  equation  is 


6,  ,j’j  Y  +  6I12~4F  +  auX  +  a,  jr'A 
+  0|]i  A  4*  A  4*  fljjS  X 


(4.S1) 


This  realization  is  shown  in  Figure  37  on  page  70. 
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Figure  35.  The  Augmented  11R  Filter;  L  ■  2  with  Ka=  1,  Km  =  0 
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figure  36.  The  Augmented  UR  Filter:  L  *  3  with  Ka®  0,  Km®  3 


69 


Figure  37.  The  Augmented  1IR  Filter.  L  “  3  nith  Ka=  i,  Km*  I 


8.  Example  for  L  “  4: 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation 
(4.76)  are  Ka  =  0,  Km  —  4  and  Ka  —  1,  Km  —  ?.  or  Ka  -  2,  Km  =  0. 
a.  Ka  —  0f  Km  =  4 


* 


L  =  2Ka  4  Km 
=  2(0} +  4 

XSS  4 


For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  0,  Km  =  4  . 

This  realizes  the  augmented  HR  filter  for  L  =»  4,  but  with  an  additional  8 
delays.  The  equation  is 


2-®y= 


-8 


6,  Y 4  £»,22*“5y  4-  a, , Ar  4  anz~'X 
4  a^X  4  aXiz-3X  t  4  a162~5A*  ‘ 


(4.82) 


This  realization  is  shown  in  Figure  38  on  page  72. 
6,  Ka  =  /,  Am  °  2 


L  =  2Aa  4  Am 
=  2(1)  4  2 

«$r  4 


For  these  values,  we  use  pipeline  modules  with  delays  A'a  «  1 .  Km  =  2  . 

This  realizes  the  augmented  UR  filter  for  L  =  4  but  with  an  additional  15 
delays.  The  equation  is 


6||2^)?4  5 Ft  *X 

4  sq^'A'  4  awr~'.V4  4  a)62~sA' 

M 


(4.S3) 


This  realization  is  shown  in  Figure  39  on  page  73. 

(.  Ka  =  2,  Km  =  0 

L  *=  2 An  4  Km 
=  2(2)40 
«  4 

For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  2,  Km  =  0 . 

This  realizes  the  augmented  UR  filter  for  L  =  4,  but  with  an  additional  22 
delays.  Tire  equation  is 
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V  >:\' 


figure  39,  The  AtigmenteJ  IIR  Filter:  L  -  4  with  Ka=  1,  Km  =2 
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•22 


(4.84) 


z~22Y  =  z 


b^Y+b^Y+a^+a^X 
+  d\^Z,  X  +  &\iZ  X  G]$Z  X  +  fljgZ  X 


This  realization  is  shown  in  Figure  40  on  page  75. 


Table  5.  STATE  SPACE  REALIZATION  SUMMARY 


L 

Ka-add  units 

Knvmultipliers 

D(  output  delayed) 

Figure 

1 

0 

1 

2 

31 

2 

0 

2 

4 

32 

2 

1 

0 

9 

33 

3 

0 

3 

6 

34 

3 

1 

1 

12 

35 

4 

0 

4 

8 

36 

4 

1 

2 

15 

37 

4 

2 

0 

22 

38 

F.  SUMMARY 

This  chapter  develops  an  alternate  method  for  high  speed  implementation  of  state 
space  digital  filters.  The  proposed  method  implements  \L  instead  of  A  as  seen  in 
equation  (4.33).  The  proposed  method  allows  arbitrary  A  in  the  implementation,  and  it 
can  take  full  advantage  of  the  optimum  state  space  filters  along  with  the  benefit  of  A4 
realization.  The  delay  and  complexity  of  each  realization  are  summarized  in  Table  5. 
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V.  CONCLUSIONS 


A.  SUMMARY 

This  thesis  describes  pipeline  realizations  of  the  HR  notch  filter.  The  transfer 
function  and  the  coefficients  that  are  appropiate  for  HR  notch  filter  have  been  pre¬ 
sented.  Realization  of  a  high  speed  HR  notch  filter  using  a  pipeline  technique  by  Loomis 
and  Sinha  has  been  presented.  An  alternate  technique  in  using  state  space  realization 
by  R.  Gnanasekaran  has  also  been  presented. 

B.  COMPARISON 

Table  4  in  Chapter  111  summarizes  the  Loomis, Sinha  pipeline  recursive  design  and 
Table  5  in  Chapter  IV  summarizes  Gnanasekaran's  state  space  realization.  It  was  shown 
in  Chapter  III  section  I  that  we  have  only  one  choice  of  Ka  and  Km  for  each  value  of 
p  and  in  Chapter  IV  section  E  that  we  have  1  choice  of  Ka  and  Km  for  L=  1,  2  choices 
of  Ka  and  Km  for  L=2  and  for  L=3,  and  3  choices  of  Ka  and  Km  for  L  =  4. 

In  the  Loomis/Sinha  pipeline  recursive  design  we  note  that  for  p=  2,  Ka=  1,  Km=  2, 
output  delay  =  5,  the  filter  is  unstable.  However,  in  Gnanaskaran's  state  space  realization 
for  L  =  4,  Ka=  1,  Km®  2,  output  delay®  15,  the  filter  is  stable.  Thus  we  see  that  in  cases 
where  Loomis/Sinha  realization  are  unstable,  a  state  space  realization  always  exists  and 
is  stable.  However,  it  will  require  more  hardware  and  have  a  greater  overall  delay.  All 
of  these  results  are  summarized  in  Table  6.  In  that  table,  the  column  labelled  D  refers 
to  the  overall  delay  and  the  columns  labelled  #Mlt  and  #Add  refer  to  the  number  of 
multipliers  and  number  of  adders  as  shown  in  the  appropriate  figure. 

Note  that  for  a  given  value  of  Ka  and  Km,  the  number  of  adders  and  multipliers 
required  for  the  state  space  realization  is  larger  than  that  required  for  the  scalar.  This 
disadvantage  is  offset  by  the  guaranteed  stability  of  the  state  space  realization. 
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Table  6.  SUMMARY  OF  REALIZATIONS 


oooooooooooooooooooonoooooonooooooooooonooooonno 


APPENDIX  A.  PROGRAM  AND  DATA  OF  FREQUENCY  RESPONSE 
A.  PROGRAM  FOR  FREQUENCY  RESPONSE  CALCULATION 

DIGFREQ.  FOR  VERSION:  2/03/88 


PURPOSE:  THIS  PROGRAM  COMPUTES  THE  FREQUENCY  RESPONSE  OF 

DISCRETE  SYSTEMS.  THE  PROGRAM  CONSISTS  OF  A  MAIN 
PROGRAM  THAT  CONTROLS  THE  INPUT/OUTPUT  AND  THE 
SUBROUTINES  dfresp  AND  coeff.  SUBROUTINE  dfresp 
COMPUTES  THE  FREQUENCY  RESPONSE  OF  EACH  SYSTEM. 
SUBROUTINE  coeff  ALLOWS  THE  USER  THE  OPTION  OF 
GENERATING  THE  FILTER  COEFFICIENTS  OF  THE  SYSTEMS 
TO  BE  ANALYZED  BY  WRITING  THE  APPROPRIATE  EQUATIONS. 

IF  THE  USER  ELECTS  TO  GENERATE  THE  COEFFICIENTS  BY 
USING  THE  SUBROUTINE  coeff,  THE  EQUATIONS  MUST  BE 
WRITTEN  INTO  THE  SUBROUTINE  USING  STANDARD  FORTRAN  77 
STATEMENTS.  THE  COEFFICIENTS  MUST  BE  STORED  IN  THE 
ARRAYS  b( )  AND  c()  WHICH  CORRESPOND  RESPECTIVELY  TO  THE 
NUMERATOR  AND  DENOMINATOR  TERMS  OF  THE  SYSTEM  EQUATION. 
THE  USER  CAN  SELECT  ONE  OF  TWO  OPERATING  MODES:  BATCH 
OR  TEST.  IN  BATCH  MODE  THE  AMOUNT  OF  INTERFACE  WITH 
THE  USER  IS  MINIMIZED  AND  IT  IS  ASSUMED  THAT  THE  INPUT 
DATA  HAS  BEEN  STORED  IN  THE  DEFAULT  FILE  ' DIGFREQ.  IN' . 

IN  TEST  MODE  THE  USER  IS  PROMPTED  FOR  THE  NAME  OF  THE 
INPUT  FILE  OR  HAS  THE  OPTION  TO  PERFORM  A  TRIAL  RUN  BY 
USING  THE  INPUT  DATA  STORED  IN  THE  FILE  ' DIGFREQ. TST'. 

IT  IS  RECOMMENDED  THAT  FIRST-TIME  USERS  SELECT  THE  TEST 
MODE  AND  MAKE  A  TRIAL  RUN  WITH  THE  PRESTORED  INPUT  DATA. 
THE  TEST  MODE  ECHOES  THE  INPUT  DATA  ONTO  THE  MONITOR  TO 
ALLOW  VERIFICATION  OF  ITS  ACCURACY.  THIS  PROGRAM  WILL 
COMPUTE  THE  FREQUENCY  RESPONSE  OF  UP  TO  THREE  SYSTEMS. 
FOR  EACH  SYSTEM,  THE  USER  HAS  THE  OPTION  OF  HAVING  THE 
OUTPUT  EXPRESSED  IN  DECIBELS  (db).  THE  OUTPUT  OF  THIS 
PROGRAM  IS  STORED  IN  TABULAR  FORM  IN  THE  FILE 
'DIGFREQ.  OUT'  AND  IN  A  FORM  SUITABLE  FOR  PLOTTING 
IN  THE  FILE  ’DIGFREQ.  DAT' . 


INPUT 


THIS  PROGRAM  ASSIMES  THAT  EACH  DISCRETE  SYSTEM  IS  MODELED  BY  THE 
EQUATION:  H(z)  =  num/den  WHERE: 

num  =  b(0)*z**T.  +  b(l)*z**(L-l)  +  ...  +  b(L-l)*z  +  b(L) 

den  =  c(0)*z**N  +  c(l)*z**(N-l)  +  ...  +  c(N-l)*z  +  c(N) 

L  =  A  NON-NEGATIVE  INTEGER,  THE  DEGREE  OF  THE  NUMERATOR 
POLYNOMIAL. 

N  «  A  NON-NEGATIVE  INTEGER,  THE  DEGREE  OF  THE  DENOMINATOR 
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00000000000000000000000000000000000000 


C  POLYNOMIAL. 

C  b(0). . .  b(L)  =  REAL  COEFFICIENTS  OF  THE  NUMERATOR  TERMS. 

C  c(0). . . c(N)  =  REAL  COEFFICIENTS  OF  THE  DENOMINATOR  TERMS. 

C 

C  THE  INPUT  PARAMETERS  SHOULD  BE  STORED  IN  A  FILE  NAMED 
C  'DIGFREQ. IN'.  ALL  OF  THE  READ  STATEMENTS  USED  BY  THIS  PROGRAM 

C  REQUIRE  FORMATTED  INPUT.  PARTICULAR  ATTENTION  SHOULD  BE  PAID 

C  TO  THE  FORMATS,  ESPECIALLY  THE  USE  OF  THE  DECIMAL  POINT  TO 
C  DENOTE  ’REAL'  NUMBERS.  THE  INPUT  PARAMETERS  REQUIRED  BY  THE 
C  PROGRAM  ARE  LISTED  BELOW. 

C 

C 


NAME 

TYPE 

RANGE  (ARRAYS) 

RESTRICTIONS 

numsys 

INTEGER 

1  <=  numsys  <=  3 

L 

INTEGER 

0  <=  L  <=  128 

N 

INTEGER 

0  <=  N  <=  128 

dsorce 

CHARACTER 

'F'  OR  'S’ 

yscal 

CHARACTER 

'STD'  OR  'LOG' 

thetaO 

REAL 

dlthta 

REAL 

numpts 

INTEGER 

1  <=  numpts  <=  101 

b() 

REAL 

0, 

1 ,  2 ,  . . .  ,  L 

•  0  <=  L  <=  128 

c() 

REAL 

o, 

1,  2,  ....  N 

0  N  <=  128 

WHERE: 

numsys  =  THE  NUMBER  OF  DISTINCT  SYSTEMS  H(z)  TO  BE  ANALYZED. 

THIS  INTEGER  VALUE  MUST  OCCUR  AT  THE  TOP  OF  THE  INPUT 
FILE.  IT  DELINEATES  THE  NUMBER  OF  SYSTEMS  TO  BE  READ  BY 
THE  PROGRAM  AND  ANALYZED.  FOR  EACH  SYSTEM  (1,  ....  nurasys) 
THE  PARAMETERS  BELOW  MUST  APPEAR  IN  THE  INPUT  FILE. 

L  =  AN  INTEGER  VALUE  SPECIFYING  THE  DEGREE  OF  THE  NUMERATOR 
POLYNOMIAL. 

N  =  AN  INTEGER  VALUE  SPECIFYING  THE  DEGREE  OF  THE  DENOMINATOR 
POLYNOMIAL. 

dsorce  =  THE  CHARACTER  STRING  * F’  OR  'S'  DENOTING  WHETHER  THE 

SYSTEM  COEFFICIENTS  ARE  TO  BE  READ  FROM  THE  INPUT  FILE  (F) 
OR  GENERATED  (S)  THROUGH  USE  OF  THE  SUBROUTINE  coeff. 

yscal  =»  A  CHARACTER  STRING  SPECIFYING  THE  DESIRED  MAGNITUDE  OPTION: 
'STD'  WILL  PRODUCE  STANDARD  MAGNITUDE  OUTPUT; 

'LOG'  WILL  PRODUCE  MAGNITUDE  EXPRESSED  IN  DECIBELS  (db). 

thetaO  =  THE  STARTING  VALUE  OF  THETA  (RADS)  AS  IN  Z=EXP( J*THETA). 

dlthta  =  THE  INCREMENT  OF  THETA  (RADIANS). 

numpts  =  THE  NUMBER  OF  FREQUENCY  POINTS  FOR  WHICH  THE  OUTPUT  IS 
TO  BE  COMPUTED. 

b()  =  THE  NUMERATOR  COEFFICIENTS  IN  ORDER  b(0),  b(l),  ....  b(L). 

IF  dsorce  “  'F'  IS  SELECTED  THEN  THE  USER  MUST  SUPPLY  THE 
L+l  NUMERATOR  COEFFICIENTS  IN  THE  FILE.  IF  dsorce  «  'S' 
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THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE  NUMERATOR 
COEFFICIENTS  BY  WRITING  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  THE  SPACE  PROVIDED  IN  SUBROUTINE  coeff.  IF  THIS  METHOD 
OF  DATA  GENERATION  IS  ELECTED  THE  PROGRAM  MUST  BE  RECOMPILED 
BEFORE  EXECUTION. 

c()  =  THE  DENOMINATOR  COEFFICIENTS  IN  ORDER  c(0),  c(l),  . ..,  c(N). 
IF  dsorce  =  'F*  IS  SELECTED  THEN  THE  USER  MUST  SUPPLY  THE 
N+l  DENOMINATOR  COEFFICIENTS  IN  THE  FILE.  IF  dsorce  =  ’S’ 
THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE  DENOMINATOR 
COEFFICIENTS  BY  WRITING  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  THE  SPACE  PROVIDED  IN  SUBROUTINE  coeff.  IF  THIS  METHOD 
OF  DATA  GENERATION  IS  ELECTED  THE  PROGRAM  MUST  BE  RECOMPILED 
BEFORE  EXECUTION. 


NOTE:  THE  INPUT  FORMAT  STATEMENTS  OCCUR  IN  THE  MAIN  PROGRAM 

FOLLOWING  THE  CAPTION:  ********  INPUT  FORMAT  ********. 
THE  FORM  OF  THE  INPUT  DATA  FILE  IS: 


LINE  n 

ENTRIES 

FORMAT 

1 

nurasys 

il 

2 

L,N,dsorce,yscal 

i3,tll,i3,t21,al,t31,a3 

3 

dlthta,thetaO,numpts 

2f 10.  0,i3 

NOTE  1 

b(k),  k=0 , 1 , . . . , L 

6f 10. 0 

NOTE  2 
NOTE  3 

c(k) ,  k=Q , 1 , . . . ,N 

6fl0. 0 

WHERE: 

NN  =  1  +  (L/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER). 

ND  =  1  +  (N/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER). 

NOTES  1. 

THE  NEXT  NN  LINES  ARE  ONLY  REQUIRED  IF  dsorce  a  *FT .  IF 

dsorce  =  'S'  THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE 
L+l  NUMERATOR  COEFFICIENTS  IN  THE  SUBROUTINE  coeff. 

THE  USER  MUST  PROVIDE  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  SUBROUTINE  coeff  TO  GENERATE  THE  VALUES  FOR  b(). 

2.  THE  NEXT  ND  LINES  ARE  ONLY  REQUIRED  IF  dsorce  =  'F'.  IF 
dsorce  =  'S'  THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE 
N+l  DENOMINATOR  COEFFICIENTS  IN  THE  SUBROUTINE  coeff. 

THE  USER  MUST  PROVIDE  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  SUBROUTINE  coeff  TO  GENERATE  THE  VALUES  FOR  c(). 

3.  FOR  numsys  >  1  THE  FORMAT  OF  LINES  2. . .  IS  REPEATED. 

4.  THE  FORMAT  f 10. 0  USED  FOR  INPUT  DATA  PERMITS  THE  DECIMAL 
POINT  TO  BE  PLACED  ANYWHERE  IN  THE  FIELD  OF  10  COLUMNS 
AND  ALSO  ALLOWS  THE  EXPONENTIAL  FORMAT  TO  BE  USED  (EG. 
3146.2  «  3. 1462E+03). 


C******************************  OUTPUT  ******************************* 
C 

C  THE  OUTPUT  DATA  CREATED  BY  THE  PROGRAM  IS  STORED  IN  TABULAR  FORM 
C  IN  THE  FILE  'DIGFREQ. OUT' .  ADDITIONALLY,  THE  OUTPUT  DATA  IS 
C  WRITTEN  INTO  THE  FILE  'DIGFREQ. DAT1  TO  FACILITATE  PLOTTING  BY 
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C  A  SEPARATE ,  USER  SUPPLIED  PROGRAM.  THE  FORMAT  OF  THE  DATA  IN 
C  'DIGFREQ.DAT'  IS:  el2.6,  2x,  el2. 6.  THE  FIRST  ENTRY  CORRESPONDS 
C  TO  THE  ORDINATE  VALUE  (THETA)  AND  THE  SECOND  ENTRY  THE  ABSCISSA 
C  VALUE  (MAGNITUDE  OR  PHASE).  ADDITIONAL  HEADER  INFORMATION  IS 
C  WRITTEN  INTO  THE  DATA  FILE  TO  ALLOW  FOR  CONTROL  AND  LABELING  OF 
C  EACH  PLOT. 

C 

C 


EXAMPLE 


THE  INPUT  PARAMETERS  FOR  THE  SYSTEM  DESCRIBED  BELOW  ARE  STORED  IN 
THE  SAMPLE  INPUT  FILE  'DIGFREQ. TST'  AND  CAN  BE  USED  FOR  A  TRIAL 
RUN  IN  THE  TEST  MODE. 


SYSTEM:  H(z)=z/(z-0. 5) 

GOAL:  TO  OBTAIN  THE  FREQUENCY  RESPONSE  FOR  THIS  SYSTEM  FROM 

THETA  =  0.  0  TO  THETA  =  3. 14159  (PI  RADIANS)  IN  STEPS 
OF  dlthta  =  PI/10.  0 

FOR  THE  SYSTEM  DESCRIBED  ABOVE  THE  INPUT  FILE  IS: 

1 

001  001  F  STD 

0.314159  0.0  011 

1.0  0.0 

1.0  -0.5 


THE  RESULTING  OUTPUT  DATA  FILE:  'DIGFREQ. OUT*  IS: 

INPUT  DATA  FOR  SYSTEM  #  1 

INPUT  DATA  SOURCEF I LE:  DIGFREQ. TST 
DEGREE  OF  NUMERATOR  »  1 

DEGREE  OF  DENOMINATOR  =  1 

dsorce  =  F 

NUMBER  OF  FREQUENCY  POINTS  =  11  MAGNITUDE  OPTION  «  STD 
STARTING  VALUE  OF  THETA  =  . OOOOOOE+OO 
INCREMENT  OF  THETA  =  .314159E+00 

THE  NUMERATOR  COEFFICIENTS  b(0) ,b(l). . .  b(L)  ARE 

. 1000E+01  . 0000E+00 


THE  DENOMINATOR  COEFFICIENTS  c(0) ,c( 1). . .  c(N)  ARE 
. 1000E+01  5000E+00 


OUTPUT  DATA  FOR  SYSTEM  #  1 

THETA  MAGNITUDE  PHASE 

(RADIANS)  (DEGREES) 
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.  OOOOOOE+OO 
.  3141S9E+00 
.  628318Z+00 
.  942477E+00 
. 125664E+01 
. 157080E+01 
. 188495E+01 
.  219911E+01 
. 251327E+01 
. 282743E+01 
. 314159E+01 


.  200000E+01 
. 182897E+01 
. 150588E+01 
. 122S86E+01 
. 103088E+01 
.  894428E+00 
. 800894E+00 
. 737654E+00 
. 696900E+00 
. 674Q38E+00 
. 666667E+00 


. OOOOOOE+OO 
164149E+02 
-. 262677E+02 
29807E+02 
293546E+02 
265651E+02 
223862E+02 
-. 173608E+02 
-. 118186E+02 
597793E+01 
484184E-04 


END  OF  RUN,  SYSTEM  #  1 


FOR  ILLUSTRATIVE  PURPOSES  THE  COEFFICIENTS  b()  AND  c()  COULD 
HAVE  BEEN  GENERATED  BY  SPECIFYING  dsorce  =  'S*  AND  WRITING  THE 
APPROPRIATE  FORTRAN  STATEMENTS  INTO  SUBROUTINE  coeff.  THE 
STATEMENTS  THAT  COULD  BE  USED  TO  ACCOMPLISH  THIS  ARE  WRITTEN  INTO 
THE  SUBROUTINE  BUT  ARE  'COMMENTED  OUT'. 

MAIN  PROGRAM 


character  infile*12,  mode*l,  ylabl*13,  dsorce*! ,  yscal*3 
real  mh(101),  ph(101),  thetav(lOl),  c( 0:128),  b(0: 128) 

C  PROMPT  USER  FOR  MODE:  BATCH  OR  TEST. 

write(*, 1115) 
read(,v,  1117)  mode 

if  ((mode.  eq.  *Y*  ).or.  (mode.  eq.  'y' ))  then 
mode  =  ' Y' 
write(,1f,  1118) 
read(*,1119)  infile 
else 

infile  =  'DIGFREQ. IN' 
end  if 

C  UNIT=1  DEFINED  AS  INPUT  FILE.  UNITS=2,3  DEFINED  AS  OUTPUT  FILES. 

open(unit=l,files=infile,status='old' ,iostat®ierr,err-999) 
open(unit=2 , f ile=' DIGFREQ.  OUT' ) 
open(  un  i  t=3 ,  f  i  1  e=> '  DIGFREQ .  DAT '  ) 

C  READ  INPUT  PARAMETERS  AND  CONDUCT  ERROR  CHECKS. 

read( 1,1000)  numsys 
numplts  =  numsys+2 
vrite(3,2000)  numplts 
if((numsys.  It.  l).or.  (numsys.  gt.  3))  then 
writef*, 1122)  numsys 

stop  'Error,  numsys  must  be  in  the  range:  1  <»  numsys  <-  3.  ' 
end  if 
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do  10  nsys=l,  numsys 

data  mh/ 101*0.0/,  ph/101*0.0/,  thetav/101*0. 0/ 
data  b/ 129*0.  0/,  c/ 129*0.0/ 

read( 1,1001)  L,  N,  dsorce,  yscal 
read( 1,1002)  dlthta,  thetaO,  numpts 

lf((L.  It.  0),  or.  (L.  gt.  128))  then 
writer*, 1124)  nsys,  L 

stop  'Error,  L  must  be  in  the  range:  0  <=  L  <=  128. ' 
elseif  ( (N.  It.  0).  or.  (N.  gt.  128)  )  then 
write(*, 1125)  nsys,  N 

stop  'Error,  N  must  be  in  the  range:  0  <=  N  <=  128. * 
endif 

if (  ( dsorce.  eq. '  F ' ) .  or.  (  dsorce.  eq. '  f ' ) )  then 
dsorce  =  'F' 

elseif ((dsorce. eq. *S' ).  or.  (dsorce. eq. 's' ))  then 
dsorce  -  'S' 
else 

writer*, 1018)  dsorce 

stop  'The  allowed  values  for  dsorce  are:  "S"  or  "F".' 
endif 

if( (numpts. It. l).or. (numpts. gt. 101))  then 
writer*, 1127)  nsys,  numpts 

stop  Error,  numpts  must  be  in  the  range:  1  <=  numpts  <«  101. ' 
endif 

if((yscal. eq.  'STD'). or.  (yscal. eq. 'std'))  then 
yscal  =  'STD' 
ylabl  «  ’  MAGNITUDE  ' 

elseif((yscal,  eq,  'LOG* ). or. (yscal. eq. 'log'))  then 
yscal  «  'LOG' 
ylabl  o  'MAGNITUDE(db)' 
else 

writer*, 1128)  yscal 

stop  'Error,  yscal  must  be  the  string:  ''LOG1'  or  "STD".’ 
endif 

C  FOR  dsorce  =  'F'  READ  THE  COEFFICIENTS  b()  AND  c( )  FROM  THE  INPUT 
C  FILE.  FOR  dsorce  *  'S'  CALL  coeff  TO  GENERATE  THE  COEFFICIENTS. 

if(dsorce. eq. 'F' )  then 
read( 1,1003)  (b(k),k*0,L) 
rcad(  1,1003)  (c(k) ,k=*0,N) 
else 

call  coeff(L,N,nsys,b,c) 
endif 

C  WRITE  INPUT  DATA  INTO  THE  OUTPUT  FILE:  DIGFREQ.  OUT. 

write(2,1008)  nsys 
write(2,10l0)  infile 
write(2,1110)  I- 
writo(2,llll)  N 
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write(2,1019) 
write(2,1112) 
write(2,1113) 
write(2,1114) 
write(2,1004) 
write(2, 1005) 
write(2,1006) 
write(2, 1005) 
write(2, 1009) 
write(2,1126) 
write(2,1007) 


dsorce 

numpts,  yscal 

thetaO 

dlthta 

(b(k),k=0,L) 

(c(k) ,k=0,N) 

nsys 

ylabl 


C  FOR  TEST  MODE  ECHO  ALL  INPUTS  ONTO  MONITOR  (UNIT  =  *). 

if (mode.  eq.  'Y* )  then 
write(*,1120)  nsys,  infile 
write(*,1110)  L 
write(*, 1111)  N 
write(*,1019)  dsorce 
write(*,1112)  numpts,  yscal 
write(*,1113)  thetaO 
write(*,1114)  dlthta 
write(*, 1004) 

write(*,10Q5)  (b(k),k=Q,L) 
write(*»1006) 

write(*,1005)  (c(k),k=0,N) 
write(*, 1123)  nsys 

pause  'END  OF  RUN,  STRIKE  <CR>  WHEN  READY  TO  CONTINUE.  ' 
end  if 


C  CALL  dfresp  TO  COMPUTE  THE  FREQUENCY  RESPONSE. 

call  dfresp( b,c,mh,ph,L,N, thetaO , dlthta , thetav , numpts , ysca 1 ) 

C  WRITE  RESULTS  INTO  OUTPUT  FILE:  DIGFREQ.DAT. 

write( 3,2001)  numpts 
vrite( 3,*)  'MAGNITUDE  RESPONSE' 
write(3,*)  'THETA  (RADS)' 
writo(3,2003)  ylabl 
do  55  npal ,  numpts 
writo(3,2010)  thetav(np),  mh(np) 

55  continue 

write(3,2001)  numpts 
write( 3,*)  'PHASE  RESPONSE' 
write(3,*)  'THETA  (RADS)' 
write(3,2003)  '  PHASE  (DEG)  ' 
do  56  np=l,  numpts 
write(3,201Q)  thetav(np),  ph(np) 

56  continue 


C  WRITE  RESULTS  INTO  OUTPUT  FILE:  DIGFREQ.  OUT. 
do  150  np28!,  numpts 
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write(2,1013)  thetav(np),  mh(np),  ph(np) 
ISO  continue 
write(2,1123)  nsys 

10  continue 

write(*,1121) 

999  close(unit=l) 
close(unit=2) 
close(unit=3) 

iffierr.  gt.  0)  then 
writef*, 1116)  ierr 
endif 

C***********  INPUT  FORMAT  t^rtr***********-*-* 


1000  format( il) 

1001  format( i3, til, i3ft21,al,t31,a3) 

1002  formate 2f 10.0,i3) 

1003  formatebflO. 0) 


1004  formate t4/ THE  NUMERATOR  COEFFICIENTS  b(0) ,be 1) _ bCL)  ARE: 

1005  formate 6e  2X,ell. 4) ,//) 

1006  0 for mate //,t4,' THE  DENOMINATOR  COEFFICIENTS  c(Q),c(l). , . c(N)' 
V  ARE: 


',/) 


1007 

1003 

1009 

1010 
1013 
1018 
1019 
1110 
mi 


{ ornate 1 6 , 1 ( RADIANS ) ’ , t38 , ' ( DEGREES) '  /) 
formate  tie,1  INPUT  DATA  FOR  SYSTEM  9  \il ,//) 


formate///, tl6, '  OUTPUT  DATA  FOR  SYSTEM  $ 
format (t4/ INTO  DATA  SOURCEFILE:  *  ,al2) 


.il ./) 


format( t4,3(el2, 6,4x1) 

£ormat(  Ix/dsorco  38  ,al ,2x,' Error ,  illegal  value  for  dsorce,  ' ) 
format( t4, 'dsorce  -  '  ,al) 
formate t4, 'DEGREE  OF  NUMERATOR  *  \i3) 
formate t4,' DEGREE  OF  DENOMINATOR  «*  ’,13) 

1112  0format(t4/NtrSES  OF  FREQUENCY  POINTS  »  '  ,i3,t39 /MAGNITUDE ' , 

1*  OPTION  «  ' ,a?) 

1113  formate t4/ STARTING  VALUE  OF  TKETA  *  \el2.6) 

1114  formate t4/ INCREMENT  OF  THETA  =  ',©12.6,/) 

1115  0 formate  lx /o0  YOU  WISH  TO  RUN  THIS  PROGRAM  IN  TEST' , 

1*  MODE  ?  (Y/N)  <CR>  :  ,) 

1116  Oforoat(///,lx, 'ERROR  OPENING  INPUT  FILE,  PROGRAM  TERMINATED.', 
1// ,  lx /ERROR  CODE:  *,i4, /////) 


1117  format(cl) 

1118  0foroat(/////, ix/TYPE  THE  NAME  OF  YOUR  DATA  FILE  FOLLOWED' 

1'  BY  <CR>. ' ,/,'  IF  YOU  DESIRE  TO  MAKE  A  TEST  RUN  USING  THE*, 

2*  SAMPLE  DATA  ALREADY  STORED',/.'  IN  THE  FILE:  DIGFREQ.TST’ , 

3  '  TYPE:  DIGFREQ.TST  <CR>\/,'  FILENAME:  ' ,  ,) 

1119  foroatfal2) 

1120  formate////, t4,’SYSTEM  9  ',11/  INPUT  DATA  SOURCEFILE:  ' ,al2) 

1121  0format(//,t4/TABULAR  OUTPUT  DATA  IS  STORED  IN  FILE:  DIGFREQ. OUT* 
l/,t4,' PLOTTING  DATA  IS  STORED  IN  FILE:  DIGFREQ.  DAT.  ') 

1122  foroat(//////, t2/The  value  of  nuosys  is:  ',11,'.') 


» 
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1123  format(/, lx, 13( f f  END  OF  RUN,  SYSTEM  #' ,il,2x, 13( ,//) 

1124  0format(//////,t2, 'The  degree(L)  of  the  numerator  for  system  ', 

1  *#  Ml,’  is  :  L  =  ',i3,\  ’) 

1125  0format(//////,t2, 'The  degree(N)  of  the  denominator  for  system' 

1,'  #  * ,il,  is  :  N  =  ' ,i3, ' .  ' ) 

1126  format(///, t8,,THETA',t21,al3,t40, 'PHASE') 

1127  formate////// ,t2, 'The  value  of  nurapts  for  system  ',il,*  is:  ',i3) 

1128  formate/////, t2.  The  value  of  yscal  is:  ',a3,'.') 

2000  formate il) 

2001  formate i3) 

2003  formate al3) 

2010  formateel2. 6,2x,eJ2. 6) 

end 


C 


SUBROUTINE:  dfresp 


C  PURPOSE: 

C 

C 

C 

C 


THIS  SUBROUTINE  COMPUTES  THE  FREQUENCY  RESPONSE  OF 
THE  SYSTEM.  ALU  FREQUENCY  CALCULATIONS  ARE  IN  RADIANS, 
HOWEVER  THE  OUTPUT  IS  CONVERTED  TO  DEGREES. 

THE  OUTPUT  FORMAT  FOR  EACH  FREQUENCY  INCREMENT  IS: 
MAGNITUDEeM)  PKASiXP)  AS  IN;  H*£XP(J'*P). 


subroutine  dfresp(b,c,mh,ph,L,N,thetaQ,dlthta,thetav,numpts, yscal) 

real  mh(mimpts),  phenumpts),  thetev(nmapts) ,  ims,  rez 

real  b(Q:L),  c(0:N) 

character  yscal*3 

complex  z,  den,  nuo,  h,  cl 

C  DEFINE  CONSTANTS. 

pi  »  4.0*atanei.O) 
ci  «  (1.0, 0.0) 

C  ITERATE  FROM  thotaO,  IN  INCREMENTS  OF  dlthta. 

do  100  np=«l ,  numpts 
num  »  ei*b(0) 
den  =  ci^ctO) 

thotav(np)  -  thotaO  +  enn-O^dlthta 
rez  =  cos(thetavenp)) 

ios  ®  sinethetav(np)) 

&  a  cojplxerez.ioz) 

C  CALCULATE  NUMERATOR  FOR  GIVEN  VALUE  OF  THETA,  IF  L  >  0. 

if(L.  gt.  0)  then 
do  50  k=l,  L 
nua  =  2*nua  +  c  k) 

50  continue 


end  if 


oooo 


C  CALCULATE  DENOMINATOR  FOR  GIVEN  VALUE  OF  THETA,  IF  N  >  0. 


if(N. gt. 0)  then 
do  70  k=l,  N 

den  =  z*den  +  ci*c(k) 

70  continue 

end  if 

h  =  num/den 

C  CONVERT  COMPLEX  VALUE  V  INTO  MAGNITUDE(mh)  AND  PHASE(ph)  TERMS. 

C  IF  yscal  «  'LOG*  THEN  CONVERT  MAGNITUDE  TO  DECIBELS  (db). 

C  DIVIDE  BY  ZERO  AVOIDED  BY  ’if’  STATEMENTS. 

mh(np)  =  cabs(h) 

if (yscal. eq. 'LOG' )  then 
if(mh(np).  gt.  0. 00001)  then 
mh(np)  =  20. Q*loglQ(mh(np)) 
else 

tnh(np)  =  -100. 0 

endif 

endif 

if(abs(real(h)). It. 1.  Oe-15)  then 

if(abs(aitoag(h)).  le.  1.  Oe-15)  ph(np)=0.  0 
if(aimsg(h).  gt.  1.  Oe-15)  ph(np)=9Q.  0 
if(airaag(h).  It.  -1,  Oe-15)  ph(np)»-90.  0 
else 

ph(np)  =  (160. 0/pi )*atan2(aimag(h),real(h)) 
endif 

100  continue 


return 

end 


C 


SUBROUTINE:  coeff 


PURPOSE:  THIS  SUBROUTINE  ALLOWS  THE  USER  TO  GENERATE  THE 

NUMERATOR  AND  DENOMINATOR  COEFFICIENTS  THAT  DESCRIBE 
EACH  SYSTEM  TO  BE  ANALYZED.  IF  dsorce  «  ’S'  THEN 
THE  MAIN  PROGRAM  WILL  CALL  THIS  SUBROUTINE. 


subrout  ine  coe  f  f ( I- ,  N ,  nsys ,  b ,  c ) 
real  b(0:L),  c(0:N) 


pi  «  4. 0*atau( 1.  0) 
el  a  L 
en  »  N 


C  DEVELOP  THE  EQUATIONS  TO  GENERATE  VALUES  FOR  THE  ARRAYS  b( )  AND  c( ) 
C  IN  THIS  SPACE.  THE  STATEMENTS  TYPED  IN  MUST  FOLLOW  STANDARD 
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FORTRAN  77  RULES  AND  MAY  USE  FORTRAN  77  INTRINSIC  FUNCTIONS  SUCH  AS: 
SINQ,  C0S(),  ABS( ). . .  AN  EXAMPLE  IS  SHOWN  BELOW.  NOTE  THAT  THE 
VALUE  nsys  CAN  BE  USED  TO  DISTINGUISH  BETWEEN  SYSTEMS  IF  MORE  THAN 
ONE  SYSTEM  (numsys  >  1)  IS  TO  BE  ANALYZED. 

***  EXAMPLE  *** 

if  (nsys.  eq.  1)  then 
do  2  i=0,  L 

b(i)  =  cos(i*pi/(2.0*el)) 

2  continue 

do  3  i=0,  N 

c(i)  =  cos(2. 0*i*pi/(3. 0*en)) 

3  continue 
endif 


return 

ond 


B.  MAGNITUDE,  PHASE  Or  SECOND  ORDEii  HR  NOTCH  FILTER:  P-0 

INPUT  DATA  FOR  SYSTEM  it  l 


INPUT  DATA  SOURCEFILE.  digireq. tst 
DEGREE  OF  NUMERATOR  =  2 

DEGREE  OF  DENOMINATOR  *  2 

dsorce  *  F 

NUMBER  OF  FREQUENCY  POINTS  «*  100  MAGNITUDE  OPTION  «*  STD 
STARTING  VALUE  OF  THETA  -  . OOOOOOE+OC 
INCREMENT  OF  THETA  •  . 280600E-01 

THE  NUMERATOR  COEFFICIENTS  b(0) ,b( 1). . .  b(L)  ARE: 

. 9900E+00  -. 3299E+00  . 9900E+00 


THE  DENOMINATOR  COEFFICIENTS  c(0),c(l). ..c(N)  ARE: 
. 100QE+01  -.  3299E+00  .9801E+00 


SS 


OUTPUT  DATA  FOR  SYSTEM  #  1 


THETA 

(RADIANS) 

. OGOOOCE+OG 
. 280600E-01 
.S61200E-01 
.8' 18J0E-01 
.  112240E+0') 
,140r00E+O0 
. 168360E+00 
, 196420E+00 
. 2244S0E+00 
. 252540E+00 
. 280600E+00 
. 306660E+00 
. 336720E+0Q 
. 364780E+00 
. 392840E+00 
. 420900E+00 
. 448960E+00 
•477020E+00 
. 505080E+00 
. 533140E+00 
. 561200E+00 
. 58926QE+OQ 
. 6I7320E+00 
. 645380E+00 
. 67344OE+Q0 
. 701500E+00 
. 729560E+00 
. 75762QE+00 
. 78S680E+00 
. 613740E+00 
. 841800E+00 
. 869860E+00 
. 897920E+00 
. 925980E+00 
. 9S4040E+00 
. 9821Q0E+00 
. 101016E+01 
. 103822E+01 
. 1O6628E+0I 
. 109434E+01 
. I12240E+01 
.115046E-H)1 
.  117852E+01 
.  1206S8E+OI 
.  123464E+01 
.  126270E+0I 
.  129076E+01 


MAGNITUDE 


. 999939E+00 
. 999939E+00 
. 999939E+00 
. 999939E+00 
. 999938E+00 
. 999938E+00 
. 999937E+00 
. 999936E+00 
. 999935E+00 
. 999934E+00 
. 999933E+00 
. 999931E+00 
, 999930E+0Q 
. 999928E+00 
. 999926E+00 
. 999923E+00 
. 99992UE+OU 
. 999917E+00 
. 999914E+00 
. 999910E+00 
. 999906E+00 
, 999902E+00 
. 999896E+00 
. 99989 1E+00 
. 999884E+O0 
. 999676E+00 
. 99986SE+00 
. 9998S8E+00 
. 999847E+00 
.999833E+00 
. 999S20E+00 
. 999803E+00 
. 999783E+00 
. 999759E+00 
.999731E+00 
. 999697E+00 
. 999655E+00 
. 999603E+00 
. 999537E+O0 
. 999453E+00 
. 999343E+00 
. 999193E+00 
. 998985E+00 
. 998681E+00 
.998213E+00 
.  997440E+00 
.  996024E+00 


PHASE 

(DEGREES) 

. OOOOOOE+OO 
193947E-01 
388289E-01 
S83436E-01 
779784E-01 
977754E-01 
117779E+00 
138030E+00 
158581E+00 
179474E+00 
200769E+00 
222516E+00 
244780E+00 
267620E+00 
291112E+00 
31S332E+00 
34U3bbh+0U 
366300E+00 
393251E+00 
421327E+00 
450663E+OO 
-.481405E+00 
5I3723E+00 
547804E+00 
». 583872E+00 
622182E+00 
663033E+00 
706764E+00 
753788E+00 
804587E+00 
-.859743E+00 
919964E+00 
966100E+00 
10592 1E+01 
11406 1E+01 
123199E+01 
133548E+01 
1453S9E+01 
15909SE+01 
175176E+01 
194344E+01 
217626E+01 
246S60E+O1 
283559E+01 
332633E+01 
*>.  400969E+01 
502355E+01 
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.  131882E+01 
.  134688E+01 
. 137494E+01 
.  140300E+01 
. 143106E+01 
. 145912E+01 
.  148718E+01 
. 151524E+01 
.  154330E+01 
. 157136E+01 
. 159942E+01 
. 162748E+01 
. 16S5S4E+01 
. 168360E+01 
. 171166E+01 
. 173972E+01 
. 176778E+01 
. 179584E+01 
. 182390E+01 
. 185196E+01 
. 188002E+01 
. 190808E+01 
. 193614E+01 
. 196420E+01 
. 199226E+01 
,  ?o?rmE+m 
. 204838E+01 
. 207644E+01 
. 210450E+01 
. 213256E+01 
. 216062E+01 
. 218868E+01 
, 221674E+01 
. 224480E+01 
. 227286E+01 
. 230092E+01 
. 232898E+01 
. 23S704E+01 
. 238510E+01 
. 241316E+01 
. 244122E+01 
. 246928E+01 
. 249734E+01 
. 252540E+01 
.25S346E+01 
.  256152E+01 
. 2609S8E+01 
.263764E+01 
. 266S70E+01 
, 269376E+01 
. 272182E^01 
.274988E+01 
. 277794E+01 


. 992994E+00 
. 984530E+00 
.  94285 7E+00 
.  381903E-01 
.  939926E+00 
.  984111E+00 
.  992867E+00 
. 995971E+00 
.  997413E+00 
.  998198E+00 
.  998672E+00 
.  998979E+00 

•  • ■  ;  999190E+00 
.  999341E+00 
.  999452E+00 
.  999537E+00 
.  999603E+00 
.  999655E+00 
.  999698E+00 
.  999732E+00 
.  999761E+00 
. 999785E+00 
.  999805E+00 
. 999822E+00 
.  999837E+00 
qqQAsnE+nn 
. 999861E+00 
.  999871E+00 
.  999880E+00 
. 999888E+00 
.  999894E+00 
. 999901E+00 
. 999906E+00 
. 999911E+00 
.  999915E+00 
. 999919E+00 
. 999923E+00 
. 999926E+00 
. 999929E+00 
. 99f >32E+00 
. 999934E+00 
. 999937E+00 
. 999939E+00 
. 999941E+00 
. 999942E+00 
. 999944E+00 
.  999945E+00 
. 999947E+O0 
.  999948E+00 
. 999949E+00 
. 999950E+00 
. 99995 1E+00 
. 999952E+00 

END  OF  RUN.  SYSTEM 


. 671241E+01 
. 100259E+02 
. 194059E+02 
. 877665E+02 
. 200014E+02 
. 102595E+02 
. 687161E+01 
. 516115E+01 
. 413001E+01 
. 344006E+01 
. 294546E+01 
. 257310E+01 
. 228226E+01 
. 204852E+01 
. 185629E+01 
. 169521E+01 
. 155808E+01 
. 143978E+01 
. 133653E+01 
. 124550E+01 
.  116455E+01 
.  109198E+01 
.  102647E+01 
.  966960E+00 
.912579E-HX) 

. 816531E+OQ 
. 773792E+00 
. 734015E+00 
. 696843E+00 
. 661967E+00 
. 629193E+00 
. 598246E+00 
. 56895OE+Q0 
.541144E+00 
. 514683E+00 
. 489434E+00 
. 465265E+00 
. 442146E+00 
.  419918E+00 
.  398521E+00 
.  377889E+00 
.  3S7941E+00 
.  338639E+00 
.  319915E+00 
.  301724E4-00 
.  284024E-KH3 
.  266771E-K)0 
.  249929E+00 
.  2334S9E+0O 
. 217332E+00 
. 201519E+00 
. 185987E+00 

#  1  . 
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C.  MAGNITUDE,  PHASE  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P=  1 


INPUT  DATA  FOR  SYSTEM  #  1 


INPUT  DATA  SOURCEFILE:  digf req. tst 
DEGREE  OF  NUMERATOR  =  3 

DEGREE  OF  DENOMINATOR  =  3 

dsorce  =  F 

NUMBER  OF  FREQUENCY  POINTS  =  100  MAGNITUDE  OPTION  =  STD 
STARTING  VALUE  OF  THETA  =  .  OOOOOOE+OO 
INCREMENT  OF  THETA  =  . 280600E-01 

THE  NUMERATOR  COEFFICIENTS  b(0) ,b( 1). . .  b(L)  ARE: 

. 99Q0E+QQ  ». 3299E-02  .8811E+00  .  3266E+QQ 


THE  DENOMINATOR  COEFFICIENTS  c(0) ,c( 1). . . c(N)  ARE: 
. 100QE+01  , 00QQE+00  .8712E+00  .3234E+00 


OUTPUT  DATA  FOR  SYSTEM  #  1 


THETA 

(RADIANS) 

.  000000E+00 
. 280600E-01 
. 56120QE-01 
. 841S00E-01 
. 112240E+00 
. 140300E+00 
.  168360E+00 
. 196420E+00 
. 224480E+00 
. 252540E+00 
. 280600E+00 
. 308660E+00 
. 336720E+00 
. 364780E+00 
. 39284QE+00 
. 42090CE+00 
. 448960E+00 
. 477020E+00 


MAGNITUDE 


. 999940E+00 
.  999940E+00 
.  999939E+00 
.  999939E+00 
.  999939E+00 
.  999938E+00 
. 999937E+00 
. 999936E+00 
. 99993SE+00 
. 999934E+00 
.999933E+00 
. 999932E+00 
. 999930E+00 
. 999928E+O0 
. 999926E+00 
. 999923E+00 
. 39992 1E+00 
. 999918E+00 


PHASE 

(DEGREES) 

. OOOOOOE+OO 
193950E-01 
3$8293E*01 
-.583432E-01 
77979SE-01 
-.977766E-01 
117780E+00 
13803 1E+00 
158S79E+00 
179475E+00 
**.  200768E+00 
222515E+00 
244780E+00 
267616E+00 
291112E+00 
315332E+00 
340364E+00 
366303E+00 


505080E+00 

533140E+00 

561200E+00 

589260E+00 

617320E+00 

645380E+00 

673440E+00 

701500E+00 

729560E+00 

7S7620E+00 

785680E+00 

813740E+00 

841800E+00 

869860E+00 

897920E+00 

925980E+00 

9S4040E+00 

982100E+00 

1G1016E+01 

103822E+01 

106628E+01 

109434E+01 

112240E+01 

115046E+01 

117852E+01 

120658E+01 

123464E+01 

126270E+01 

129076E+01 

131882E+01 

134688E+01 

137494E+01 

140300E+01 

143106E+01 

145912E+01 

14671SE+01 

1S1524E+01 

1S4330E+01 

1S7136E+01 

1S9942E+01 

162748E+01 

16S554E+01 

163360E+01 

171166E+01 

173972E+01 

176778E+01 

179584E+01 

182390E+C1 

185196E+01 

188002E+01 

190808E+01 

193614E+01 

196420E+01 

199226E+01 

202032E+01 

204638E+01 


. 999914E+00 
.  999911E+00 
.  999906E+00 
.  999902E+00 
. 999897E+00 
. 999891E+00 
.  999884E+00 
.  999877E+00 
. 999868E+00 
.  999859E+00 
.  999848E+00 
.  999835E+00 
. 999820E+00 
. 999803E+00 
. 999783E+00 
. 999759E+00 
. 999731E+00 
. 999697E+00 
. 99965SE+00 
. 999603E+Q0 
. 999538E+00 
. 999454E+00 
. 999343E+00 
. 999194E+0Q 
. 998985E+00 
. 998681E+00 
. 998214E+00 
. 997441E+00 
. 996025E+00 
.  992995E+0Q 
. 984532E+00 
. 942860E+Q0 
. 382032E-01 
.  939922E+00 
. 984110E+00 
. 992866E+Q0 
. 995970E+00 
. 997413E+00 
. 998197E+00 
. 998671E+00 
. 998979E+00 
. 999189E+00 
. 999340E+00 
. 999452E+00 
. 999537E+00 
. 999603E+00 
. 9996SSE+00 
, 999697E+00 
. 999732E+00 
. 999761E+00 
. 99978SE+00 
. 999805E+00 
. 999822E+00 
. 999837E+00 
. 9998S0E+00 
. 99986 1E+00 


393250E+00 
421327E+00 
450665E+00 
481406E+00 
513721E+00 
». 547805E+00 
583875E+00 
622183E+00 
663032E+00 
706764E+00 
753788E+00 
804590E+00 
859746E+00 
919964E+00 
986101E+00 
105921E+01 
114061E+01 
123199E+01 
133S48E+01 
145389E+01 
159095E+01 
17S176E+01 
194344E+01 
217625E+01 
246560E+01 
283S59E+01 
332632E+01 
400968E+01 
502853E+01 
*. 671239E+01 
1002S9E+02 
1940S8E+02 
877594E+02 
. 200013E+02 
.  102595E-K52 
.  667158E+01 
.  516U2E+01 
. 412999E+01 
. 344004E+01 
. 294544E+01 
. 257309E+01 
. 22822SE+01 
, 204850E+01 
.  18S627E+01 
. 169S20E+01 
.  1S5807E+01 
.  143977E+01 
. 133652E+01 
. 124549E+01 
. 1164S4E+01 
. 109197E-HJ1 
. 102646E+01 
. 966947E+00 
. 912566E+0G 
.  862621E-HJ0 
. 816S22E+00 
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. 207644E+01 
. 210450E+01 
. 213256E+01 
. 216062E+01 
. 218868E+01 
. 221674E+01 
. 224480E+01 
. 227286E+01 
, 230092E+01 
. 232898E+01 
. 235704E+01 
. 238510E+01 
. 241316E+01 
. 244122E+01 
. 246928E+01 
. 249734E+01 
. 252540E+01 
. 255346E+01 
. 258152E+01 
,  260958E+01 
. 263764E+01 
. 266570E+01 
. 269376E+01 
. 272182E+01 
. 274988E+01 

*  ^  n  n+t  #  w%  t  r\  <4 
•  <./// 


. 999871E+00 
. 999880E+00 
. 99988 7E+00 
. 999894E+00 
. 999900E+00 
. 999906E+00 
. 999911E+00 
. 999915E+00 
. 999919E+00 
. 999923E+00 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 999934E+00 
. 999936E+00 
. 999939E+00 
. 999940E+00 
. 999942E+00 
. 999944E+00 
. 99994SE+00 
. 999947E+00 
. 999948E+00 
.  999949E+00 
. 999950E+00 
.99995 1E+00 
.  5S3S32E4-00 


. 773783E+00 
. 734003E+00 
. 696833E+00 
. 661978E+00 
. 629187E+00 
. 598235E+00 
. 568938E+00 
. 541135E+00 
. 514669E+00 
. 489423E+00 
. 465277E+00 
. 442138E+00 
. 419907E+00 
. 398514E+00 
. 377878E+00 
. 357937E+00 
. 338629E+00 
. 319905E+00 
. 301716E+00 
. 284016E+00 
. 266763E+00 
. 249920E+00 
. 233452E+00 
. 217326E+00 
. 201509E+00 


END  OF  RUN,  SYSTEM  V  1 


D.  MAGNITUDE,  PHASE  OF  THE  AUGMENTED  HR  NOTCH  FILTER:  P» 


INPUT  DATA  FOR  SYSTEM  V  1 


INPUT  DATA  50URC“FILE:  di$freq.tst 
DEGREE  OF  NUMERA*.  3R  ■  4 

DEGREE  OF  DENOMINATOR  *  4 

dsorce  =  F 

NUMBER  OF  FREQUENCY  POINTS  *  100  MAGNITUDE  OPTION  «  STD 
STARTING  VALUE  OF  THETA  -  . OOOOOOE+OO 
INCREMENT  OF  THETA  =  .  280600E-01 

THE  NUMERATOR  COEFFICIENTS  b(0),b(l). . .  b(L)  ARE: 

. 9900E+00  -.3299E-02  . 1861E-01  .6141E+00  -.8625E+00 


THE  DENOMINATOR  COEFFICIENTS  c(0) ,c(l). . . c(N)  ARE: 

. 1000E+01  . OOOOE+OO  .OOOOE+OO  .6108E+00  -.8539E+00 


OUTPUT  DATA  FOR  SYSTEM  #  1 


THETA 

(RADIANS) 

. OOOOOOE+OO 
. 280600E-01 
. 561200E-01 
. 841800E-01 
.  11224QE+00 
.  140300E+00 
.  168360E+00 
. 196420E+00 
.  224480E+00 
. 252540E+00 
. 280600E+00 
.  308660E+0Q 
.  336720E+Q0 
. 364780E+00 
. 392S4CE+CO 
. 420900E+00 
. 448960E+00 
. 477020E+00 
. 505030E+00 
. 533140E+00 
.561200E+00 
. S89260EfOO 
. 617320E+00 
. 645380E+00 
. 67344GE+00 
.  7Q15O0E+OO 
•  «:i$60e+00 
. 75?62QE+:>f> 
. 785680E+Q0 
. 813740E+00 
.641800E+00 
. S69860E+00 
. 897920E+00 
. 925980E+00 
,  9S4040E+00 
. 982100E+00 
. 1Q1016E+01 
. 103822E+01 
. 106628E+01 
. 109434E+01 
. 112240E+01 
. 11S046E+01 
. 117652E+01 
. 120658E+01 
. 123464E+01 


MAGNITUDE 


. 999940E+00 
. 999940E+00 
. 999940E+00 
. 999940E+00 
. 999939E+00 
. 999938E+00 
. 999938E+00 
. 999937E+00 
. 999936E+00 
. 999934E+00 
. 999933E+00 
. 999932E+00 
. 999930E+0Q 
. 999928E+00 
. 992226E+00 
. 999923E+00 
. 999920E+00 
. 999917E+00 
. 999914E+00 
» 999910E+0Q 
. 999906B+00 
. 999901E+0Q 
. 999896E+00 
. 99989 lE+OO 
. 9998S4E+00 
. 999876E+00 
. 99966SE+00 
. 999SS8E+00 
. 999847E+00 
. 999835E+00 
. 999820E+00 
.  999803E+00 
. 999783E+00 
. 999759E+00 
. 999731E+00 
. 999697E+00 
. 9996S5E+00 
. 999603E+00 
. 999537E+00 
. 999453E+00 
. 999343E+00 
. 999193E+00 
. 998985E+00 
. 99868 1E+00 
. 998213E+00 


PHASE 

(DEGREES) 

.  OOOOOOE+OO 
194012E-01 
388423E-01 
*. 583613E-01 
780022E-01 
978035E-01 
1178Q2E+Q0 
■. 138056E+00 
1S8606E+00 
I79500E+00 
200794E+00 
222545E+00 
24480SE+00 
267638E+00 
291136E+CC 
31535 1E+00 
3403S6E+00 
36632 1E+00 
393270E+00 
421345E+0Q 
-.450679E+00 
461422E+O0 
*. 513734E+00 
54781SE+00 
-.583885E+00 
*. 622197E+00 
663044E+00 
706773E+00 
753796E+00 
804595E+00 
85975 1E+00 
919967E+00 
*. 986102E+00 
10592 1E+01 
114062E+01 
123199E+01 
133S48E+01 
145388E+01 
159094E+01 
17517SE+01 
194343E+01 
217625E+01 
<*.  246558E+01 
283557E+01 
332630E+01 
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*- 


.  126270E+01 
. 129076E+01 
. 131882E+01 
.  134688E+01 
. 137494E+01 
. 140300E+01 
. 143106E+01 
. 145912E+01 
. 148718E+01 
. 151524E+01 
. 154330E+01 
. 157136E+01 
. 159942E+01 
. 162748E+01 
.  165554E+01 
. 168360E+01 
. 171166E+01 
. 173972E+01 
. 176776E+01 
. 179584E+01 
. 182390E+01 
. 185196E+01 
,  J.88002E+01 
. 190808E+01 
. 193614E+G1 

1  rtA  A 

. 199226E+01 
. 202032E+01 
. 204836E+01 
. 207644E+01 
. 2104S0E+01 
.253256E+01 
. 216062E+01 
. 21S868E+01 
.221674E+01 
.2244802+01 
. 227286E+01 
.230092E+01 
* 2329962+01 
. 23S704E+01 
. 23SSIOE+01 
. 2413162+01 
>  2441222401 
. 246928E+01 
. 249734E+01 
. 252540E401 
.2SS346E401 
.2S81S2E+01 
. 260958E+01 
.263764E401 
. 2665702401 
. 269376E+01 
. 2721822401 
. 274988E401 
. 277794E+01 


.  997440E+00 
.  996024E+00 
. 992994E+00 
. 984530E+00 
.  942856E+00 
. 3818872-01 
.  939924E+00 
. 984111E+00 
.  992867E+00 
. 995970E400 
. 997413E400 
. 998198E400 
. 998671E+00 
. 998979E+00 
. 999189E+00 
. 999340E+00 
. 999452E+00 
. 999S37E+00 
. 999603E+00 
. 999655E+00 
. 999697E+00 
. 999732E+Q0 
. 999761E+00 
. 999785E+00 
. 999805E+00 
. 399322Et0C 
. 999837E+00 
. 999850E+00 
. 999861E+00 
. 999Q71E+Q0 
. 999880E+00 
. 999887E+00 
. 999694E+00 
. 9999Q0E+Q0 
. 999906E+00 
.  9999112+00 
. 999915E+0Q 
. 999919E+00 
. 999923E400 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 9999342+00 
. 999936E+0Q 
.999938E+00 
. 999940E+00 
. 999942E+00 
. 999944E+00 
. 99994SE+00 
. 999947E+00 
. 999948E+00 
. 999949E+00 
. 999950E+00 
. 99995 1E+00 
. 999952E+00 


400966E+01 
502849E+01 
671234E+01 
100258E+02 
-. 194057E+02 
877335E+02 
. 200011E+02 
. 102594E+02 
. 68715 1E+01 
. 516108E+01 
. 412995E+01 
. 34400 1E+01 
. 294542E+01 
. 257306E+01 
. 228223E+01 
. 204849E+01 
. 185626E+01 
. 169518E+01 
. 155806E+01 
. 143976E+01 
. 133650E+01 
. 124548E+01 
. 116454E+01 
. 109196E+Q1 
. 102646E+01 
. 5685452400 
. 912S60E+Q0 
. 862616E+00 
. 816516E+QQ 
.  773779E+00 
.  734001E+0Q 
. 696832E+00 
♦ 661977E+00 
. 6291822+00 
. S98231E+00 
.  568938E+00 
.541134E400 
. 514669E+00 
.  489424E+00 
.4652792+00 
.4421362+00 
.4199102+00 
. 398511E400 
. 377878E+Q0 
. 357935E+00 
. 338632E+00 
. 319908E+00 
. 301717E+00 
. 2B4017E+00 
. 266763E+00 
. 249920E+00 
. 2334S0E+00 
. 217325E+00 
. 201511E+00 
. 1859832+00 
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END  OF  RUN,  SYSTEM  it  1 


E.  MAGNITUDE,  PHASE  OF  THE  AUGMENTED  HR  NOTCH  FILTER:  P  =  3 

INPUT  DATA  FOR  SYSTEM  it  1 


INPUT  DATA  SOURCEFILE:  digfreq. tst 
DEGREE  OF  NUMERATOR  =  5 

DEGREE  OF  DENOMINATOR  -  5 

dsorce  =  F 

NUMBER  OF  FREQUENCY  POINTS  =  100  MAGNITUDE  OPTION  =  STD 
STARTING  VALUE  QF  THETA  =  . 0Q0Q00E+00 
INCREMENT  OF  THETA  =  . 280600E-01 

THE  NUMERATOR  COEFFICIENTS  b(0) ,b( 1). . . b(L)  ARE: 

.  9900E+00  3299E-Q2  -  1861E-01  . 9374E-02  -.6610E+00 

6Q47E+00 


THE  DENOMINATOR  COEFFICIENTS  c(0) ,c( 1). . . c(N)  ARE: 

. 1000E+01  .  0Q00E+QQ  .0000E+00  . OOOOE+OO  6524E+00 

*. S967E+00 


OUTPUT  DATA  FOR  SYSTEM  0  1 


THETA 

(RADIANS) 

. OOOOOOE+OO 
. 28060QE-01 
. 561200E-01 
.84180QE-01 
. 112240E+00 
.  140300E*00 
. 168360E+00 
. 196420E+00 
.  224480E+00 
. 252540E+00 


MAGNITUDE 


. 999937E+00 
. 999938E+00 
. 999938E+00 
. 999938E+00 
. 999938E+00 
. 999938E+00 
. 999937E+00 
. 999936E+00 
. 99993SE+00 
. 999934E+00 


PHASE 

(DEGREES) 

.  OOOOOOE+OO 
194494E-0I 
3889705-01 
584018E-01 
78028SE-01 
978157E-01 
U7817E+00 
138065E+00 
158610F-+00 
179497E+00 
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. 280600E+00 
.  308660E+00 
.  336720E+00 
. 364780E+00 
. 392840E400 
. 420900E+00 
. 448960E+00 
. 477020E+00 
. S0S080E+00 
. 533140E+00 
. 561200E+00 
. 589260E+00 
. 617320E+00 
. 645380E+00 
. 673440E+00 
. 701500E+00 
. 729560E+Q0 
. 757620E+00 
. 785660E+00 
*  813740E+00 
. 841800E+00 
. S69860E+00 
. 897920E+00 
. 92596QE+Q0 
. 954040E+00 
. 9S21QCE+QQ 
. 101O16E+Q1 
. 103822E+01 
. 106628E+01 
. 109434E+01 
.  H2240E+01 
.  11S046E+01 
. U7352E+01 
.  120638E401 
. 123464E+01 
.  126270E+01 
.  129076E+01 
. 131S82E+01 
. 1346BUE+01 
. 137494E+01 
.  140300E+01 
. 143106E+01 
. 145912E+01 
. 14S718E+01 
.  151S24E+01 
. 1S4330E+01 
. 1S7136E+OI 
. 159942E401 
, 162748E+01 
. 165554E+01 
. 168360E+O1 
. 171166S+01 
. 173972E+01 
. 176778E+01 
. 179S84E*01 
. 182390E+01 


.  999933E+00 
.  999931E+00 
. 999929E+00 
. 999927E+00 
. 999925E+00 
. 999923E+00 
. 999920E+00 
. 999917E+00 
. 999914E+00 
. 999910E+00 
. 999906E+00 
. 999902E+Q0 
. 999896E+00 
. 999891E+00 
. 999884E+00 
. 999876E+00 
. 999868E+00 
.  999858E-KJ0 
. 999847E+00 
. 99983SE+Q0 
. 999820E+00 
. 999803E+00 
. 999783E+00 
. 999759E+00 
. 999731E4Q0 
. 999697E+00 
. 999655E+00 
. 9996O3E+O0 
. 999S37E+00 
. 99945 3E+00 
.  999343&4-00 
. 999193E+00 
. 998985E+00 
.99868 1E+00 
. 998213E+0Q 
. 997440E+00 
. 996024E+00 
. 992994E+00 
. 984530E+00 
. 942356E+00 
. 381933E-01 
. 939924E+0O 
.984UOE+0O 
. 992867E+00 
>  99597OE+0O 
. 9974I3E400 
. 998198E+0O 
. 998671E400 
.998979E+0O 
. 999189E+0O 
. 999340E+OO 
. 9994S2E+00 
. 999537E400 
. 999603E+00 
. 9996S3E+00 
. 999697E40O 


200787E+00 
222538E+00 
244798E+00 
267635E+00 
291127E+00 
315344E+00 
340375E+00 
366312E+00 
393263E+00 
-.421337E+00 
45067 1E-H50 
-.461411E+00 
-.513727E+00 
547810E+00 
583876E+00 
622189E+00 
663036E+Q0 
706766E+00 
753794E+00 
804590E+00 
859744E+00 
*. 919960E+00 
986096E+0Q 
105921E+01 
114061E+01 
123199E+01 
133547E401 
14538BE+01 
1S9094E-HJ1 
17S175E+01 
194343E+01 
-.217624E+01 
246558E+01 
•>  263557E401 
*. 332630E+01 
*. 40O966E+OS 
•.502848E4O1 
-.671234E401 
1002S8E402 
».  19405  7E-HJ2 
*. 877273E+02 
.200011E+O2 
.  102394EKi 
. 687151E+C1 
.  516.108E+01 
. 41299SE+01 
. 34400 1E+01 
.294542E401 
.2S7306E+01 
. 228223E+01 
. 204848E401 
. 185626E+01 
. 169S18E+01 
. 1S5805E+01 
. 143975E+01 
. 133650E401 
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.  185196E+01 
. 188002E+01 
. 190808E+01 
. 193614E+01 
. 196420E+01 
. 199226E+01 
. 202032E+01 
. 204838E+01 
. 207644E+01 
. 2104S0E+01 
. 213256E+01 
, 216062E+01 
. 218868E+01 
. 221674E+01 
. 224480E+01 
. 227286E+01 
. 230092E+01 
. 232898E+01 
. 235704E+G1 
. 238510E+01 
. 241316E+01 
. 244122E+01 
. 246928E+01 
. 249734E+01 
. 2S2540E+01 
. 2S5346E+Q1 
. 258132E+01 
. 2609S6E+01 
. 263764E+01 
. 266570E+01 
. 269376E+01 
.272162E+01 
.274986E+01 
. 277794E+01 


.  999732E+00 
. 999761E+00 
. 99978SE+00 
. 999805E+00 
. 999822E+C0 
. 999837E+00 
. 999850E+00 
. 99986 1E+00 
. 99987 1E+00 
. 999880E+00 
. 999887E+00 
. 999894E+00 
. 999900E+00 
. 999906E+00 
. 999911E+00 
. 999915E+00 
.  999919B-H30 
. 999922E+00 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 999934E+00 
. 999936E+00 
. 999938E+00 
. 999940E+00 
. 999942E+0Q 
. 999943E+00 
.  999945E+0G 
•999946E+00 
. 999946E+00 
. 999949E+00 
. 999950E+00 
.  99995 lK-fOO 
. 999952E+00 


.  124548E+01 
.  116453E+01 
.  109196E+Q1 
.  102645E+01 
.  966937E+00 
. 912553E+00 
. 862608E+00 
. 816508E+00 
. 773770E+00 
. 733993E+00 
. 696823E+00 
. 661968E+GQ 
. 629177E+00 
. 598225E+00 
. 568932E+00 
. 541125E+O0 
, 514664E+00 
. 489416E+00 
. 465272E+00 
. 442129E+00 
. 419903E+00 
. 398506E+00 
. 377S72E400 
. 357931E+00 
. 338631E+00 
.  31991  Vi >'-00 
. 3Q1721E+GO 
.264018E+00 
. 26877 1E+00 
. 249936E+00 
- 233465E+O0 
. 217347E+00 
. 201S33E+00 
. 186002E+00 


END  OF  RUN,  SYSTEM  0  X 
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APPENDIX  B.  PROGRAM  AND  DATA  OF  DIFFERENCE  EQUATION 
A.  PROGRAM  FOR  DIFFERENCE  EQUATION  CALCULATION 

DIFFEQ.  FOR  VERSION:  2/03/88 


PURPOSE:  THIS  PROGRAM  COMPUTES  THE  ITERATIVE  SOLUTION  TO  A 
LINEAR,  TIME -INVARIANT  (LTI)  DIFFERENCE  EQUATION. 

THE  DIFFERENCE  EQUATION  MUST  BE  IN  THE  FORM: 
y(ns)  =  a(l)*y(ns-l)  +  ...  +  a(N)*y(ns-N)  + 

b(0)*x(ns)  +  b(l)*x(ns-l)  +  ...  +  b(L)*x(ns-L). 
THE  PROGRAM  CONSISTS  OF  A  MAIN  PROGRAM  AND  TWO 
SUBROUTINES.  SUBROUTINE  diffeq  IS  CALLED  BY  THE  MAIN 
PROGRAM  TO  ITERATIVELY  SOLVE  THE  DIFFERENCE  EQUATIONS  AND 
SUBROUTINE  xgen  ALLOWS  THE  USER  THE  OPTION  OF  GENERATING 
THE  INPUT  SEQUENCE  x( )  BY  WRITING  THE  APPROPRIATE 
EQUATIONS.  IF  THE  USER  ELECTS  TO  GENERATE  THE  SEQUENCE 
x( )  BY  USING  xgen  THEN  THE  PROGRAM  MUST  BE  COMPILED 
AGAIN  BEFORE  EXECUTION.  THE  USER  HAS  THE  OPTION  OF 
SELECTING  ONE  OF  TOO  OPERATING  MODES:  BATCH  OR  TEST. 

IN  BATCH-  MODE  THE  AMOUNT  OF  INTERFACE  WITH  THE  USER 
IS  MINIMIZED  AND  IT  IS  ASSUMED  THAT  THE  INPUT  DATA 
HAS  BEEN  STORED  IN  THE  DEFAULT  FILE  'DIFFEQ. IN'.  IN 
TEST  MODE  THE  USER  IS  PROMPTED  FOR  THE  NAME  OF  THE 
INPUT  FILE  OR  HAS  THE  OPTION  OF  PERFORMING  A  TEST  RUN 
USING  THE  INPUT  DATA  STORED  IN  THE  FILE  'DIFFEQ.  TST' . 

IT  IS  RECOMMENDED  THAT  FIRST-TIME  USERS  SELECT  THE 
TEST  MODE  AND  MAKE  A  TRIAL  RUN  WITH  THE  PRESTORED 
INPUT  DATA.  THE  TEST  MODE  ECHOES  PORTIONS  OF  THE 
INPUT  DATA  ONTO  THE  MONITOR  TO  ALLOW  VERIFICATION  OF 
ITS  ACCURACY.  BOTH  BATCH  AND  TEST  MODES  ALLOW  THE 
USER  TO  SOLVE  UP  TO  FOUR  DIFFERENCE  EQUATIONS  IN  A 
SINGLE  PROGRAM  EXECUTION.  THE  OUTPUT  OF  THE  PROGRAM 
'DIFFEQ. FOR'  IS  STORED  IN  THE  ARRAY  y().  THE  OUTPUT  IS 
STORED  IN  TABULAR  FORM  IN  THE  OUTPUT  FILE  'DIFFEQ. FOR' 

AND  IN  A  FORM  SUITABLE  FOR  PLOTTING  IN  THE  FILE 
'DIFFEQ. DAT' . 


i’f-V  ■■'(•Mclfin'tirMtMt'/r, 


INPUT 


THIS  PROGRAM  ASSUMES  THAT  EACH  DIFFERENCE  EQUATION  IS  IN  THE 
CANONICAL  FORM: 

y(ns)  «  a(l)*y(ns-l)  +  ...  +  a(N)*y(n$-N)  + 

b(0)*x(ns)  +  b(l)*x(ns-l)  +  ...  v  b(L)*x(ns-L) 

L  *  A  NON-NEGATIVE  INTEGER,  THE  NUMBER  OF  INPUT  DELAYS. 

N  =  A  NON-NEGATIVE  INTEGER,  THE  NUMBER  OF  OUTPUT  DELAYS. 


99 


UOOOOOOOOOOOOOOOOOOCJOOOOOOOOOOCJOOOOOOOOOOOOOOOOOCJOOC200CJCJ 


a(l). . .  a(N)  =  REAL  COEFFICIENTS  OF  THE  OUTPUT  TERMS. 
b(0)...b(L)  =  REAL  COEFFICIENTS  OF  THE  INPUT  TERMS. 

THE  INPUT  PARAMETERS  SHOULD  BE  STORED  IN  A  FILE  NAMED 
'DIFFEQ.IN'.  ALL  OF  THE  READ  STATEMENTS  USED  BY  THIS  PfOGRAM 
REQUIRE  FORMATTED  INPUT.  PARTICULAR  ATTENTION  SHOULD  BE  PAID 
TO  THE  FORMATS,  ESPECIALLY  THE  USE  OF  THE  DECIMAL  POINT  TO 
DENOTE  'REAL'  NUMBERS.  THE  INPUT  PARAMETERS  REQUIRED  BY  THE 
PROGRAM  ARE  LISTED  BELOW. 


NAME 

TYPE 

RANGE  (ARRAYS) 

RESTRICTIONS 

numsys 

INTEGER 

1 

<=  numsys  <=  4 

L 

INTEGER 

0  <=  L  <=  128 

N 

INTEGER 

0  <=  N  <=  128 

nstop 

INTEGER 

0 

<=  nstop  <-  999 

xsorce 

CHARACTER 

'F'  OR  'S' 

b(k) 

REAL 

0,1,  ... ,  L 

0  <=  L  <=  128 

a(k) 

REAL 

1,2,  ...  ,  N 

0  <=  N  <=  128 

y(k) 

REAL 

-N,  ....  -1 

1  <=  N  <*  128 

x(ns) 

REAL 

0,1,  ....  nstop 

0 

<=  nstop  <=  999 

WHERE: 


numsys  =  THE  NUMBt*  OF  SYSTEMS  TO  BE  EVALUATED. 

THIS  INTEGER  VALUE  MUST  OCCUR  AT  THE  TOP  OF  THE  INPUT 
FILE.  IT  DELINEATES  THE  NUMBER  OF  SYSTEMS  TO  BE  READ  BY 
THE  PROGRAM  AND  ANALYSED.  FOR  EACH  SYSTEM  ( 1. .. numsys) 
THE  PARAMETERS  BELOW  MUST  APPEAR  IN  THE  INPUT  FILE. 

L  =  AN  INTEGER  VALUE  THAT  SPECIFIES  THE  MAXIMUM  NUMBER  OF  DELAYS 
. IN  THE  INPUT  SEQUENCE. 

N  =  AN  INTEGER  VALUE  THAT  SPECIFIES  THE  MAXIMUM  NUMBER  OF  DELAYS 
IN  THE  OUTPUT  SEQUENCE. 

us top  a  AN  IirPEGER  VALUE  THAT  SPECIFIES  THE  LARGEST  TIME  INDEX 
(ns)  FOR  WHICH  THE  DIFFERENCE  EQUATION  IS  TO  BE  SOLVED. 

xsorce  «  A  CHARACTER  VALUE  OF  *F'  OR  'S'  DENOTING  WHETHER  THE 

INPUT  SEQUENCE  x()  IS  TO  BE  READ  FROM  THE  INPUT  FILE  (F) 
OR  TO  BE  GENERATED  (S)  USING  THE  SUBROUTINE  xgen.  THIS 
LATER  OPTION  IS  ATTRACTIVE  WHEN  nstop  IS  A  LARGE  NUMBER 
AND  THE  INPUT  SEQUENCE  x()  CAN  BE  READILY  DESCRIBED  BY  AN 
ANALYTICAL  EXPRESSION.  IF  xsorce  *  ’S'  THE  USER  MUST 
PROVIDE  THE  APPROPRIATE  FORTRAN  STATEMENTS  IN  THE  SPACE 
PROVIDED  IN  SUBROUTINE  xgen  AND  THE  PROGRAM  MUST  BE 
RECOMPILED  BEFORE  EXECUTION. 

b(k)  a  REAL  COEFFICIENTS  OF  THE  INPUT  SEQUENCE  x(ns*k)  IN  THE 
ORDER:  b(0),  b(l),  ....  b(L). 

a(k)  »  REAL  COEFFICIENTS  OF  THE  OUTPUT  SEQUENCE  y(ns-k)  IN  THE 

ORDER:  a(l),  a(2),  ....  a(N).  IF  N  «  0  THEN  THE  EQUATION 
IS  NON-RECURSIVE  AND  NO  a(k)  COEFFICIENTS  SHOULD  BE  IN 
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THE  INPUT  FILE. 


y(k)  =  THE  INITIAL  CONDITIONS  FOR  THE  OUTPUT  SEQUENCE  IF  THE 
DIFFERENCE  EQUATION  IS  RECURSIVE,  I.E. ,  N  >  0.  THIS 
PROGRAM  CALCULATES  THE  SOLUTION  TO  THE  DIFFERENCE 
EQUATION  FROM  ns  =  0  TO  ns  =  nstop  THEREFORE  THE  INITIAL 
CONDITIONS  y(-N)  TO  y(-l)  MUST  BE  PROVIDED  IN  THE  INPUT 
FILE  IN  THE  ORDER:  y(-N),  y(-N+l),  ....  y(-l).  IF  N  **  0 
THEN  THE  EQUATION  IS  NON-RECURSIVE  AND  NO  INITIAL 
CONDITIONS  SHOULD  BE  GIVEN  IN  THE  INPUT  FILE. 

x(ns)  a  THE  INPUT  SEQUENCE.  IF  xsorce  =  'F*  THEN  THE  INPUT 

SEQUENCE  x(0) ,  ....  x(nstop)  MUST  BE  PROVIDED  BY  THE  USER 
IN  THE  INPUT  FILE.  IF  xsorce  =  'S'  THEN  THE  USER  HAS 
ELECTED  TO  GENERATE  THE  INPUT  SEQUENCE  BY  PROVIDING  TIE 
APPROPRIATE  FORTRAN  STATEMENTS  IN  THE  SUBROUTINE  xgev. 

NOTE:  THE  INPUT  FORMAT  STATEMENTS  OCCUR  IN  THE  MAIN  PROGRAM 

FOLLOWING  THE  CAPTION:  ********  INPUT  FORMAT  ********* 
THE  FORM  OF  THE  INPUT  DATA  FILE  IS: 


LINE  # 


ENTRIES 


FORMAT 


1 

2 

NEXT  NB  LINES 
NEXT  NA  LINES 
NEXT  NY  LINES 
NEXT  NX  LINES 


numsys 

L,N,nstop,xsorce 
b(k),  k=0,l,...  ,L 
a(k),  k=l,. . .  ,N 
y(k),  k=  -N,. . .  ,-l 
x(ns),  ns»  0,..., nstop 


il 

13, til, i3,t21,i3,t31,al 
6f  10.  0 
6f  10.  0 
6f  10.  0 
6f 10. 0 


WHERE:  NB  «  1  +  (L/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 


NA  =  0  IF  N  =  0  OR 

NA  “  1  +  ((N-U/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 


NY  *  0  IF  N  »  0  OR 

NY  a  1  +  C(N*l)/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 
NX  a  0  IF  xsorce  =  ’S’  OR 

NX  «  1  +  (ns top/ 6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 
IF  xsorce  =  'F' 


*NOTE:  FOR  nufflsys  >  1  THE  FORMAT  OF  LINES  2. . .  IS  REPEATED. 

THE  FORM.  T  f 10. 0  USED  FOR  INPUT  DATA  PERMITS  THE  DECIMAL 
POINT  TO  BE  PLACED  ANYWHERE  IN  THE  FIELD  OF  10  COLUMNS 
AND  ALSO  ALLOWS  THE  EXPONENTIAL  FORMAT  TO  BE  USED  (EG. 
3146.2  «  3. 1462E+03). 


***************  ***  *********  w*  OUTPUT 


**<**********************i)i**** 


THE  INPUT  DATA  AS  WELL  AS  THE  OUTPUT  DATA  ARE  STORED  IN  TABULAR 
FORM  IN  THE  FILE  ' DIFFEQ. OUT* .  ADDITIONALLY,  THE  INPUT  SEQUENCE 
AND  THE  OUTPUT  SEQUENCE  ARE  WRITTEN  INTO  THE  FILE  'DIFFEQ. DAT'  TO 
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C  FACILITATE  PLOTTING  BY  A  SEPARATE.  USER  SUPPLIED  PROGRAM.  THE 
C  FORMAT  OF  THE  DATA  IN  'DIFFEQ.  DAT*  IS:  el2. 6,  2X,  el2. 6.  THE  FIRST 
C  ENTRY  CORRESPONDS  TO  THE  ORDINATE  VALUE,  AND  THE  SECOND  ENTRY,  THE 
C  ABSCISSA  VALUE.  ADDITIONAL  HEADER  INFORMATION  IS  WRITTEN  INTO 
C  'DIFFEQ. DAT1  TO  ALLOW  FOR  CONTROL  AND  LABELING  OF  EACH  PLOT. 

C 

C 

EXAMPLE  *****irk***ir****ickit*ititit*in^irk-k* 


THE  INPUT  PARAMETERS  FOR  THE  SYSTEM  DESCRIBED  BELOW  ARE  STORED  IN 
THE  SAMPLE  INPUT  FILE  'DIFFEQ. TST'  AND  CAN  BE  USED  FOR  A  TRIAL 
RUN  IN  TEST  MODE. 

DIFFERENCE  EQUATION: 

y(ns)  =  1.  2*y(ns-l)  +  1.5*x(ns) 

GOAL:  TO  OBTAIN  THE  SOLUTION  TO  THIS  DIFFERENCE  EQUATION  FOR 

ns  =  0  TO  ns  *  10,  GIVEN:  x(0)...x(10)  =  100.0  AND 
THE  INITIAL  ''ONDITION  y(-l)  =  25.0. 


THE  INPUT  FILE  IS: 

1 

000  001  010 

1.5 
1.2 
25.0 

100.0  100.0  100.0 

100.  0  100. 0  100. 0 


F 

1C0.0  100.0  100.0 

100. 0  100. 0 


THE  RESULTING  OUTPUT  FILE  'DIFFEQ. -UT'  IS: 

INPUT  DATA  FOR  PROBLEM  0  1 

PROBLEM  it  1  INPUT  DATA  SO'JRCtFILE:  DIFFEQ  TST 

THE  NUMBER  OF  INPUT  DELAYS:  L  =  0 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =  1 

THE  VALUE  OF  r.stop  IS:  10 

THE  COEFFICIENTS  b(0),  b(l),  ....  b(L) 

.  1SOOOOE4-C1 


THE  COEFFICIENTS  a(l),  ....  a(N>  ARE: 
. 120000E+01 


OUTPUT  DATA  FOR  PROBLEM  0  1 

ns  x(ns)  y(ns) 

-1  .  OOOOv.  DE+00  . 250000E+02 
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c 

0 

.  100000E+03 

.  180000E+03 

C 

1 

.  100000E+03 

.  366000E+03 

C 

2 

.  100000E+03 

.  589200E+03 

C 

3 

.  100000E+03 

.  857040E+03 

C 

4 

.  100000E+03 

.  117845E+04 

C 

5 

.  100000E+03 

.  156414E+04 

C 

6 

.  100000E+03 

.  20269 7E+04 

C 

7 

.  100000E+03 

.  258236E+04 

C 

8 

.  100000E+03 

.  324883E+04 

c 

9 

.  100000E+03 

.  404860E+04 

c 

10 

.  100000E+03 

.  500832E+04 

c 

C . END  OF  PROBLEM  #  1  . 

C 

C 

C  FOR  ILLUSTRATIVE  PURPOSES  THE  INPUT  SEQUENCE  x()  COULD  HAVE  BEEN 
C  GENERATED  BY  SPECIFYING  xsorce  =  'S’  AND  WRITING  THE  APPROPRIATE 
C  FORTRAN  STATEMENTS  INTO  SUBROUTINE  xgen.  THE  STATEMENTS  THAT 
C  COULD  BE  USED  TO  ACCOMPLISH  THIS  ARE  WRITTEN  INTO  THE  SUBROUTINE 
C  BUT  ARE  ’COMMENTED  OUT*. 

C 

c 

PROGRAM 


character  infile,v12,  mode*l,  xsorce*l 

real  a( 1:128),  b(0:128),  y( -128: 850),  x( -128: 850),  ii 

C  PROMPT  USER  FOR  MODE:  BATCH  OR  TEST. 

write(*, 1115) 
read(*,1117)  mode 

if  ((mode.  eq.  *  Y* ).  or.  (mode.  eq.  'y  * ))  then, 
mode  =  *Y' 
write(*, 1118) 
read(*,lll9)  infile 
else 

infile  =  ' DIFFEQ. IN' 
end  if 

C  UNIT=1  DEFINED  AS  INPUT  FILE.  UNITS«2,3  DEFINED  AS  OUTPUT  FILES. 

open(unit=l,flle-infile,status='old' ,iostat«ierr,err=999) 
open( unit-2 .files' DIFFEQ. OUT' ) 
open(unit“3 .file*' DIFFEQ. DAT ' ) 

C  READ  INPUT  PARAMETERS  AND  CONDUCT  ERROR  CHECKS. 

read( 1,1000)  numsys 

if((numsys. le.  0).  or.  (numsys. gt. 4))  then 
write(*,1126)  numsys 

stop  'The  allowed  values  for  numsys  are  1  <=  numsys  <=  4. ' 
end  if 

numplts  =  2*numsys 
write( 3,2000)  numplts 
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do  10  nprob=l,  numsys 

read( 1,1001)  L,  N,  ns top,  xsorce 

if ( C L.  It.  0).  or.  (L.  gt.  128))  then 
write(*,1124)  nprob,  'L  ,  L 

stop  'The  allowed  values  for  "L"  are:  o  <=  L  <=  128.' 
endif 

if((N.  It.  0).  or.  (N.  gt.  128))  then 
write(*,1124)  nprob,  'Nl ,  N 

stop  'The  allowed  values  for  "N"  are*.  0  <=  N  <=  128.' 
endif 

if((nstop.  It.  0).or.  (nstop.  gt.  850))  then 
write(*, 1127)  nprob,  nstop 

stop  'The  allowed  values  for  nstop  are:  0  <=  nstop  <=  850. ' 
endif 

if((xsorce.  eq.  'F').or.  (xsorce.  eq.  'f' ))  then 
xsorce  =>  'F* 

elseif( (xsorce. eq.  'S').or.  (xsorce. eq. 's' ))  then 
xsorce  =  'S' 
else 

write(*,1126)  nprob,  xsorce 

stop  'The  allowed  values  for  "xsorce"  are:  "F"  or  "S".  ' 
endif 

C  INITIALIZE  EACH  ARRAY  TO  ZERO  BEFORE  EACH  RUN. 

data  a/ 128*0.0/,  b/ 129*0.0/ 
data  y/979*0. 0/,  x/979*0.0/ 

C  READ  THE  COEFFICIENTS  b(),  a()  AND  THE  INITIAL  CONDITIONS 
C  y(-N)...y(-l). 

read( 1,1002)  (b(k),  k=0,L) 
if(N. gt. 0)  then 
icstart  =  -N 

read(l, 1002)  (a(k),  k=l,N) 
read( 1,1002)  (y(k),  k“icstart,-l) 
endif 

C  FOR  xsorce  =  *F'  READ  THE  ARRAY  x()  FROM  THE  INPUT  FILE. 

C  FOR  xsorce  «  'S'  CALL  xgen  TO  GENERATE  THE  ARRAY  x(). 

if (xsorce.  eq.  'F* )  then 
read(l, 1002)  (x(k),  k*0, nstop) 
else 

call  xgen(x, ns top, nprob) 
endif 

C  FOR  TEST  MODE  ECHO  INPUT  PARAMETERS  ONTO  MONITOR  (UNIT  =  *). 
if(mode. eq. 'Y')  then 
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write(*,1007) 

write(*,1120)  nprob,  infile 
write(*,1110)  'INPUT'  V,  L 
write(*,1110)  'OITPUT1,  'Nf,  N 
write(*,1112)  nstcp 
write(*, 1004) 

write(*,1005)  (b(k),k=0,L) 
if(N. eq.  0)  then 
write(*,1131) 
cXs& 

write(*,1006) 

write(*, 1005)  (a(k),k=l,N) 
endif 

write(*,1123)  nprob 

pause  'END  OF  RUN,  STRIKE  <CR>  WHEN  READY  TO  CONTINUE. ' 
endif 

C  WRITE  INPUT  DATA  INTO  FILE:  DIFFEQ. OUT. 

write(2,1008)  'INPUT',  nprob 
write(2, 1120)  nprob,  infile 
write(2,1110)  'INPUT'  'L*.  L 
write(2, 1110)  'OUTPUT' ,  'N  ,  N 
write(2,1112)  nstop 
write(2, 1004) 

write(2, 1005)  (b(k),k=0,L) 
if ( N. eq. 0)  then 
write(2,1131) 
else 

write(2,1006) 

write(2,1005)  (a(k),k=l,N) 
endif 

C  WRITE  THE  INPUT  SEQUENCE  INTO  FILE:  DIFFEQ. DAT. 

write(3,2001)  nstop  +  1 
write(3,*)  'INPUT  SEQUENCE  x(ns)’ 
write(  3,*v)  'SAMPLE  0  (ns)' 
write(3,*)  'x(ns)' 
do  54  ns=0,  nstop 
ii  E  ns 

writo(3 ,2010)  ii,  x(ns) 

54  continue 

C  CALL  diffeq  TO  COMPUTE  THE  SOLUTION  TO  THE  DIFFERENCE  EQUATION. 

call  diffeq(N,L,a,b,x,y,nstop) 

C  WRITE  RESULTS  INTO  FILE:  DIFFEQ,  DAT. 

write(3,2001)  N  +  nstop  +  l 
writ  <:.*)  'OUTPUT  SEQUENCE  Y(n s)' 
writ  *,*)  'SAMPLE  0  (ns)’ 
writ  vj,*)  'y(ns) ' 
do  55  ns=  -N,  nstop 
ii  -  ns 
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write(3,2010)  ii,  y(ns) 

55  continue 

C  WRITE  RESULTS  INTO  FILE:  DIFFEQ. OUT. 

write(2,1008)  'OUTPUT',  nprob 
write(2, 1129) 
do  102  ns=  -N,  nstop 

write(2,1130)  ns,  x(ns),  y(ns) 

102  continue 

wxite(2,1123)  nprob 

10  continue 

write(*,1121)  . 

999  close(unit=l) 
close(unit»2) 
close(unit=3) 


if(ierr. gt.  0)  then 
write(*,1116)  ierr 
endif 

q********  INPUT  FORMAT  ****** *** 

1000  format( il) 

1001  format( i3,tll , i3,t21,i3,t31,al) 

1002  formatebflO.  0) 

1004  format(t4, 'THE  COEFFICIENTS  b(0),  bei),  ....  beL)  ARE:',/) 

1005  formate6Clx,el2.6)) 

1006  formate //,t4, 'THE  COEFFICIENTS  a(l) . aeN)  ARE:',/) 

1007  formate//////////) 

1008  format(///.tl6,a6» '  DATA  FOR  PROBLEM  //  ’  il,//) 

1110  format( t4,cTHE  NUMBER  OF  *,a6,'  DELAYS:  1  ,al,  »  ’,i3) 

1112  format(t4,’THE  VALUE  OF  natop  IS:  \i3) 

1115  Oformate lx, 'DO  YOU  WISH  TO  RUN  THIS  PROGRAM  IN  TEST', 

1'  MODE  ?  CY/N)  <CR>  :  ',  ,) 

1116  0format( /// , lx , 1  ERROR  OPENING  INPUT  FILE,  PROGRAM  TERMINATED.  ' , 
1//, lx, 'ERROR  CODE:  ',i4 ,/////) 

1117  formateal) 

1118  Oformate///// » lx, 'TYPE  THE  NAME  OF  YOUR  DATA  FILE  FOLLOWED'. 

1'  BY  <CR>. ' ,/,'  IF  YOU  DESIRE  TO  MAKE  A  TEST  RUN  USING  THE1, 

2'  SAMPLE  DATA  ALREADY  STORED',/,'  IN  THE  FILE:  DIFFEQ. TST’, 

3  '  TYPE:  DIFFEQ. TST  <CR>' ,/,'  FILENAME:  ',  ,) 

1119  formate al2) 

1120  format(///.t4, 'PROBLEM  //  '  .il, '  INPUT  DATA  SOURCEF I LE:  \al2) 

1121  Oformate//, '  TABULAR  OUTPUT  DATA  IS  STORED  IN  FILE:  DIFFEQ. OUT. ’ 
1,/,'  PLOTTING  DATA  IS  STORED  IN  FILE:  DIFFEQ.  DAT. ' ) 

1122  format( i3) 

1123  formate/, lx, loe '-'),'  END  OF  PROBLEM  #' ,12, 2x,16(* -'),//) 

1124  Oform&t(//,'  For  problem  #',i2,'  the  value  for  ',al,'  is:  ',i3, 
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1*.  This  value  is  not  allowed.’) 

1126  format(//,'  numsys  =  ',i4.'.  This  value  is  not  allowed.’) 

1127  Oformat(//,'  For  problem  #  ,12,’  the  value  for  ’’nstop’’  is: 
li3,'.  This  value  is  not  allowed.') 

1128  Oformat(//,'  For  problem  #’,i2,'  the  value  for  ’ ’xsorce' '  is: 
lal,’.  This  value  is  not  allowed.’) 

1129  format( t6, 'ns' ,tl6, ’xens) ’ ,t35, 'y(ns)' ) 

1130  formate t4,i4,tll,el4.  6,t30,el4.  6) 

1131  formate//  THIS  SYSTEM  IS  NON-RECURSIVE,  I.E. ,  N  =  0. ') 

2000  format( il) 

2001  formate  13) 

2010  format(el2. 6,2x,el2. 6) 

end 


SUBROUTINE:  diffeq 


C  PURPOSE:  THIS  SUBROUTINE  COMPUTES  THE  SOLUTION  TO  A  DIFFERENCE 
C  EQUATION.  ALL  PARAMETERS  DESCRIBING  THE  EQUATION,  AND 

C  THE  INPUT  AND  OUTPUT  SEQUENCES  AND  y()  ARE  PASSED 

C  TO  THE  SUBROUTINE  BY  THE  MAIN  PROGRAM. 

subroutine  diffeq(N,L,a,b,x,y,nstop) 

real  xe-128:nrtop) ,  ye -128:  nstop) ,  a(l:N)»  b(0:L) 


do  500  ns=0,  nstop 
y(ns)  =  0.  0 
do  501  k=0.  max(K, L) 

yens)  =  y(ns)  +  a(k)*yens-k)  +  b(k)*x(ns-k) 
501  continue 
500  continue 

return 

end 


C  SUBROUTINE:  xgen 


C  PURPOSE:  THIS  SUBROUTINE  ALLOWS  THE  USER  TO  GENERATE  VALUES  FOR 
C  THE  ARRAY  x().  IF  xsorce  *  'S*  THE  MAIN  PROGRAM  WILL 
C  CALL  THIS  SUBROUTINE.  IF  xsorce  =  *F*  THIS  SUBROUTINE 
C  WILL  NOT  BE  CALLED  BY  THE  MAIN  PROGRAM. 


subroutine  xgen(x,nstop,nprob) 
real  x( -128: nstop) 


C  DEVELOP  THE  ALGORITHM  FOR  GENERATING  VALUES  OF  x( )  IN  THIS  SPACE. 
C  THE  STATEMENTS  TYPED  IN  MUST  FOLLOW  STANDARD  FORTRAN  77  RULES  AND 
C  MAY  USE  FORTRAN  77  INTRINSIC  FUNCTIONS  SUCH  AS:  SINC),  COSe),  ... 
C  NOTE  THAT  THE  VALUE  nprob  CAN  BE  USED  IN  A  LOGICAL  ? IF*  STATEMENT 
C  TO  MATCH  THE  GENERATING  FUNCTIONS  TO  THE  CORRESPONDING  SYSTEM 
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EQUATION  READ  FROM  THE  INPUT  FILE  IF  MORE  THAN  ONE  SYSTEM  OF 
EQUATIONS  EXIST.  AN  EXAMPLE  OF  AN  ALGORITHM  GENERATING  VALUES 
FOR  x()  IS: 


**  EXAMPLE  *** 

if(nprob. eq. 1)  then 
do  1  k=0,  nstop 
x(k)  =  100.0 
1  continue 
end  if 

x(0)  »  1.0 
do  1  k=l, nstop 
x(k)  =  0.0 
1  continue 


return 

end 


B.  INPUT,  OUTPUT  DATA  OF  SECOND  ORDER  HR  NOTCH  FILTER:  P  =  0 


INPUT  DATA  FOR  PROBLEM  #  1 


PROBLEM  #  1  INPUT  DATA  S0URCEFILE:  diffeq.tst 

THE  NUMBER  OF  INPUT  DELAYS:  L  **  2 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  *  2 

THE  VALUE  OF  nstop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  ....  b(L)  ARE: 

. 990000E+00  329930E+00  .  990000E+00 


THE  COEFFICIENTS  a(l),  ...,  a(N)  ARE: 
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329930E+00  980100E+00 


OUTPUT  DATA  FOR  PROBLEM  #  1 


ns 

x(ns) 

y(ns) 

-2 

. 000000E+00 

. OOOOOOE+OO 

-1 

. OOOOOOE+OO 

. OOOOOOE+OO 

0 

. 100000E+01 

. 990000E+00 

1 

. 000000E+00 

329930E-02 

2 

. OOOOOOE+OO 

. 186125E-01 

3 

. OOOOOOE+OO 

. 937446E-02 

4 

. OOOOOOE+OO 

151492E-G1 

5 

.  OOOOOOE+OO 

141861E-01 

6 

. OOOOOOE+OO 

. 101673E-01 

7 

. OOOOOOE+OO 

. 172583E-01 

8 

. OOOOOOE+QO 

427095E-02 

9 

. OOOOOOE+OO 

183239E-01 

10 

. OOOOOOE+OO 

185966E-02 

11 

. OOOOOOE+OO 

. 173457E-01 

12 

. OOOOOOE+OO 

. 754553E-02 

13 

. OOOOOOE+OO 

145111E-01 

14 

.  OOOOOOE+OO 

121830E-01 

15 

. OOOOOOE+OO 

. 102028E-Q1 

16 

. OOOOOOE+OO 

. 153068E-01 

17 

. OOOOOOE+OO 

494956E-02 

13 

. OOOOOOE+OO 

166352E-01 

19 

. OOOOOOE+OO 

637380E-03 

20 

. OOOOOOE+OO 

. 160938E-01 

21 

. OOOOOOE+OO 

. 593453E-02 

22 

. OOOOOOE+OO 

138156E-01 

23 

. OOOOOOE+OO 

103746E-01 

24 

. OOOOOOE+OO 

. 101178E-01 

25 

. OOOOOOE+OO 

. 135063E-01 

26 

. OOOOOOE+OO 

546028E-02 

27 

. OOOOOOE+OO 

150390E-01 

28 

. OOOOOOE+OO 

. 389786E-03 

29 

.OOOOOOE+OO 

. 148664E-01 

30 

.  OOOOOOE+OO 

. 4S2349E-02 

31 

.  OOOOOOE+OO 

130801E-01 

32 

.  OOOOOOE+OO 

874898E-02 

33 

. OOOOOOE+OO 

. 993321E-02 

34 

. OOOOOOE+OO 

. 118521E-01 

35 

.  OOOOOOE+OO 

-.582S16E-02 

36 

.  OOOOOCE+OO 

135382E-01 

37 

.  ooooo;' i+oo 

. 124259E-02 

38 

.  OOOOOOE+OO 

. 136787E-01 

39 

.  OOOOOOE+OO 

. 329516E-02 

40 

. OOOOOOE+OO 

123194E-0i 

41 

.  OOOOOOE+OO 

-.7294UE-02 

42 

.  OOOOOOE+OO 

. 966765E-02 

43 

. OOOOOOE+OO 

. 103366E-01 

44 

.  OOOOOOE+OO 

606425E-02 

45 

. OOOOOOE+OO 

121336E-01 
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46 

. OOOOOOE+OO 

. 194032E-02 

47 

. OOOOOOE+OO 

. 125324E-01 

48 

. OOOOOOE+OO 

. 223309E-02 

49 

. OOOOOOE+OO 

115462E-01 

50 

. OOOOOOE+OO 

599809E-02 

51 

. OOOOOOE+OO 

. 933748E-02 

52 

. OOOOOOE+OO 

. 895944E-02 

53 

. OOOOOOE+OO 

619567E-02 

54 

. OOOOOOE+OO 

108253E-01 

55 

. OOOOOOE+OO 

. 250079E-02 

56 

. OOOOOOE+OO 

. 114349E-01 

57 

. OOOOOOE+OO 

. 132171E-02 

58 

. OOOOOOE+OO 

107713E-01 

59 

. OOOOOOE+OO 

484919E-02 

60 

. OOOOOOE+OO 

. 895708E-02 

61 

. OOOOOOE+OO 

. 770790E-02 

62 

. OOOOOOE+OO 

623577E-02 

63 

.  OOOOOOE+OO 

961188E-02 

64 

.  OOOOOOE+OO 

, 294043E-02 

65 

.  OOOOOOE+OO 

. 103907E-01 

66 

.  OOOOOOE+OO 

. 546304E-03 

67 

.  OOOOOOE+OO 

100037E-01 

68 

.  OOOOOOE+OO 

383596E-Q2 

69 

.  OOOOOOE+OO 

. 853905E-02 

70 

.  OOOOOOE+OO 

. 657691E-02 

71 

.  OOOOOOE+OO 

-.619920E-02 

72 

.  OOOOOOE+OO 

849133E-02 

73 

.  OOOOOOE+OO 

. 327429E-02 

74 

.  OOOOOOE+OO 

. 940264E-02 

75 

.  OOOOOOE+OO 

106918E-03 

76 

.  OOOOOOE+OO 

*. 925080E-02 

77 

. OOOOOOE+OO 

294733E-02 

78 

. OOOOOOE+OO 

. 80943QE-02 

79 

. OOOOOOE+OO 

. 555923E-02 

80 

. OOOOOOE+OO 

609907E-02 

81 

.  OOOOOOE+OO 

746086E-02 

82 

.  OOOOOOE+OO 

. 351613E-02 

63 

.  OOOOOOE+OO 

. 847247E-02 

S4 

.  OOOOOOE+OO 

-.650839E-03 

65 

,  OOOOOOE+OO 

851860E-02 

86 

.  OOOOOOE+OO 

217265E-02 

87 

.  OOOOOOE+OO 

. 763226E-02 

88 

.  OOOOOOE+OO 

. 464753E-02 

89 

.  OOOOOOE+OO 

*. 594701E-02 

90 

.  OOOOOOE+OO 

-.651714E-02 

91 

.  OOOOOOE+OO 

. 367847E-02 

92 

. OOOOOOE+OO 

. 760109E-02 

93 

.  OOOOOOE+OO 

109744E-02 

94 

.  OOOOOOE+OO 

781190E-02 

95 

.  OOOOOOE+OO 

150178E-02 

96 

.  OOOOOOE+OO 

. 716096E-02 

97 

.  OOOOOOE+OO 

. 38345  IE-02 

98 

.  OOOOOOE+OO 

575334E-02 

99 

.  OOOOOOE+OO 

565640E-02 

100 

.  OOOOOOE+OO 

.377263E-02 

101 

.  OOOOOOE+OO 

. 678855E-02 

no 


102 

. OOOOOOE+OO 

145781E-02 

103 

.  OOOOOOE+OO 

713443E-02 

104 

. OOOOOOE+OO 

925062E-03 

105 

. OOOOOOE+OO 

.668725E-02 

106 

.  OOOOOOE+OO 

.  3U298E-02 

107 

. OOOOOOE+OO 

552711E-02 

108 

. OOOOOOE+OO 

487459E-02 

109 

.  OOOOOOE+OO 

. 330885E-02 

110 

.  OOOOOOE+OO 

. 603424E-02 

111 

. OOOOOOE+OO 

174217E-02 

112 

. OOOOOOE+OO 

648895E-02 

113 

.  OOOOOOE+OO 

433395E-03 

114 

.  OOOOOOE+OO 

. 621683E-02 

115 

. OOOOOOE+OO 

. 247589E-02 

116 

. OOOOOOE+OO 

527624E-02 

117 

. OOOOOOE+OO 

416741E-02 

118 

. OOOOOOE+OO 

♦ 379629E-02 

119 

. OOOOOOE+OO 

. 533699E-02 

120 

. OOOOOOE+OO 

195991E-02 

121 

. OOOOOOE+OO 

587742E-02 

122 

. OOOOOOE+OO 

182246E-04 

123 

. OOOOOOE+OO 

. 575444E-02 

124 

. OOCOOOE+OO 

. 191643E-02 

125 

. COQQOOE+OO 

«. 500764E-02 

126 

.  OOOOOOE+OO 

353046E-02 

127 

.  OOOOOOE+OO 

. 374319E-Q2 

128 

.  OOOOOOE+OO 

. 469519E-02 

129 

.  OOOOOOE+OO 

211961E-02 

130 

.  OOOOOOE+OO 

*. 530108E-02 

131 

.  OOOOOOE+OO 

*  328445E-03 

132 

.  OOOOOOE+OO 

. 530396E-02 

133 

.  OOOOOOE+OO 

.  142802E-02 

134 

.  OOOOOOE+OO 

*. 472726E-02 

135 

.  OOOOOOE+OO 

295927E-02 

136 

.  OQOOOQE+OO 

. 365683E*02 

137 

.  OOOOOOE+OO 

•410688E-02 

138 

.  OOOOOOE+OO 

222908E-02 

139 

. OOOOOOE+OO 

476059E-02 

140 

*  OOOOOOE+OO 

•614057E-03 

141 

. OOOOOOE+OO 

. 466845E-02 

142 

. OOOOOOE+OO 

. 10044  IE-02 

143 

.  OOOOOOE+OO 

444019E-02 

144 

. OOOOOOE+OO 

244937E-02 

145 

. OOOOOOE+OO 

. 3S4370E-02 

146 

. OOOOOOE+OO 

.356981E-02 

147 

. OOOOOOE+OO 

229540E-02 

148 

. CQOOOOE+OO 

42560°E-02 

149 

. OUOOOOE+OO 

. 84550SE-03 

ISO 

.  OOOOOOE+OO 

. 445035E*02 

1S1 

. OOOOOOE+OO 

. 639621E-03 

152 

. OOOOOOE+OO 

-.415076E-02 

153 

. OOOOOOE+OO 

199635E-02 

154 

. OOOOOOE+OO 

. 340950E-02 

155 

. OOOOOOE+OO 

. 308152E-02 

1S6 

. OOOOOOE+OO 

232497E-02 

157 

. OOOOOOE+OO 

378727E-02 

ill 


158 

. OOOOQOE+OO 

. 102916E-02 

159 

. 000000E+00 

. 405146E-02 

160 

. OOOOOOE+OO 

. 328015E-03 

161 

. OOOOOOE+OO 

386261E-02 

162 

. OOOOOOE+OO 

159588E-02 

163 

. OOOOOOE+OO 

. 325922E-02 

164 

. OOOOOOE+OO 

. 263944E-02 

165 

.  OOOOOOE+OO 

*. 232353E-02 

166 

.  OOOOOOE+OO 

33535 IE-02 

167 

.  OOOOOOE+OO 

. 117087E-02 

168 

.  OOOOOOE+OO 

.367308E-02 

169 

.  OOOOOOE+OO 

. 642924E-04 

170 

.  OOOOOOE+OO 

357878E-02 

171 

.  OOOCOOE+QO 

124376E-02 

172 

.  OOOOOOE+OO 

. 309721E-02 

173 

.  OOOOOOE+OO 

. 224087E-02 

174 

.  OOOOOOE+OO 

229624E-02 

175 

.  OOOOOOE+OO 

295387E-02 

176 

,  OOOOOOE+OO 

. 127597E-02 

177 

.  OOOOOOE+OO 

. 331607E-02 

178 

.  OOOOOOE+OO 

156510E-03 

179 

.  OOOOOOE+OO 

330172E-C2 

180 

.  OOOOOOE+OO 

-.93S942E-03 

181 

.  OOOOOOE+OO 

. 292722E-02 

182 

.  OOOOOOE+OO 

. 18S310E-02 

183 

.  OOOOOOE+OO 

-.224766E-02 

164 

.  OOOOOOE+OO 

258720E-02 

185 

.  OOOOOOE+OO 

. 134936E-02 

186 

.  OOOOOOE+OO 

.298091E-02 

187 

.  OOOOOOE+OO 

339014E-03 

188 

.  OOOOOOE+OO 

303344E-02 

189 

.  OOOOOOE+OO 

668554E-03 

190 

.  OOOOOOE+OO 

. 275250E»02 

191 

.  OOOOOOE+OO 

. 1S6336S-02 

192 

.  OOOOOQE+QQ 

*. 218192E‘02 

193 

.  OOOOOOE+OO 

-.225215E-02 

194 

.  OQOOOOE+OO 

. 139S44E-02 

195 

.  OOOOOOE+OO 

.266773E-02 

196 

.  OOOOOOE+OO 

-.487S10E-03 

197 

.  OOOOOOE+OO 

-.  277549E-02 

198 

.  OOOOOOE+OO 

-.437908E-03 

199 

.  OOOOOOE+OO 

. 257578E“02 

200 

OOOOOOE+OO 

. 127902E-02 

201 

.  OOOOOOE+OO 

2102S3E-02 

202 

.  OOOOOOE+OO 

194725E-02 

203 

.  OOOOOOE+OO 

. 141623E-02 

204 

.  OOOOOOE+OO 

. 237642E-02 

205 

. OOOOOOE+OO 

605957E-03 

206 

.  OOOOOOE+OO 

252905E-02 

207 

.  ooooooi- -mo 

240512E-03 

208 

.OOOOOOE+OO 

. 239937E-02 

209 

.  OOOOOOE+OO 

. 102735E-02 

210 

.  OOOOOOE+OO 

201267E-02 

211 

.  OOOOOOE+OO 

*. 167095E-02 

212 

.  OOOOOOE+OO 

. 142132E-02 

213 

.  OOOOOOE+OO 

. 210663E-02 
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214 

. 000000E+00 

697996”  03 

215 

. OOOOOOE+OO 

229500L  62 

216 

. OOOOOOE+OO 

730823E-04 

217 

. OOOOOOE+OO 

. 222522E-02 

218 

. OOOOOOE+OO 

. 805795E-03 

219 

. OOOOOOE+OO 

191508E-02 

220 

. OOOOOOE+OO 

142160E-02 

221 

.  OOOOOOE+OO 

. 140794E-02 

222 

.  OOOOOOE+OO 

. 185783E-02 

223 

. OOOOOOE+OO 

766969E-03 

224 

. OOOOOOE+OO 

207391E-02 

225 

. OOOOOOE+OO 

. 674610E-04 

226 

. OOOOOOE+OO 

. 205490E-02 

227 

. OOOOOOE+OO 

.611853E-03 

228 

. OOOOOOE+OO 

181214E-02 

229 

. OOOOOOE+OO 

119756E-02 

230 

. OOOOOOE+OO 

. 138096E-02 

231 

. OOOOOOE+OO 

. 162935E-02 

232 

. OOOOOOE+OO 

815913E-03 

233 

. OOOOOOE+OO 

186612E-02 

234 

. OOOOOOE+OO 

. 183989E-03 

235 

. OOOOOOE+OO 

. 188968E-02 

236 

. OOOOOOE+OO 

. 443136E-03 

237 

. OOOOOOE  vO 

170587Z-02 

238 

. OOOOOOE. 00 

997137E-03 

239 

. OOOOOOE+OO 

. 134294E-02 

240 

. OOOOOOE+OO 

.  142.037E-02 

241 

. OOOOOOE+OO 

847595E-03 

242 

. OOOOOOE+OO 

167175E-02 

243 

.  OOOOOOE+.,0 

. 279167E-03 

244 

. OOOOOOE+OO 

. 173059E-02 

245 

. OOOOOOE+OO 

. 297362E-03 

246 

. OOOOOOE+OO 

159804E-02 

247 

. OOOOOOE+OO 

818687E-03 

248 

.  OOOOOOE+OO 

. 129613E-02 

249 

.  OOOOOOE+OO 

. 123003E-U2 

250 

.  OOOOOOE+OO 

864516E-03 

251 

.  OOOOOOE+OO 

149078E-02 

252 

.  OOOOOOE+OO 

. 355459E-03 

253 

.  OOOOOOE+OO 

. 157839E -02 

254 

.  OOOOOOE+OO 

. 172373E-03 

255 

.  OOOOOOE+OO 

149011E-02 

256 

.  OOOOOOE+OO 

660574E-03 

257 

.  OOOOOOE+OO 

. 124251E-02 

256 

.  OOOOOOE+OO 

. 105737E-02 

259 

.OOOOOOE  1-00 

868928E-03 

260 

.  OOOOOOE+OO 

132301E-02 

261 

. OOOOOOF+OO 

. 415134E-Q3 

262 

. OOOOOOE+OO 

. 143365E-02 

263 

.  OOOOOOE+OO 

. 661320E-04 

264 

.  OOOOOOE+OO 

138330E-02 

265 

.  OOOOOOE+OO 

521209E-03 

266 

.  OOOOOOE+OO 

. 118381E-02 

267 

. OOOOOOE+OO 

. 901413E-03 

260 

.  OOOOOOE+OO 

662852E-03 

269 

.  OOGOOOE+OO 

116816E-02 

113 


fe£as3*at«2tei^& sft& 


-JS&KS 


a&iy&aBfefc i 


270 

. OOOOOOE+OO 

. 460272E-03 

271 

. OOOOOOE+OO 

. 129677E-02 

272 

.  OOOOOOE+OO 

232701E-04 

273 

.  OOOOOOE+OO 

127864E-02 

274 

.  OOOOOOE+OO 

399054E-03 

275 

.  OOOOOOE+OO 

. 112153E-02 

276 

.  OOOOOOE+OO 

. 761140E-03 

277 

.  OOOOOOE+OO 

848092E-03 

278 

.  OOOOOOE+OO 

102580E-02 

279 

.  OOOOOOE+OO 

. 492771E-03 

280 

.  OOOOOOE+OO 

. 116797E-02 

281 

.  OOOOOOE+OO 

976162E-04 

282 

.  OOOOOOE+OO 

117693E-02 

283 

.  OOOOOOE+OO 

292632E-03 

284 

. OOOOOOE+OO 

. 105697E-02 

285 

. OOOOOOE+OO 

. 635534E-03 

286 

. OOOOOOE+OO 

826250E-03 

287 

. OOOOOOE+OO 

895491E-03 

288 

.  ooooocz+oo 

. 514358E-03 

289 

.  OOOOOOE+OO 

. 104737E-02 

290 

.  OOOOOOE+OO 

158563E-03 

291 

.  Onnnn0ii+00 

107885E-02 

292 

. OOOOOOE+OO 

200536E-03 

293 

. OOOOOOE+OO 

. 991213E-03 

294 

. OOOOOOE+OO 

. 523576E-03 

295 

.  OOOOOOE+OO 

798744E-03 

296 

. OOOOOOE+OO 

776687E-03 

297 

. OOOOOOE+OO 

. 526597E-03 

298 

. OOOOOOE+OO 

. 934971E-03 

299 

. OOOOOOE+OO 

-.207643E-03 

300 

. OOOOOOE+OO 

984873E-03 

301 

.  OOOOOOE+OO 

121428E-03 

302 

.  OOOOOOE+OO 

. 925211E-03 

303 

.  OOOOOOE+OO 

. 424267E-03 

304 

. OOOOOOE+OO 

766821E-03 

305 

. OOOOOOE+OO 

668821E-03 

306 

. OOOOOOE+OO 

. 530897E-03 

307 

. OOOOOOE+OO 

.  830670E-03 

308 

.  OOOOOOE+OO 

246269E-03 

309 

.  OOOOOOE+OO 

895391E-03 

310 

. OOOOOOE+OO 

540483E-04 

311 

. OOOOOOE+OO 

. 859741E-03 

312 

. OOOOOOE+OO 

. 336627E-03 

313 

. OOOOOOE+OO 

731569E-03 

314 

. OOOOOOE+OO 

571295E-03 

315 

. OOOOOOE+OO 

. 528523E-03 

316 

. OOOOOOE+OO 

. 734302E-03 

317 

. OOOOOOE+OO 

275737E-03 

318 

. OOOOOOE+OO 

810663E-03 

319 

. OOOOOOE+OO 

. 278823E-05 

320 

. OOOOOOE+OO 

. 795451E-03 

321 

. OOOOOOE+OO 

. 259710E-03 

322 

. OOOOOOE+OO 

693935E-03 

323 

.  OOOOOOE+OO 

483492E-03 

324 

.  OOOOOOE+OO 

. 520607E-03 

325 

.  OOOOOOE+OO 

. 645634E-03 
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326 

. 000000E+00 

297233E-03 

327 

. OOOOOOE+OO 

730852E-03 

328 

.  OOOOOOE+OO 

. 501878E-04 

329 

. OOOOOOE+OO 

. 732867E-03 

330 

. OOOOOOE+OO 

. 192606E-03 

331 

. OOOOOOE+OO 

654736E-03 

332 

.  OOOOOOE+OO 

404790E-03 

333 

. OOOOOOE+OO 

. 508155E-03 

334 

. OOOOOOE+OO 

. 564390E-03 

335 

.  OOOOOOE+OO 

311833E-03 

336 

. OOOOOOE+OO 

656042E-03 

337 

.  OOOOOOE+OO 

. 891797E-04 

338 

. OOOOOOE+OO 

. 672410E-03 

339 

.  OOOOOOE+OO 

.  134443E-03 

340 

. OOOOOOE+OO 

614672E-03 

341 

. OOOOOOE+OO 

334566E-03 

342 

. OOOOOOE+OO 

. 492056E-03 

343 

.  OOOOOOE+OO 

. 490253E-03 

344 

. OOOOOOE+OO 

320515E-03 

345 

. OOOOOOE+OO 

586244E-03 

346 

. OOOOOOE+OO 

. 120718E-03 

347 

. OOOOOOE+OO 

. 614406E-03 

348 

. OOOOOOE+OO 

. 843958E-04 

349 

. OOOOOOE+OO 

574335E-03 

350 

.  OOOOOOE+OO 

272207E-03 

351 

.  OOOOOOE+OO 

. 473097E-03 

352 

.  OOOOOOE+OO 

. 422879E-03 

353 

.  OOOOOOE+OO 

324162E-03 

354 

.  OOOOOOE+OO 

521414E-03 

355 

.  OOOOOOE+OO 

. 145681E-03 

356 

.  OOOOOOE+OO 

. 559102E-03 

357 

.  OOOOOOE+OO 

. 416829E-04 

358 

. 0oncC0E+00 

534224E-03 

359 

. OOOOOOE+OO 

217110E-03 

360 

. OOOOOOE+OO 

. 451961E-03 

361 

.  OOOOOOE+OO 

. 361905E-03 

362 

.  OOOOOOE+OO 

323564E-03 

363 

.  OOOOOOE+OO 

461457E-03 

364 

.  OOOOOOE+OO 

. 164877E-03 

365 

.  OOOOOOE+OO 

. 506671E-03 

366 

.  OOOOOOE+OO 

. 557033E-05 

367 

.  OOOOOOE+OO 

494751E-03 

368 

.  OOOOOOE+OO 

168693E-03 

369 

.  OOOOOOE+OO 

. 429248E-03 

370 

.  OOOOOOE+OO 

. 306958E-03 

371 

.  OOOOOOE+OO 

319432E-03 

372 

.  OOOOOOE+OO 

406239E-03 

373 

.  OOOOOOE+OO 

. 179045E-03 

374 

. OOOOOOE+OO 

.  457227E-03 

375 

.  OOOOOOE+OO 

246287E-04 

376 

.  OOOOOOE+OO 

456254E-03 

377 

. OOOOOOE+OO 

126393E-03 

378 

.  OOOOOOE+OO 

. 405474E-03 

379 

.  OOOOOOE+OO 

. 257656E-03 

380 

. OOOOOOE+OO 

312396E-03 

381 

. OOOOOOE+OO 

355598E-03 
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382 

. 000000E+00 

383 

. OOOOOOE+OO 

384 

. OOOOOOE+OO 

385 

. OOOOOOE+OO 

386 

. OOOOOOE+OO 

387 

. OOOOOOE+OO 

338 

. OOOOOOE+OO 

389 

. OOOOOOE+OO 

390 

. OOOOOOE+OO 

391 

. OOOOOOE+OO 

392 

. OOOOOOE+OO 

393 

. OOOOOOE+OO 

394 

. OOOOOOE+OO 

395 

. OOOOOOE+OO 

396 

. OOOOOOE+OO 

397 

. OOOOOOE+OO 

398 

. OOOOOOE+OO 

399 

. OOOOOOE+OO 

400 

. OOOOOOE+OO 

.  END  OF 

.  188857E-03 
.  410831E-03 
495535E-04 
419005E-03 
896748E-04 
.  381080E-03 
.  213620E-03 
303017E-03 
309343E-03 
.  194925E-03 
. 367499E-03 
697972E-04 
383214E-03 
580256E-04 
. 356444E-03 
. 174472E-03 
291787E-03 
267270E-03 
. 197800E-03 

PROBLEM  #  1  . 


C.  INPUT,  OUTPUT  DATA  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P  « 1 


INPUT  DATA  FOR  PROBLEM  it  1 


PROBLEM  it  1  INPUT  DATA  SOURCEFILE:  diffeq. tst 

THE  NUMBER  OF  INPUT  DELAYS:  L  =  3 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =  3 

THE  VALUE  OF  nstop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  b(L)  ARE: 

.  990000E+00  329900E-02  .  88U46E+00  .  32663 1E+00 


THE  COEFFICIENTS  a(l),  ....  a(N)  ARE: 

. OOOOOOE+OO  •.  871246E+00  323364E+00 

OUTPUT  DATA  FOR  PROBLEM  it  1 
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ns 

x(ns) 

y(ns) 

-3 

. OOOOOOE+OO 

. OOOOOOE+OO 

-2 

. OOOOOOE+OO 

. OOOOOOE+OO 

-1 

. OOOOOOE+OO 

. OOOOOOE+OO 

0 

. 100000E+01 

. 990000E+00 

1 

. OOOCOOE+OQ 

329900E-02 

2 

. OOOOOOE+OO 

. 186127E-01 

3 

. OOOOOOE+OO 

. 937460E-02 

4 

. OOOOOOE+OO 

151494E-01 

S 

. OOOOOOE+OO 

141862E-01 

6 

. OOOOOOE+OO 

. 101675E-01 

7 

. OOOOOOE+OO 

. 172585E-01 

8 

. OOOOOOE+OO 

427105E-02 

9 

. OOOOOOE+OO 

183242E-01 

10 

. OOOOOOE+OO 

185964E-02 

11 

.  OOOOOOE+OO 

. 173460E-01 

12 

. OOOOOOE+OO 

. 754559E-02 

13 

. OOOOOOE+OO 

145113E-01 

14 

. OOOOOOE+OO 

121831E-01 

15 

. OOOOOOE+OO 

. 102029E-01 

16 

. OOOOOOE+OO 

. 153069E-01 

17 

. OOOOOOE+OO 

494966E-02 

18 

. OOOOOOE+OO 

166353E-01 

19 

. OOOOOOE+OO 

637331E-03 

20 

. OOOOOOE+OO 

. 160940E-01 

21 

. OOOOOOE+OO 

. 593455E-02 

22 

. OOOOOOE+OO 

138158E-01 

23 

. OOOOOOE+OO 

103747E-01 

24 

. OOOOOOE+OO 

. 101179E-01 

25 

.  OOOOOOE+OO 

. 135064E-01 

26 

.  OOOOOOE+OO 

546039E-02 

27 

.  OOOOOOE+OO 

150392E-01 

28 

.  OOOOOOE+OO 

. 38985 7E-03 

29 

.  OOOOOOE+OO 

. 148685E-01 

30 

.  OOOOOOE+OO 

.452347E-02 

31 

.  OOOOOOE+OO 

130802E-01 

32 

.  OOOOOOE+OO 

874900E-02 

33 

.  OOOOOOE+OO 

. 993335E-02 

34 

.  OOOOOOE+OO 

. 118522E-01 

35 

. OOOOOOE+OO 

582528E-02 

36 

. OOOOOOE+OO 

135383E-01 

37 

. OOOOOOE+OO 

. 124268E-02 

38 

. OOOOOOE+OO 

. 136788E-0 1 

39 

. OOOOOOE+OO 

. 329511E-02 

40 

. OOOOOOE+OO 

123195E-01 

41 

. OOOOOOE+OO 

729410E-02 

42 

. OOOOOOE+OO 

. 966777E-02 

43 

. OOOOOOE+OO 

.  103386E-01 

44 

. OOOOOOE+OO 

606436E-02 

45 

. OOOOOOE+OO 

121337E-01 

46 

. OOOOOOE+OO 

. 194041E-02 

47 

. OOOOOOE+OO 

. 125324E-01 

48 

. OOOOOOE+OO 

. 223302E-02 

49 

. OOOOOOE+OO 

-.115463E-01 

50 

. OOOOOOE+OO 

599805E-02 
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51 

. OOOOOOE+OO 

. 933758E-02 

52 

. OOOOOOE+OO 

. 895943E-02 

53 

. OOOOOOE+OO 

619577E-02 

54 

. OOOOOOE+OO 

108253E-01 

55 

. OOOOOOE+OO 

. 250088E-02 

56 

. OOOOOOE+OO 

. 114350E-01 

57 

. OOOOOOE+OO 

. 132163E-02 

58 

. OOOOOOE+OO 

107714E-01 

59 

. OOOOOOE+OO 

484913E-02 

60 

. OOOOOOE+OO 

. 895716E-02 

61 

. OOOOOOE+OO 

. 770786E-02 

62 

. OOOOOOE+OO 

623586E-02 

63 

. OOOOOOE+OO 

961187E-02 

64 

. OOOOOOE+OO 

. 294052E-02 

65 

. OOOOOOE+OO 

. 103908E-01 

66 

.  OOOOOOE+OO 

. 546215E-03 

67 

. OOOOOOE+OO 

100038E-01 

68 

. OOOOOOE+OO 

383588E-02 

69 

. OOOOOOE+OO 

. 853911E-02 

70 

. OOOOOOE+OO 

. 657685E-02 

71 

. OOOOOOE+OO 

619928E-02 

72 

. OOOOOOE+OO 

849130E-02 

73 

.  OOOOOOE+OO 

. 327438E-02 

74 

.  OOOOOOE+OO 

. 940263E-02 

75 

. OOOOOOE+OO 

107009E-03 

76 

. OOOOOOE+OO 

-.925082E-02 

77 

. OOOOOOE+OO 

294724E-02 

78 

, OOOOOOE+OO 

. 809434E-02 

79 

. OOOOOOE+OO 

. 555915E-02 

80 

. OOOOOOE+OO 

609913E-02 

81 

.  OOOOOOE+OO 

746081E-02 

82 

.  OOOOOOE+OO 

. 351621E-02 

83 

.  OOOOOOE+OO 

. 847244E-02 

84 

.  OOOOOOE+OO 

650931E-03 

85 

.  OOOOOOE+OO 

851860E-02 

86 

.  OOOOOOE+OO 

217256E-02 

87 

.  OOOOOOE+OO 

. 763228E-02 

86 

.  OOOOOOE+OO 

. 464744E-02 

89 

.  OOOOOOE+OO 

594707E-02 

90 

.  OOOOOOE+OO 

651707E-02 

91 

.  OOOOOOE+OO 

. 367854E-02 

92 

.  OOOOOOE+OO 

. 760104E-02 

93 

.  OOOOOOE+OO 

109753E-02 

94 

.  OOOOOOE+OO 

781188E-02 

95 

. OOOOOOE+OO 

150168E-02 

96 

. OOOOOOE+OO 

. 716097E-02 

97 

. OOOOOOE+OO 

.383442E-02 

98 

. OOOOOOE+OO 

575338E-02 

99 

. OOOOOOE+OO 

565632E-02 

100 

. OOOOOOE+OO 

. 377269E-02 

101 

. OOOOOOE+OO 

. 678848E-02 

102 

. OOOOOOE+OO 

145789E-02 

103 

. OOOOOOE+OO 

-.713439E-02 

104 

. OOOOOOE+OO 

924967E-03 

10S 

. OOOOOOE+OO 

. 668724E-02 

106 

. OOOOOOE+OO 

. 311288E-02 

118 


107 

. 000000E+00 

552713E-02 

108 

. OOOOOOE+OO 

487450E-02 

109 

. OOOOOOE+OO 

. 380890E-02 

110 

. OOOOOOE+OO 

. 603416E-02 

111 

. OOOOOOE+OO 

174225E-02 

112 

. OOOOOOE+OO 

648890E-02 

113 

. OOOOOOE+OO 

-.433302E-03 

114 

. OOOOOOE+OO 

- 621681E-02 

115 

. OOOOOOE+OO 

. 247579E-02 

116 

. OOOOOOE+OO 

527625E-02 

117 

. OOOOOOE+OO 

416731E-02 

118 

. OOOOOOE+OO 

.  379633E-02 

119 

. OOOOOOE+OO 

. 533691E-02 

120 

. OOOOOOE+OO 

195998E-02 

121 

. OOOOOOE+OO 

587736E-02 

122 

. OOOOOOE+OO 

181364E-04 

123 

. OOOOOOE+OO 

. 57544 IE-02 

124 

. OOOOOOE+OO 

. 191633E-02 

125 

. OOOOOOE+OO 

500764E-02 

126 

. OOOOOOE+OO 

353036E-O2 

127 

. OOOOOOE+OO 

. 374322E-02 

128 

. OOOOOOE+OO 

. 469510E-02 

129 

. OOOOOOE+OO 

211967E-02 

130 

. OOOOOOE+OO 

530101E-02 

131 

. OOOOOOE+OO 

. 328528E-03 

132 

. OOOOOOE+OO 

. 530391E-02 

133 

. OOOOOOE+OO 

. 142793E-02 

134 

. OOOOOOE+OO 

472725E-02 

135 

. OOOOOOE+OO 

295917E-02 

136 

. OOOOOOE+OO 

. 365685E-02 

137 

. OOOOOOE+OO 

. 410679E-02 

138 

. OOOOOOE+OO 

222913E-02 

139 

. OOOOOOE+OO 

476052E-Q2 

140 

. OOOOOOE+OO 

. 614133E-03 

141 

.  OOOOOOE+OO 

. 486840E-02 

142 

.  OOOOOOE+OO 

. 100432E-02 

143 

.  OOOOOOE+OO 

444016E-02 

144 

. OOOOOOE+OO 

244927E-02 

145 

.  OOOOOOE+OO 

. 35437  IE-02 

146 

.  OOOOOOE+OO 

. 356971E-02 

147 

,  OOOOOOE+OO 

229544E-02 

148 

. OOOOOOE+OO 

*, 425600E-02 

149 

.  OOOOOOE+OO 

. 845577E-03 

150 

.  OOOOOOE+OO 

. 445029E-02 

151 

.  OOOOOOE+OO 

. 639532E-03 

152 

. OOOOOOE+OO 

415072E-02 

153 

.  OOOOOOE+OO 

199625E-02 

154 

. OOOOOOE+OO 

. 340950E-02 

155 

.  OOOOOOE+OO 

. 308142E-02 

156 

. OOOOOOE+OO 

232500E-02 

157 

.  OOOOOOE+OO 

37  87 19E-02 

158 

.  OOOOOOE+OO 

. 1Q2922E-02 

159 

. OOOOOOE+OO 

. 405139E-02 

160 

.  OOOOOOE+OO 

. 327932E-03 

161 

. OOOOOOE+OO 

386257E-02 

162 

. OOOOOOE+OO 

159578E-02 
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163 

.  OOOOOOE+OO 

. 325921E-02 

164 

. 000000E+00 

. 263934E-02 

165 

. OOOOOOE+OO 

232355E-02 

166 

. 000000E+00 

335342E-02 

167 

. OQOOOOE+OO 

. 117092E-02 

168 

. OOOOOOE+OO 

. 367301E-02 

169 

. OOOOOOE+OO 

. 642165E-04 

170 

. OOOOOOE+OO 

357873E-02 

171 

. OOOOOOE+OO 

124367E-02 

172 

. OOOOOOE+OO 

. 309719E-02 

173 

. OOOOOOE+OO 

. 224077E-02 

174 

. OOOOOOE+OO 

229626E-02 

175 

. OOOOOOE+OO 

295378E-02 

176 

. OOOOOOE+OO 

. 127602E-02 

177 

. OOOOOOE+OO 

. 331600E-02 

178 

. OOOOOOE+OO 

156579E-03 

179 

. OOOOOOE+OO 

330167E-02 

180 

. OOOOOOE+OO 

935856E-03 

181 

. OOOOOOE+OO 

. 292720E-02 

182 

. OOOOOOE+OO 

. 188300E-02 

183 

. OOOOOOE+OO 

224769E-02 

184 

. OOOOOOE+OO 

258711E-02 

185 

. OOOOOOE+OO 

. 134939E-02 

186 

. OOOOOOE+OO 

. 298083E-02 

187 

. OOOOOOE+OO 

339076E-03 

188 

. OOOOOOE+OO 

303338E-02 

189 

. OOOOOOE+OO 

668474E-03 

190 

. OOOOOOE+OO 

. 275246E-02 

191 

, OOOOOOE+OO 

. 156329E-02 

192 

. OOOOOOE+OO 

21819 IE -02 

193 

. OOOOOOE+OO 

225206E-02 

194 

. OOOOOOE+OO 

. 139547E-02 

195 

. OOOOOOE+OO 

, 266765E-02 

196 

. OOOOOOE+OO 

-.487564E-03 

197 

.  OOOOOOE+OO 

277542E-02 

198 

. OOOOOOE+OO 

437834E-03 

199 

.  OOOOOOE+OO 

. 257574E-02 

200 

. OOOOOOE+OO 

. 127893E-02 

201 

. OOOOOOE+OO 

210252E-02 

202 

. OOOOOOE+OO 

1947 17E-02 

203 

. OOOOOOE+OO 

. 141825E-02 

204 

. OOOOOOE+OO 

. 237634E-02 

205 

.  OOOOOOE+OO 

606003E-03 

206 

. OOOOOOE+OO 

252899E-02 

207 

. OOOOOOE+OO 

240445E-03 

208 

. OOOOOOE+OO 

. 239933E-02 

209 

. OOOOOOE+OO 

. 102727E-02 

210 

. OOOOOOE+OO 

201266E-02 

211 

. OOOOOOE+OO 

167086E-02 

212 

. OOOOOOE+OO 

. 142134E-02 

213 

. OOOOOOE+OO 

. 21065SE-02 

214 

. OOOOOOE+OO 

698037E-03 

215 

. OOOOOOE+OO 

229493E-02 

216 

. OOOOOOE+OO 

730215E-04 

217 

. OOOOOOE+OO 

. 222517E-02 

218 

. OOOOOOE+OO 

. 805719E-03 
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219 

. OOOOOOE+OO 

191506E-02 

220 

. OOOOOOE+OO 

142152E-02 

221 

. OOOOOOE+OO 

. 140795E-02 

222 

. OOOOOOE+OO 

. 185775E-02 

223 

. OOOOOOE+OO 

-.767001E-03 

224 

. OOOOOOE+OO 

207384E-02 

225 

. OOOOOOE+OO 

. 675152E-04 

226 

. OOOOOOE+OO 

. 205485E-02 

227 

. OOOOOOE+OO 

. 611783E-03 

228 

. OOOOOOE+OO 

181211E-02 

229 

. OOOOOOE+OO 

119748E-02 

230 

. OOOOOOE+OO 

. 138096E-02 

231 

. OOOOOOE+OO 

. 162927E-02 

232 

. OOrtOOOE+OO 

815938E-03 

233 

. OOOOOOE+OO 

186605E-02 

234 

. OOOOOOE+OO 

. 184036E-03 

235 

. OOOOOOE+OO 

. 188963E-02 

236 

. OOOOOOE+OO 

. 443071E-03 

237 

. OOOOOOE+OO 

170584E-02 

238 

. OOOOOOE+OO 

997063E-03 

239 

. OOOOOOE+OO 

. 134294E-02 

240 

. OOOOOOE+OO 

. 142030E-02 

241 

. OOOOOOE+OO 

847613E-03 

242 

. OOOOOOE+OO 

167168E-02 

243 

. OOOOOOE+OO 

. 279207E-03 

244 

. OOOOOOE+OO 

. 173054E-02 

245 

. OOOOOOE+OO 

. 297304E-03 

246 

. OOOOOOE+OO 

159801E-02 

247 

. OOOOOOE+OO 

818618E-03 

248 

. OOOOOOE+OO 

. 129612E-02 

249 

. OOOOOOE+OO 

. 122996E-02 

250 

. OOOOOOE+OO 

864528E-03 

251 

. OOOOOOE+OO 

149071E-02 

252 

. OOOOOOE+OO 

. 355493E-03 

253 

. OOOOOOE+OO 

. 157833E-02 

254 

.  OOOOOOE+OO 

. 172321E-03 

255 

,  OOOOOOE+OO 

149007E-02 

256 

.  OOOOOOE+OO 

660510E-03 

257 

.  OOOOOOE+OO 

. 124250E-02 

258 

.  OOOOOOE+OO 

. 105730E-02 

259 

.  OOOOOOE+OO 

868935E-03 

260 

.  OOOOOOE+OO 

«. 132295E-02 

261 

.  OOOOOOE+OO 

. 415163E-03 

262 

,  OOOOOOE+OO 

. 14336CE-02 

263 

.  OOOOOOE+OO 

. 660854E-04 

264 

.  OOOOOOE+OO 

138326E-02 

265 

.  OOOOOOE+OO 

521150E-03 

266 

.  OOOOOOE+OO 

. 118379E-02 

267 

.  OOOOOOE+OO 

. 9Q1347E-03 

268 

.  OOOOOOE+OO 

862854E-03 

269 

.  OOOOOOE+OO 

116809E-02 

270 

.  OOOOOOE+OO 

. 460295E-03 

271 

.  OOOOOOE+OO 

. 129671E-02 

272 

.  OOOOOOE+OO 

233112E-04 

273 

.  OOOOOOE+OO 

127860E-02 

274 

.  OOOOOOE+OO 

399Q00E-03 
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275 

. 000000E+00 

. 11215 IE-02 

276 

. OOOOOOE+OO 

.  761079E-03 

277 

. 000000E+00 

848089E-03 

278 

. OOOOOOE+OO 

102574E-02 

279 

. OOOOOOE+OO 

. 492789E-03 

280 

. OOOOOOE+OO 

. 116792E-02 

281 

. OOOOOOE+OO 

-.976517E-04 

282 

. OOOOOOE+OO 

117689E-02 

283 

. OOOOOOE+OO 

292583E-03 

284 

. OOOOOOE+OO 

. 105694E-02 

285 

. OOOOOOE+OO 

. 635477E-03 

286 

. OOOOOOE+OO 

826244E-03 

287 

. OOOOOOE+OO 

89S433E-03 

288 

. OOOOOOE+OO 

. 514371E-03 

289 

. OOOOOOE+OO 

. 104732E-02 

290 

.  OOOOOOE+OO 

158593E-03 

291 

.  OOOOOOE+OO 

107880E-02 

292 

.  OOOOOOE+OO 

200492E-03 

293 

.  OOOOOOE+OO 

. 991185E-03 

294 

.  OOOOOOE+OO 

. 523523E-03 

295 

.  OOOOOOE+OO 

798734E-03 

296 

.  OOOOOCE+OO 

776631E-03 

297 

.  OOOOOOE+OO 

. 526606E-03 

298 

.  OOOOOOE+OO 

. 934919E-03 

299 

.  00O0Q0E+00 

207668E-03 

300 

.  OOOOOOE+OO 

984830E-03 

301 

.  OOOOOOE+OO 

121389E-03 

302 

.  OOOOOOE+OO 

. 925181E-03 

303 

.  OOOOOOE+OO 

. 424218E-03 

304 

.  OOOOOOE+OO 

766808E-03 

305 

.  OOOOOOE+OO 

668769E-03 

306 

.  OOOOOOE+OO 

. 53Q901E-03 

307 

.  OOOOOOE+OO 

. 830620E-03 

308 

.  OOOOOOE+OO 

246290E-03 

309 

.  OOOOOOE+OO 

895349E-03 

310 

.  OOOOOOE+OO 

-.540135E-04 

311 

.  OOOOOOE+OO 

. 859710E-03 

312 

.  OOOOOOE+OO 

. 336583E-03 

313 

.  OOOOOOE+OO 

731553E-03 

314 

. OOOOOOE+OO 

571245E-03 

315 

.  OOOOOOE+OO 

. 528524E-03 

316 

. OOOOOOE+OO 

.734253E-03 

317 

.  OOOOOOE+OO 

275754E-03 

318 

. OOOOOOE+OO 

810621E-03 

319 

. OOOOOOE+OO 

. 281864E-05 

320 

. OOOOOOE+OO 

. 795419E-03 

321 

. OOOOOOE+OO 

.259670E-03 

322 

. OOOOOOE+OO 

693917E-03 

323 

. OOOOOOE+OO 

483446E-03 

324 

. OOOOOOE+OO 

. 520605E-03 

325 

. OOOOOOE+OO 

.645588E-03 

326 

. OOOOOOE+OO 

297246E-03 

327 

. OOOOOOE+OO 

7308UE-03 

328 

. OOOOOOE+OO 

. 502140E-04 

329 

. OOOOOOE+OO 

. 732835E-03 

330 

. OOOOOOE+OO 

. 192569E-03 
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331 

. 000000E+00 

654717E-03 

332 

. OOOOOOE+OO 

404747E-03 

333 

. OOOOOOE+OO 

. 508149E-03 

334 

. OOOOOOE+OO 

. 564346E-03 

335 

. OOOOOOE+OO 

311842E-03 

336 

.  OOOOOOE+OO 

-.656002E-03 

337 

.  OOOOOOE+OO 

. 892021E-04 

338 

. OOOOOOE+OO 

. 672377E-03 

339 

. OOOOOOE+OO 

.  134410E-03 

340 

. OOOOOOE+OO 

614651E-03 

341 

. OOOOOOE+OO 

334527E-03 

342 

. OOOOOOE+OO 

.  492049E-03 

343 

. OOOOOOE+OO 

.49  02  HE -03 

344 

. OOOOOOE+OO 

320521E-03 

345 

. OOOOOOE+OO 

586206E-03 

346 

. OOOOOOE+OO 

. 120736E-03 

347 

. OOOOOOE+OO 

. 614374E-03 

348 

. OOOOOOE+OO 

. 843669E-04 

349 

. OOOOOOE+OO 

574313E-03 

350 

. OOOOOOE+OO 

272171E-03 

351 

. OOOOOOE+OO 

.473087E-03 

352 

. OOOOOOE+OO 

. 422840E-03 

353 

. OOOOOOE+OO 

324164E-03 

354 

. OOOOOOE+OO 

521377E-03 

355 

. OOOOOOE+OO 

. 145696E-03 

356 

. OOOOOOE+OO 

.559070E-03 

357 

. OOOOOOE+OO 

. 416575E-04 

358 

. OOOOOOE+OO 

534201E-03 

359 

.  OOOCOOE+OO 

217077E-03 

360 

.  OOOOOOE+OO 

. 451950E-03 

361 

.  OOOOOOE+OO 

. 361869E-03 

362 

.  OOOOOOE+OO 

«. 323564E-03 

363 

.  OOOOOOE+OO 

461421E-03 

364 

.  OOOOOOE+OO 

. 164889E-03 

365 

.  OOOOOOE+OO 

. 506640E-03 

366 

.  OOCCOOE+OO 

. 554827E-Q5 

367 

.  OOOOOOE+OO 

494727E-03 

368 

. OOOOOOE+OO 

168663E-03 

369 

. OOOOOOE+OO 

. 429235E-03 

370 

. OOOOOOE+OO 

.306924E-03 

371 

. OOOOOOE+OO 

319430E-03 

372 

. OOOOOOE+OO 

-.406206E-03 

373 

.  OOOOOOE+OO 

. 179054E-03 

374 

. OOOOOOE+OO 

. 457197E-03 

375 

. OOOOOOE+OO 

246476E-04 

376 

. OOOOOOE+OO 

456231E-03 

377 

. OOOOOOE+OO 

126367E-03 

378 

. OOOOOOE+OO 

. 405459E-03 

379 

. OOOOOOE+OO 

. 257625E-03 

380 

. OOOOOOE+OO 

312392E-03 

381 

. OOOOOOE+OO 

355566E-03 

382 

. OOOOOOE+OO 

. 188864E-03 

383 

. OOOOOOE+OO 

. 410802E-03 

384 

. OOOOOOE+OO 

495696E-04 

385 

. OOOOOOE+OO 

-.418981E-03 

386 

. OOOOOOE+OO 

896512E-04 
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387 

. OOOOOOE+OO 

.  381065E-03 

388 

. OOOOOOE+OO 

. 213592E-03 

389 

. OOOOOOE+OO 

3030UE-03 

390 

. OOOOOOE+OO 

309313E-03 

391 

. OOOOOOE+OO 

. 194929E-03 

392 

. OOOOOOE+OO 

. 367471E-03 

393 

.  OOOOOOE+OO 

698106E-04 

394 

.  OOOOOOE+OO 

383191E-03 

395 

.  OOOOOOE+OO 

580047E-04 

396 

.  OOOOOOE+OO 

. 356428E-03 

397 

.  OOOOOOE+OO 

. 174446E-03 

398 

.  OOOOOOE+OO 

291779E-03 

399 

.  OOOOOOE+OO 

267242E-03 

400 

.  OOOOOOE+OO 

. 197802E-03 

END  OF  PROBLEM  #  1 


D.  INPUT,  OUTPUT  DATA  OF  THE  AUGMENTED  HR  NOTCH  FILTER:  P 


INPUT  DATA  FOR  PROBLEM  #  1 


PROBLEM  #  1  INPUT  DATA  SOURCEFILE:  dif feq.  tst 

THE  NUMBER  OF  INPUT  DELAYS:  L  «*  4 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  *»  4 

THE  VALUE  OF  nstop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  b(L)  ARE: 

. 990000E+00  329900E-02  . I86124E-01  .614081E+00  -.862533E+00 


THE  COEFFICIENTS  a(l),  a(N)  ARE: 

.  OOOOOOE+OO  .OOOOOOE+OO  -.610814E+00  . 853908E+00 


OUTPUT  DATA  FOR  PROBLEM  #  1 


ns  x(ns) 

-4  .  OOOOOOE+OO 
-3  . OOOOOOE+OO 
-2  . OOOOOOE+OO 


y(ns) 

.  OOOOOOE+OO 
.  OOOOOOE+OO 
.  OOOOOOE+OO 
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-1 

. OOOOOOE+OO 

. OOOOOOE+OO 

0 

. 100000E+01 

. 990000E+00 

1 

. OOOOOOE+OO 

329900E-02 

2 

. OOOOOOE+OO 

. 186124E-01 

3 

. OOOOOOE+OO 

. 937506E-02 

4 

. OOOOOOE+OO 

151490E-01 

5 

. OOOOOOE+OO 

141858E-01 

6 

. OOOOOOE+OO 

. 101669E-01 

7 

. OOOOOOE+OO 

. 172586E-01 

8 

. OOOOOOE+OO 

427097E-02 

9 

. OOOOOOE+OO 

183234E-01 

10 

. OOOOOOE+OO 

186026E-02 

11 

. OOOOOOE+OO 

. 173461E-01 

12 

. OOOOOOE+OO 

. 754517E-02 

13 

. OOOOOOE+OO 

145102E-01 

14 

.  OOOOOOE+OO 

121837E-01 

15 

. OOOOOOE+OO 

. 102032E-01 

16 

. OOOOOOE+OO 

. 153059E-01 

17 

. OOOOOOE+OO 

494841E-02 

18 

. OOOOOOE+OO 

166360E-01 

19 

. OOOOOOE+OO 

636442E-03 

20 

. OOOOOOE+OO 

. 160924E-01 

21 

. OOOOOOE+OO 

. 593604E-02 

22 

. OOOOOOE+OO 

138169E-01 

23 

, OOOOOOE+OO 

103729E-01 

24 

. OOOOOOE+OO 

. 101 156E-01 

25 

. OOOOOOE+OO 

. 135084E-01 

26 

. OOOOOOE+OO 

546244E-02 

27 

.  OOOOOOE+OO 

150363E-01 

28 

.  OOOOOOE+OO 

. 386692E-03 

29 

.  OOOOOOE+OO 

. 148715E-01 

30 

.  OOOOOOE+OO 

. 451996E-02 

31 

.  OOOOOOE+OO 

130758E-01 

32 

.  OOOOOOE+OO 

875350E*02 

33 

.  OOOOOOE+OO 

. 993800E-02 

34 

.  OOOOOOE+OO 

. 118465E-01 

35 

. OOOOOOE+OO 

581878E-02 

36 

.  OOOOOOE+OO 

135450E-01 

37 

.  OOOOOOE+OO 

. 12S013E-02 

38 

.  OOOOOOE+OO 

. 136700E-01 

39 

. OOOOOOE+OO 

.330475E-02 

40 

.  OOOOOOE+OO 

123297E-01 

41 

.  OOOOOOE+OO 

728235E-02 

42 

.  OOOOOOE+OO 

. 965436E-02 

43 

. OOOOOOE+OO 

•103S31E-01 

44 

.  OOOOOOE+OO 

608030E-02 

45 

.  OOOOOOE+OO 

1211S3E-01 

46 

. OOOOOOE+OO 

.  192010E-02 

47 

. OuOOOOE+OO 

.  125S45E-01 

48 

. OOOOOOE+OO 

. 220828E-02 

49 

.  OOOOOOE+OO 

115183E-01 

50 

. OOOOOOE+OO 

602890E-02 

51 

. OOOOOOE+OO 

. 937158E-02 

52 

. OOQOOCE+QO 

. 892123E-02 

53 

. OOOOOOE+OO 

61S305E-02 

54 

. OOOOOOE+OO 

108724E-01 
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55 

. OOOOOOE+OO 

. 255325E-02 

56 

.  OOOOOOE+OO 

. 113763E-01 

57 

. OOOOOOE+OO 

. 138688E-02 

58 

. OOOOOOE+OO 

108436E-01 

59 

. OOOOOOE+OO 

476854E-02 

60 

. OOOOOOE+OO 

. 886717E-02 

61 

. OOOOOOE+OO 

. 780770E-02 

62 

. OOOOOOE+OO 

634675E-02 

63 

. OOOOOOE+OO 

948809E-02 

64 

. OOOOOOE+OO 

. 280269E-02 

65 

. OOOOOOE+OO 

. 105437E-01 

66 

.  OOOOOOE+OO 

. 375915E-03 

67 

.  OOOOOOE+OO 

981388E-02 

68 

. OOOOOOE+OO 

404702E-02 

69 

. OOQQOQE+QO 

. 877377E-02 

70 

. OOOOOOE+OC 

. 631545E-02 

71 

. OOOOOOE+OO 

590817E-02 

72 

. OOOOOOE+OO 

-.881493E-02 

73 

. OOOOOOE+OO 

. 363443E-02 

74 

. OOOOOOE+OO 

. 90Q161E-02 

75 

. OOOOOOE+OO 

. 339247E-03 

76 

. QOOOOOE+OO 

974710E-02 

77 

. OOOOOOE+OO 

239464E-02 

78 

. OOOOOOE+OO 

. 747933E-02 

79 

.  OOOOOOE+OO 

, 624335E-02 

80 

,  OOOOOOE+OO 

•. 6860328*02 

81 

.  OOOOOOE+OO 

661345E-02 

82 

.  OOOOOOE+OO 

. 257313E-02 

83 

.  OOOOOOE+OO 

. 952163E-02 

84 

.  OOOOOOE+OO 

•. 181849E-02 

85 

.  OOOOOOE+OO 

721896E-02 

66 

.  OOOOOOE+OO 

361872E-02 

87 

.  OOOOOOE+OO 

.924135E-02 

88 

.  OQOOQOE+OO 

. 285663E“02 

89 

.  OOOOOOE+OO 

39S398E-02 

90 

.  OOOOOOE+OO 

873480E-02 

91 

.  OOOOOOE+OO 

. 614640E-02 

92 

.  OOOOOOE+OO 

. 46S444E*02 

93 

.  OOOOOOE+OO 

. 195900E-02 

94 

.  OOOOOOE+OO 

112130E*0i 

95 

.  OOOOOOE+OO 

. 228329E-02 

96 

.  OOOOOOE+OO 

. 294S66E-02 

97 

.  OOOOOOE+OO 

.852188E-02 

98 

. OOOOOOE+OO 

109696E-01 

99 

. OOOOOOE+OO 

. 148640E-03 

100 

. OOOOOOE+OO 

26B740E-02 

101 

. OOOOOOE+OO 

. 139773E-01 

102 

. OOOOOOE+OO 

94S779E-02 

103 

. OOOOOOE+OO 

. 176843E-02 

104 

. OOOOOOE+OO 

108323E-01 

10S 

. OOOOOOE+OO 

. 177122E-01 

106 

. OOOOOOE+OO 

*. 915626E-02 

107 

. OOOOOOE+OO 

. 612659E-02 

108 

. OOOOOOE+OO 

200667E-01 

109 

. OOOOOOE+OO 

. 207174E-01 

110 

. OOOOOOE+OO 

12’824E-01 

126 


Ill 

. OOOOOOE+OO 

. 191976E-01 

112 

. OOOOOOE+OO 

-. 297913E-01 

m 

. OOOOOOE+OO 

. 254984E-01 

114 

. OOOOOOE+OG 

-. 226412E-01 

115 

. OOOOOOE+OO 

. 345899E-01 

i  16 

.OOOOOOE+OO 

410138E-01 

117 

. OOOOOOE+OO 

. 356029E-01 

118 

. OOOOOOE+OO 

-. 404615E-01 

119 

. OOOOOOE+OO 

. 545884E-01 

120 

. OOOOOOE+OO 

-. 567688E-01 

121 

. OOOOOOE+OO 

. 551160E-01 

122 

. OOOOOOE+OO 

678938E-01 

123 

. OOOOOOE+OO 

. 812886E-01 

124 

. OOOOOOE+OO 

821409E-C1 

125 

. OOOOOOE+OO 

. 885345E-01 

126 

. OOOOOOE+OO 

-. 107627E+00 

127 

. OOOOOOE+OO 

. 119586E+00 

128 

. OOOOOOE+OO 

-. 124219E+00 

129 

. OOOOOOE+OO 

. 14134 1E+00 

130 

. OOOOOOE+OO 

-. 164948E+00 

131 

. OOOOOOE+OO 

. 177990E+00 

132 

. OOOOOOE+OO 

-. 192404E+00 

133 

. OOOOOOE+OO 

. 221445E+00 

134 

. OOOOOOE+OO 

-. 249570E+00 

135 

. OOOOOOE+OO 

. 269510E+00 

136 

. OOOOOOE+OO 

-. 299557E+00 

137 

.gooooof foo 

. 341534E+00 

138 

. OOOOOOE+OO 

-. 377730E+00 

139 

. OOOOOOE+OO 

. 413111E+00 

140 

. OOOOOOE+OO 

-. 464408E+00 

141 

. OOOOOOE+OO 

522361E+00 

142 

. OOOOOOE+OO 

-. 574880E+00 

143 

. OOOOOOE+OO 

. 636425E+00 

144 

. OOOOOOE+OO 

-. 715627E400 

145 

. OOOOOOE+OO 

. 797194E+00 

146 

. OOOOOOE+OO 

-. 879632E+00 

147 

. OOOOOOE+OO 

. 980564E+00 

148 

. OOOOOOE+OO 

-.  109802E+01 

149 

. OOOOOOE+OO 

.  121802E+01 

150 

.OOOOOOE+OO 

135007E+01 

151 

. OOOOOOE+OO 

. 150800E+01 

152 

. OOOOOOE+OO 

-.  168159E+01 

153 

. OOOOOOE+OO 

.  186472E+01 

154 

. OOOOOOE+OO 

-.  207394E+01 

155 

. OOOOOOE+OO 

. 231483E+01 

156 

. OOOOOOE+OO 

-. 257492E+01 

157 

. OOOOOOE+OO 

. 285909E+01 

158 

. OOOOOOE+OO 

318488E+01 

159 

. OOOOOOE+OO 

. 354945E+01 

160 

. OOOOOOE+OO 

-.  394512E+01 

161 

. OOOOOOE+OO 

. 438677E+01 

162 

.  OOOOOOE+OO 

-.488765 1+01 

163 

. OOOOOOE+OO 

. 544063E+01 

164 

. OOOOOOE+OO 

-. 604826E+01 

165 

. OOOOOOE+OO 

. 673134E+01 

166 

. OOOOOOE+OO 

-. 74968 1E+01 
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167 

. OOOOOOE+OO 

.  834016E+01 

168 

. OOOOOCE+OO 

927626E+01 

16)i 

. OOOOOOE+OO 

. 103271E+02 

170 

. OOOOOOE+OO 

1 14S59E+02 

171 

. OOOOOOE+OO 

. 127878E+02 

172 

. OOOOOCE+OO 

142290E+02 

173 

. OOOOOOE+OO 

. 158402E+02 

174 

. OOOOOOE+OO 

176274E+02 

175 

. OOOOOOE+OO 

. 1961G9E+02 

176 

. OOOOOOE+OO 

218257E+02 

177 

. OOOOOOE+OO 

. 242932E+02 

178 

. OOOOOOE+OO 

270308E+02 

179 

. OOOOOOE+OO 

. 300773E+02 

180 

. OOOOOOE+OO 

334757E+02 

181 

. OOOOOOE+OO 

. 372549E+02 

182 

. OOOOOOE+OO 

414534E+02 

183 

. OOOOOOE+OO 

.  461307E+02 

184 

.  OOOOOOE+OO 

513410E+C2 

185 

.  OOOOOOE+OO 

. 571326E+02 

186 

.  OOOOOOE+OO 

635747E+02 

187 

.  OOOOOOE+OO 

. 707512E+02 

188 

.  OOOOOOE+OO 

787379E+02 

189 

.  OOOOOOE+OO 

. 876183E+02 

190 

.  OOOOOOE+OO 

975028E+02 

191 

.  OOOOOOE+OO 

. 108509E+03 

192 

.  OOOOOOE+OO 

120753E+03 

193 

.  OOOOOOE+OO 

. 134374E+03 

194 

.  OOOOOOE+OO 

149537E+03 

195 

.  OOOOOOE+OO 

. 166415E+03 

196 

.  OOOOOOE+OO 

185190E+03 

197 

.  OOOOOOE+UO 

. 206083E+03 

198 

.  OOOOOOE+OO 

229340E+03 

199 

.  OOOOOOE+OO 

. 255219E+03 

200 

.  OOOOOOE+OO 

284013E+03 

201 

.  OOOOOOE+OO 

.  316059E+0?- 

202 

.  OOOOOOE+OO 

351727E+03 

203 

. OOOOOOE+OO 

. 391413E+03 

204 

. OOOOOOE+OO 

435575E+03 

205 

. OOOOOOE+OO 

. 484725E+03 

206 

.  OOOOOOE+OO 

539423E+03 

207 

. OOOOOOE+OO 

. 600286E+03 

208 

.  OOOOOOE+OO 

668018E+03 

209 

. OOOOOOE+OO 

. 743396E+03 

210 

. OOOOOOE+OO 

». 827281E+03 

211 

. OOOOOOE+OO 

. 920624E+03 

212 

. OOOOOOE+OO 

102450E+04 

213 

.  OOOOOOE+OO 

. 114011E+04 

214 

. OOOOOOE+OO 

126875E+04 

215 

. OOOOOOE+OO 

.  141191E+04 

216 

. OOOOOOE+OO 

157123E+04 

217 

. OOOOOOE+OO 

. 174852B+04 

218 

. OOOOOOE+OO 

194581E+04 

219 

. OOOOOOE+OO 

. 216537E+04 

220 

. OOOOOOE+OO 

240970E+04 

221 

. OOOOOOE+OO 

. 268160E+04 

222 

. OOOOOOE+OO 

298418E+04 
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223 

. 000000E+00 

. 332090E+04 

224 

. 000000E+00 

369562E+04 

225 

. 000000E+00 

. 411262E+04 

226 

. 000000E+00 

457667E+04 

227 

. OOOOOOE+OO 

. 509308E+04 

228 

. 000000E+00 

566777E+04 

229 

. OOOOOOE+OO 

. 630729E+04 

230 

. OOOOOOE+OO 

701898E+04 

231 

. OOOOOOE+OO 

. 761098E+04 

232 

. OOOOOOE+OO 

86923?E+04 

233 

. OOOOOOE+OO 

. 967314E+04 

234 

. OOOOOOE+OO 

107646E+05 

235 

. OOOOOOE+OO 

. 119793E+05 

236 

. OOOOOCE+OO 

133309E+05 

237 

. OOOOOOE+OO 

. 14835 1E+05 

238 

. OOOOOOE+OO 

16509 1E+05 

239 

. OOOOOOE+OO 

. 183719E+05 

240 

. OOOOOOE+OO 

204449E+05 

241 

. OOOOOOE+OO 

. 227518E+05 

242 

. OOOOOOE+OO 

25319 1E+05 

243 

. OOOOOOE+OO 

. 281760E+05 

244 

. OOOOOOE+OO 

313552E+05 

245 

. OOOOOOE+OO 

. 348932E+05 

246 

. OOOOOOE+OO 

388304E+05 

247 

. OOOOOOE+OO 

. 432119E+05 

248 

. OOOOOOE+OO 

48087 7E+05 

249 

. OOOOOOE+OO 

. 535138E+05 

250 

. OOOOOOE+OO 

595520E+05 

251 

. OOOOOOE+OO 

. 662716E+05 

232 

. OOOOOOE+OO 

737494E+05 

253 

. OOOOOOE+OO 

. 820710E+05 

254 

. OOOOOOE+OO 

913316E+05 

255 

. OOOOOOE+OO 

. 101637E+06 

256 

. OOOOOOE+OO 

113105E+06 

257 

.  OOOOOOE+OO 

. 125868E+06 

258 

.  OOOOOOE+OO 

140070E+06 

259 

.  OOOOOOE+OO 

. 155875E+06 

260 

. OOOOOOE+OO 

173463E+06 

261 

. OOOOOOE+OO 

. 193036E+06 

262 

. OOOOOOE+OO 

214818E+06 

263 

.  OOOOOOE+OO 

. 239057E+06 

264 

.  OOOOOOE+OO 

266031E+06 

265 

.  OOOOOOE+OO 

. 296049E+06 

266 

. OOOOOOE+OO 

329454E+06 

267 

.  OOOOOOE+OO 

. 366628E+C6 

268 

. OOOOOOE+OO 

407997E+06 

269 

. OOOOOOE+OO 

. 454033E+06 

270 

. OOOOOOE+OO 

505265E+06 

271 

. OOOOOOE+OO 

. 562277E+06 

272 

. OOOOOOE+OO 

625722E+06 

273 

. OOOOOOE+OO 

. 696325E+06 

274 

. OOOOOOE+OO 

774896E+06 

275 

. OOOOOOE+OO 

. 862332E+06 

276 

. OOOOOOE+OO 

959634E+06 

277 

. OOOOOOE+OO 

. 106792E+07 

278 

. OOOOOOE+OO 

118841E+07 
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279 

.  OOOOOOE+OO 

. 13225 1E+07 

280 

. OOOOOOE+OO 

147174E+07 

281 

. OOOOOOE+OO 

. 163780E+07 

282 

. OOOOOOE+OO 

182260E+07 

283 

. OOOOOOE+OO 

. 202826E+07 

284 

.  OOOOOOE+OO 

225712E+07 

285 

. OOOOOOE+OO 

. 251180E+07 

286 

. OOOOOOE+OO 

279523E+07 

287 

.  OOOOOOE+OO 

. 311063E+07 

288 

.  OOOOOOE+OO 

346162E+07 

289 

.  OOOOOOE+OO 

. 385221E+07 

290 

.  OOOOOOE+OO 

428688E+07 

291 

.  OOOOOOE+OO 

. 477059E+07 

292 

.  OOOOOOE+OO 

530889E+07 

293 

.  OOOOOOE+OO 

. 590792E+07 

294 

.  OOOOOOE+OO 

657455E+07 

295 

.  OOOOOOE+OO 

. /31639E+07 

296 

.  OOOOOOE+OO 

814194E+07 

297 

.  OOOOOOE+OO 

, 906065E+07 

298 

.  OOOOOOE+OO 

100830E+08 

299 

.  OOOOOOE+OO 

.  112207E+08 

300 

.  OOOOOOE+OO 

124868E+08 

301 

.  OOOOOOE+OO 

. 138958E+08 

302 

.  OOOOOOE+OO 

154637E+08 

303 

,  OOOOOOE+OO 

. 172086E+08 

304 

.  OOOOOOE+OO 

1915Q4E+08 

305 

, OOOOOOE+OO 

. 213112E+08 

306 

.  OOOOOOE+OO 

237159E+08 

307 

.  OOOOOOE+OO 

. 263919E+08 

308 

.  OOOOOOE+OO 

293698E+08 

309 

.  OOOOOOE+OO 

. 326838E+08 

310 

.  OOOOOOE+OO 

363717E+08 

311 

. OOOOOOE+OO 

. 404757E+08 

312 

. OOOOOOE+OO 

450429E+08 

313 

. OOOOOOE+OO 

. 501253E+08 

314 

. OOOOOOE+OO 

557812E+08 

315 

. OOOOOOE+OO 

. 620754E+08 

316 

. OOOOOOE+OO 

690797E+08 

317 

. OOOOOOE+OO 

. 768744E+08 

318 

. OOOOOOE+OO 

855485E+08 

319 

. OOOOOOE+OO 

. 952015E+08 

320 

. OOOOOOE+OO 

105944L+09 

321 

. OOOOOOE+OO 

. 117898E *09 

322 

. OOOOOOE+OO 

«. 131201E-09 

323 

.  OOOOOOE+OO 

. 146005E+09 

324 

.  OOOOOOE+OO 

162480E+09 

325 

.  OOOOOOE+OO 

. 18Q813E+09 

326 

.  OOOOOOE+OO 

201216E+09 

327 

.  OOOOOOE+OO 

. 223920E+09 

328 

.  OOOOOOE+OO 

249186E+09 

329 

.  OOOOOOE+OO 

. 277303E+09 

330 

.  OOOOOOE+OO 

308593E+09 

331 

.  OOOOOOE+OO 

v 343413E+09 

332 

.  OOOOOOE+OO 

382163E+09 

333 

.  OOOOOOE+OO 

. 425284E+G9 

334 

. OOOOOOE+OO 

473272E+09 

130 


335 

. OOOOOOE+OO 

. 526674E+09 

336 

. OOOOOOE+OO 

586101E+09 

337 

. OOOOOOE+OO 

. 652235E+09 

338 

. OOOOOOE+OO 

725830E+09 

339 

. OOOOOOE+OO 

.  807730E+09 

340 

. OOOOOOE+OO 

89887 1E+09 

341 

. OOOOOOE+OO 

. 100030E+10 

342 

. OOOOOOE+OO 

111316E+10 

343 

. OOOOOOE+OO 

. 123877E+10 

344 

.  OOOOOOE+OO 

137855E+10 

345 

.  OOOOOOE+OO 

. 153410E+10 

346 

. OOOOOOE+OO 

170720E+10 

347 

. OOOOOOE+OO 

. 189983E+10 

348 

. OOOOOOE+OO 

211420E+10 

349 

. OOOOOOE+OO 

. 235276E+10 

350 

. OOOOOOE+OO 

261823E+10 

351 

. OOOOOOE+OO 

. 291366E+10 

352 

.OOOOOOE+OO  ' 

324243E+10 

353 

. OOOOOOE+OO 

. 360829E+10 

354 

. OOOOOOE+OO 

401544E+10 

355 

. OOOOOOE+OO 

. 446852E+10 

356 

. OOOOOOE+OO 

497273E+10 

357 

. OOOOOOE+OO 

. 553384E+10 

358 

. OOOOOOE+OO 

615825E+10 

359 

. OOOOOOE+OO 

. 685312E+10 

360 

. OOOOOOE+OO 

762640E+10 

361 

. OOOOOOE+OO 

. 848693E+10 

362 

. OOOOOOE+OO 

944456E+10 

363 

. OOOOOOE+OO 

. 105102E+11 

364 

. OOOOOOE+OO 

116962E+11 

365 

. OOOOOOE+OO 

, 130159E+11 

366 

. OOOOOOE+OO 

144846E+11 

367 

. OOOOOOE+OO 

. 161190E+11 

368 

.  OOOOOOE+OO 

179378E+11 

369 

.  OOOOOOE+OO 

. 199618E+11 

370 

.  OOOOOOE+OO 

222142E+11 

371 

.  OOOOOOE+OO 

. 2472Q8E+11 

372 

.  OOOOOOE+OO 

275102E+U 

373 

.  OCOOdOE+OO 

. 306143E+11 

374 

.  OOOOOOE+OO 

340687E+11 

375 

.  OOOOOOE+OO 

. 379129E+11 

376 

.  OOOOOOE+OO 

*. 421908E+11 

377 

.  OOOOOOE+OO 

. 469514E+11 

378 

. OOOOOOE+OO 

522492E+11 

379 

. OOOOOOE+OO 

. 581448E+11 

380 

.  OOOOOOE+OO 

647056E+11 

381 

.  OOOOOOE+OO 

. 720067E+11 

382 

.  OOOOOOE+OO 

801317E+11 

383 

.  OOOOOOE+OO 

. 891734E+11 

384 

. OOOOOOE+OO 

992354E+11 

385 

. OOOOOOE+OO 

. 110433E+12 

386 

. OOOOOOE+OO 

122893E+12 

387 

. OOOOOOE+OO 

. 136760E+12 

388 

. OOOOOOE+OO 

152 192E+12 

389 

. OOOOOOE+OO 

. 169364E+12 

390 

. OOOOOOE+OO 

188475E+12 
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391 

. OOOOOOE+OO 

.  209741E+12 

392 

. OOOOOOE+OO 

233408E+12 

393 

. OOOOOOE+OO 

. 259745E+12 

394 

. OOOOOOE+OO 

289053E+12 

395 

. OOOOOOE+OO 

. 321669E+12 

396 

.  OOOOOOE+OO 

357964E+12 

397 

.  OOOOOOE+OO 

. 398356E+12 

398 

.  OOOOOOE+OO 

443305E+12 

399 

.  OOOOOOE+OO 

. 493325E+12 

400 

.  OOOOOOE+OO 

548990E+12 

END  OF  PROBLEM  #  1 


E.  INPUT,  OUTPUT  DATA  OF  THE  AUGMENTED  HR  NOTCH  FILTER:  P  =  3 


INPUT  DATA  FOR  PROBLEM  #  1 


PROBLEM  #  1  INPUT  DATA  SOURCEFILE:  dif feq.  tst 

THE  NUMBER  OF  INPUT  DELAYS:  L  =  5 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =  5 

THE  VALUE  OF  nstop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  b(L)  ARE: 

. 990000E+00  329900E-02  .  186124E-01  .937440E-02  -.661007E+00  -.604706E+00 


THE  COEFFICIENTS  a(l),  ....  a(N)  ARE: 

. OOOOOOE+OO  . OQOQOOE+OO  . 000000E+00  . 652382E+00  . 598659E+00 


OUTPUT  DATA  FOR  PROBLEM  0  1 


ns  x(.ns) 

-5  . OOOOOOE+OO 
-4  . OOOOOOE+OO 
-3  . OOOOOOE+OO 
-2  . OOOOOOE+OO 
-1  .OOOOOOE+OO 
0  . 100000E+01 
1  . OOOOOOE+OO 


y(ns) 

. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. 990000E+00 
. 329900E-02 
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. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
.  OOOOOOE-i^ 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
.OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
.OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 


. 186124E-01 
. 937440E-02 
151485E-01 
141854E-01 
. 101674E-01 
. 172582E-01 
427053E-02 
183231E-01 
185917E-02 
. 173458E-01 
. 754575E-02 
145103E-01 
121822E-01 
. 102031E-01 
. 153069E-01 
494890E-02 
16634 IE-01 
636683E-03 
. 16Q941E-01 
. 593505E-02 
138145E-01 
10373SE-01 
. 101184E-01 
. 13S068E-01 
545927E-02 
150377E-01 
. 390820E-03 
. 148691E-01 
. 452445E-02 
130786E-01 
8747S0E-02 
. 993428E-02 
. 118532E-01 
562363E-02 
135363E-01 
. 12441 8E-02 
. 136801E-01 
. 3296/8E-02 
123172E-01 
729197E-G2 
. 966946E-Q2 
. 103405E-01 
-.606189.  -02 
-.  121310E-01 
. 194278E-02 
. 125346E-01 
. 223574E-02 
-. 115430E-01 
599489E-02 
. 934044E-02 
. 896254E-02 
-. 619203E-02 
108213E-01 
.250464E-02 
. 114387E-01 
. 132593E-02 
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58 

. 000000E+00 

107666E-01 

59 

. OOOOOQE+OO 

484430E-02 

60 

. OOOOOOE+OO 

. 896186E-02 

61 

. OOOOOOE+OO 

. 771293E-02 

62 

. OOOOOOE+OO 

623013E-02 

63 

. OOOOOOE+OO 

960583E-02 

64 

. OOOOOOE+OO 

.  294647E-02 

65 

. OOOOOOE+OO 

. 103969E-01 

66 

.  OOOOOOE+OO 

. 552991E-03 

67 

. OOOOOOE+OO 

999640E-02 

68 

.  OOOOOOE  00 

382839E-02 

69 

.  OOOOOOE+OO 

. 854667E-02 

70 

. OOOOOOE+OO 

. 658495E-02 

71 

. OOOOOOE+OO 

619042E-02 

72 

. OOOOOOE+OO 

848201E-02 

73 

. OOOOOOE+OO 

. 328379E-02 

74 

. OOOOOOE+OO 

. 941245E-02 

75 

. OOOOOOE+OO 

963783E-04 

76 

. OOOOOOE+OO 

923947E-02 

77 

. OOOOOOE+OO 

293555E-02 

78 

. OOOOOOE+OO 

. 810639E-02 

79 

. OOOOOOE+OO 

. 557198E-02 

80 

. OOOOOOE+OO 

608536E-02 

81 

. OOOOOOE+OO 

744640E-02 

82 

. OOOOOOE+OO 

. 353107E-02 

83 

. OOOOOOE+OO 

.848803E-02 

84 

. OOOOOOE+OO 

634268E-03 

85 

. OOOOOOE+OO 

-.850096E-02 

86 

. OOOOOOE+OO 

215424E-02 

87 

. OOOOOOE+OO 

. 765135E-02 

88 

. OOOOOOE+OO 

. 466765E-02 

89 

. OOOOOOE+OO 

592558E-02 

90 

. OOOOOOE+OO 

649457E-02 

91 

. OOOOOOE+OO 

. 370195E-02 

92 

. OOOOOOE+OO 

. 762565E-02 

93 

. OOOOOOE+OO 

1Q7141E-02 

94 

. OOOOOOE+OO 

778435E-02 

95 

. OOOOOOE+OO 

147295E-02 

96 

, OOOOOOE+OO 

. 719104E-02 

97 

.  OOOOOOE+OO 

. 386619E-02 

98 

.  OOOOOOE+OO 

571978E-02 

99 

.  OOOOOOE+OO 

56211QE-02 

100 

.  OOOOOOE+OO 

. 380951E-02 

101 

.  OOOOOOE+OO 

. 682722E-02 

102 

.  OOOOOOE+OO 

141695E-02 

103 

.  OOOOOOE+OO 

709131E-02 

104 

.  OOOOOOE+OO 

879866E*03 

105 

.  OOOOOOE+OO 

. 673456E-02 

106 

.  OOOOOOE+OO 

. 316279E-02 

107 

. OOOOOOE+OO 

-.547452E-02 

108 

.  OOOOOOE+OO 

481929E-02 

109 

. OOOOOOE+OO 

.386677E-02 

110 

.  OOOOOOE+OO 

. 609505E-02 

111 

.  OOOOOOE+OO 

167805E-02 

112 

.  OOOOOOE+OO 

*. 642139E-02 

113 

.  OOOOOOE+OO 

362501E-03 
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114 

. OOOOOOE+OO 

. 629118E-02 

115 

. 000000E+00 

. 255413E-02 

116 

. OOOOOOE+OO 

519378E-02 

117 

. OOOOOOE+OO 

40807  IE-02 

118 

. OOOOOOE+OO 

. 388724E-02 

119 

. OOOOOOE+OO 

. 543255E-02 

120 

. OOOOOOE+OO 

185927E-02 

121 

. OOOOOOE+OO 

577149E-02 

122 

. OOOOOOE+OO 

. 930098E-04 

123 

. OOOOOOE+OO 

. 587123E-02 

124 

. OOOOOOE+OO 

. 203929E-02 

125 

. OOOOOOE+OO 

487829E-02 

126 

. OOOOOOE+OO 

339448E-02 

127 

. OOOOOOE+OO 

. 388597E-02 

128 

. OOOOOOE+OO 

. 484526E-02 

129 

. OOOOOOE+OO 

196167E-02 

130 

. OOOOOOE+OO 

513493E-02 

131 

. OOOOOOE+OO 

. 503001E-03 

132 

. OOOOOOE+OO 

. 548734E-02 

133 

. OOOOOOE+OO 

. 162090E-02 

134 

. OOOOOOE+OO 

45243 IE-02 

135 

. OOOOOOE+OO 

274593E-02 

136 

. OOOOOOE+OO 

. 388097E-02 

137 

. OOOOOOE+OO 

. 434249E-02 

138 

. OOOOOOE+OO 

198121E-02 

139 

. OOOOOOE+OO 

449991E-02 

140 

. OOOOOOE+OO 

. 888000E-03 

141 

. OOOOOOE+OO 

. 515634E-02 

142 

, OOOOOOE+OO 

. 130717E-02 

143 

. OOOOOOE+OO 

412173E-02 

144 

. OOOOOOE+OO 

«. 211460E-02 

145 

. OOOOOOE+OO 

. 389552E-02 

146 

. OOOOOOE+OO 

. 393967E-02 

147 

. OOOOOOE+OO 

190640E-02 

148 

. OOOOOOE+OO 

384704E-02 

149 

. OOOOOOE+OO 

. 127544E-02 

150 

. OOOOOOE+OO 

. 490226E-02 

151 

.  OOOOOOE+OO 

. 111482E-02 

152 

. OOOOOOE+OO 

365103E-02 

153 

.  OOOOOOE+OO 

147099E-02 

154 

. OOOOOOE+OO 

. 396170E-02 

155 

. OOOOOOE+OO 

. 366207E-02 

156 

. OOOOOOE+OO 

171447E-02 

157 

.  OOOOOOE+OO 

314537E-02 

158 

.  OOOOOOE+OO 

. 170392E-02 

159 

. OOOOOOE+OO 

. 476078E-02 

160 

. OOOOOOE+OO 

. 107384E-02 

161 

. OOOOOOE+OO 

». 307837E-02 

162 

. OOOOOOE+OO 

771398E-03 

163 

. OOOOOOE+OO 

. 412592E-02 

164 

. OOOOOOE+OO 

.355064E-02 

165 

. OOOOOOE+OO 

136541E-02 

166 

. OOOOOOE+OO 

234614E-02 

167 

. OOOOOOE+OO 

. 222987E-02 

168 

. OOOOOOE+OO 

. 478639E-02 

169 

. OOOOOOE+OO 

. 123486E-02 
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170 

. OOOOOOE+OO 

234799E-02 

171 

. OOOOOOE+OO 

. 501888E-04 

172 

. OOOOOOE+OO 

. 445749E-02 

173 

. OOOOOOE+OO 

. 3671Q2E-02 

174 

. OOOOOOE+OO 

792S29E-03 

175 

. OOOOOOE+OO 

137291E-02 

176 

.OOOOOOE+OO 

. 293803E-02 

177 

. OOOOOOE+OO 

. 506343E-02 

178 

. OOOOOOE+OO 

. 168066E-02 

179 

. OOOOOOE+OO 

137011E-02 

180 

. OOOOOOE+OO 

. 109482E-02 

181 

. OOOOOOE+OO 

. 506217E-02 

182 

. OOOOOOE+OO 

. 412770E-02 

183 

. OOOOOOE+OO 

. 112304E-03 

184 

. OOOOOOE+OO 

105991E-03 

185 

. OOOOOOE+OO 

. 395790E-02 

186 

. OOOOOOE+OO 

. 572336E-02 

187 

. OOOOOOE+OO 

. 254435E-02 

188 

. OOOOOOE+OO 

191513E-05 

189 

. OOOOOOE+OO 

. 251861E-02 

190 

. OOOOOOE+OO 

. 610325E-02 

191 

. OOOOOOE+OO 

. 508623E-02 

192 

. OOOOOOE+OO 

. 152195E-02 

193 

. OOOOOOE+OO 

. 164195E-02 

194 

. OOOOOOE+OO 

. 548944E-02 

195 

. OOOOOOE+OO 

. 697193E-02 

196 

. OOOOOOE+OO 

. 403781E-02 

197 

. OOOOOOE+OO 

. 198231E-02 

198 

. OOOOOOE+OO 

. 456418E-02 

199 

. OOOOOOE+OO 

. 783467E-02 

200 

. OOOOOOE+OO 

. 680801E-02 

201 

. OOOOOOE+OO 

. 371050E-02 

202 

. OOOOOOE+OO 

. 416432E-02 

203 

. OOOOOOE+OO 

. 784359E-02 

204 

. OOOOOOE+OO 

. 913172E-02 

205 

. OOOOOOE+OO 

. 649634E-02 

206 

. OOOOOOE+OO 

. 493805E-02 

207 

. OOOOOOE+OO 

. 761003E-02 

208 

. OOOOOOE+OO 

. 106530E-01 

209 

. OOOOOOE+OO 

. 970489E-02 

210 

.  OOOOOOE+OO 

. 711060E-02 

211 

.  OOOOOOE+OO 

. 792086E-02 

212 

.  OOOOOOE+OO 

. 115057E-01 

213 

.  OOOOOOE+OO 

. 127088E-01 

214 

.  OOOOOOE+OO 

. 104488E-01 

215 

.  OOOOOOE+OO 

. 942425E-02 

216 

.  OOOOOOE+OO 

. 122480E-01 

217 

.  OOOOOOE+OO 

. 15 1790E-01 

218 

.  OOOOOOE+OO 

. 144248E-01 

219 

.  OOOOOOE+OO 

. 124035E-01 

220 

.  OOOOOOE+OO 

. 136323E-01 

221 

.  OOOOOOE+OO 

.172349E-01 

222 

.  OOOOOOE+OO 

. 184976E-01 

223 

.  OOOOOOE+OO 

. 167274E-01 

224 

.  OOOOOOE+OO 

. 163189E-01 

225 

.  OOOOOOE+OO 

. 194048E-01 
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226 

.  OOOOOOE+OO 

. 223853E-01 

227 

.  000000E+00 

. 219864E-01 

228 

.  OOOOOOE+OO 

. 206602E-01 

229 

.  OOOOOOE+OO 

. 224288E-01 

230 

.  OOOOOOE+OO 

. 262206E-01 

231 

.  OOOOOOE+OO 

. 277447E-01 

232 

.  OOOOOOE+OO 

. 266407E-01 

233 

.  OOOOOOE+OO 

.270006E-01 

234 

.  OOOOOOE+OO 

. 305331E-01 

235 

.  OOOOOOE+OO 

. 337974E-01 

236 

.  OOOOOOE+OO 

. 339895E-01 

237 

.  OOOOOOE+OO 

. 335634E-01 

238 

.  OOOOOOE+OO 

. 36QS34E-01 

239 

.  OOOOOOE+OO 

. 403277E-01 

240 

.  OOOOOOE+OO 

. 424073E-01 

241 

.  OOOOOOE+OO 

. 422443E-01 

242 

.  OOOOOOE+OO 

. 436332E-01 

243 

.  OOOOOOE+OO 

. 479108E-01 

244 

.  OOOOOOE+OO 

. 513083E-01 

245 

.  OOOOOOE+OO 

.529470;: -01 

246 

.  OOOOOOE+OO 

. 537555E-01 

247 

.  OOOOOOE+OO 

. 573776E-01 

248 

.  OOOOOOE+OO 

. 624811E-01 

249 

.  OOOOOOE+OO 

. 655572E-01 

250 

.  OOOOOOE+OO 

. 667663E-01 

251 

.  OOOOOOE+OO 

. 696133E-01 

252 

.  OOOOOOE+OO 

. 751112E-01 

253 

.  OOOOOOE+OO 

. 801732E-01 

254 

.  OOOOOOE+OO 

. 828036E-01 

255 

.  OOOOOOE+OO 

. 853848E-01 

256 

.  OOOOOOE+OO 

. 906759E-01 

257 

.  OOOOOOE+OO 

. 972696E-01 

258 

.  OOOOOOE+OO 

.  102016E+00 

259 

.OOOOOOE+OO 

.  105275E+00 

260 

.  OOOOOOE+OO 

.  110272E+00 

261 

.  OOOOOOE+OO 

.  117741E+00 

262 

.  OOOOOOE+OO 

.  124785E+00 

263 

.  OOOOOOE+OO 

.  129752E+00 

264 

.  OOOOOOE+OO 

.  134963E+00 

265 

.  OOOOOOE+OO 

.  142827E+00 

266 

.  OOOOOOE+OO 

.  151894E+00 

267 

.  OOOOOOE+OO 

.  159352E+00 

268 

.  OOOOOOE+OO 

.  165725E+00 

269 

.  OOOOOOE+OO 

.  173975E+00 

270 

.  OOOOOOE+OO 

.  184598E+00 

271 

.  OOOOOOE+OO 

.  19489 1E+00 

272 

.  OOOOOOE+OO 

.  203513E+00 

273 

.  OOOOOOE+OO 

.  2127 1 1E+00 

274 

.  OOOOOOE+OO 

.  224580E+00 

275 

.  OOOOOOE+OO 

.  237655E+00 

276 

.  OOOOOOE+OO 

.  249442E+00 

277 

.  OOOOOOE+OO 

.  260604E+00 

278 

.  OOOOOOE+OO 

.  273854E+00 

279 

. OOOOOOE+OO 

.  289489E+00 

280 

.  OOOOOOE+OO 

.  305006E+00 

281 

.  OOOOOOE+OO 

.  319344E+00 
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282 

. OOOOOOE+OO 

. 334670E+00 

283 

.  000000E+00 

. 352802E+00 

284 

.  OOOOOOE+OO 

. 372286E+00 

285 

.  OOOOOOE+OO 

. 390929E+00 

286 

.  OOOOOOE+OO 

. 409511E+00 

287 

.  OOOOOOE+OO 

. 430516E+00 

288 

.  OOOOOOE+OO 

. 454081E+00 

289 

.  OOOOOOE+OO 

. 477908E+00 

290 

.  OOOOOOE+OO 

. 501191E+00 

291 

.  OOOOOOE+OO 

. 526019E+00 

292 

.  OOOOOOE+OO 

. 553967E+00 

293 

.  OOOOOOE+OO 

. 583618E+00 

294 

.  OOOOOOE+OO 

. 613072E+00 

295 

.  OOOOOOE+OO 

. 643208E+00 

296 

.  OOOOOOE+OO 

.  676304E+00 

297 

.  OOOOOOE+OO 

.  712380E+00 

298 

.  OOOOOOE+OO 

.  749346E+00 

299 

.  OOOOOOE+OO 

.  786639E+00 

300 

.  OOOOOOE+OO 

.  82627 1E+00 

301 

.  OOOOOOE+OO 

.  869620E+00 

302 

.  OOOOOOE+OO 

.  915333E+00 

303 

.  OOOOOOE+OO 

.  961792E+00 

304 

.  OOOOOOE+OO 

.  100997E+01 

305 

.  OOOOOOE+OO 

.  106198E+G1 

306 

.  OOOOOOE+OO 

.  111775E+01 

307 

.  OOOOOOE+OO 

.  117543E+01 

308 

.  OOOOOOE+OO 

.  123467E+01 

309 

.  OOOOOOE+OO 

.  129745E+01 

310 

.  OOOOOOE+OO 

.  136497E+01 

311 

.  OOOOOOE+OO 

. 143598E+01 

312 

.  OCCOOOE+OO 

.  150916E+01 

313 

.  OOOOOOE+OO 

.  158558E+01 

314 

.  OOOOOOE+OO 

.  16672.'  E+OX 

315 

.  OOOOOOE+OO 

.  175396E+01 

316 

.  OOOOOOE+OO 

.  18442 1E+01 

317 

.  OOOOOOE+OO 

. 193788E+01 

318 

.  OOOOOOE+OO 

.  203688E+01 

319 

.  OOOOOOE+OO 

.  214234E+01 

320 

.  OOOOOOE+OO 

. 225316E+01 

321 

.  OOOOOOE+OO 

. 236829E+01 

322 

.  OOOOOOE+OO 

.  248895E+01 

323 

.  OOOOOOE+OO 

.  261702E+01 

324 

.  OOOOOOE+OO 

.  275245E+01 

325 

,  OOOOOOE+OO 

.  28939 1E+01 

326 

.  OOOOOOE+OO 

.  304155F.+01 

327 

.  OOOOOOE+OO 

.  319734E+01 

328 

.  OOOOOOE+OO 

. 336236E+01 

329 

.  OOOOOOE+OO 

. 353572E+01 

330 

.  OOOOOOE+OO 

. 371672E+01 

331 

.  OOOOOOE+OO 

. 390674E+01 

332 

.  OOOOOOE+OO 

. 410766E+O1 

333 

.  OOOOOOE+OO 

. 431955E+01 

334 

.  OOOOOOE+OO 

. 454141E+01 

335 

.  OOOOOOE+OO 

. 477374E+01 

336 

.  OOOOOOE+OO 

. 501857E+01 

337 

.  OOOOOOE+OO 

. 527708E+01 
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338 

. OOOOOOE+OO 

. 554867E+01 

339 

. OOOOOOE+OO 

. 583306E+01 

340 

. OOOOOOE+OO 

. 613187E+01 

341 

. OOOOOOE+OO 

. 644709E+01 

342 

. OOOOOOE+OO 

. 677903E+01 

343 

. OOOOOOE+OO 

. 712715E+01 

344 

. OOOOOOE+OO 

. 749234E+01 

345 

. OOOOOOE+OO 

. 787687E+01 

346 

. OOOOOOE+OO 

. 828213E+01 

347 

.  OOOOOOE+OO 

. 870796E+01 

348 

. OOOOOOE+OO 

. 915460E+01 

349 

. OOOOOOE+OO 

. 962409E+01 

350 

. OOOOOOE+OO 

. 101187E+02 

351 

. OOOOOOE+OO 

. 10639 1E+02 

352 

. OOOOOOE+OO 

. 111854E+02 

353 

. OOOOOOE+OO 

. 1175912+02 

354 

. OOOOOOE+OO 

. 123628E+02 

355 

. OOOOOOE+OO 

. 129984E+02 

356 

.  oooooof.+oc 

. 136663E+02 

357 

. OOOOOOE+OO 

. 143677E+02 

358 

- 00009QE+00 

. 15I049E+02 

3S9 

. OOOOOOE+OO 

. 158810E+02 

360 

. OOCOOOE+OQ 

. 166973E+02 

361 

. OOOOOOE+OO 

. 1755477, +02 

362 

.  OOOOOOE+OO 

. 184555E+02 

363 

.  OOOOOOE+OO 

. 194032E+02 

364 

.  OOOOOOE+OO 

. 204003E+02 

365 

.  OOOOOOE+OO 

. 214484E+02 

36o 

.  OOOOOOE+OO 

.  225493E+02 

367 

.  OOOOOOE+OO 

.  237069E+G2 

368 

.  OOOOOOE+OO 

.  249247E+02 

369 

.  OOOOOOE+OO 

262054S+02 

370 

.  OOOOOOE+OO 

.  27551 1E+02 

371 

.  OOOC'OOE+OO 

.  289653E+02 

372 

.  AQOOOOE+OO 

.  304528E+02 

373 

.  OOOOOOE+OO 

.  320174E+02 

374 

.  OOOOOOE+OO 

.  336619E+C2 

375 

.  OOOOOOE+OO 

.  353902E+C2 

376 

. OOOOOOE+OO 

.  372073E+02 

377 

.  OOOOOOE+OO 

.  391164E+02 

378 

. OOOOOOE+OO 

.  411279E+02 

379 

. oodbboE+oo 

.  432400E+02 

380 

, OOOOOOE+OO 

.  454600E+02 

381 

. OOOOOOE+OO 

.  477946E+02 

382 

. OOOOOOE+OO 

.  502498E+02 

383 

. OOOOOOE+OO 

. 528306S+O2 

384 

. OOOOOOE+OO 

. S55433E+02 

385 

. OOOOOOE+OO 

. 583954E+02 

386 

.OOOOOOE+OO 

. 613940E+O2 

387 

. OOOOOOE+OO 

. 645482E+02 

386 

. OOOOOOE+OO 

. 678630E+02 

389 

. OOOOOOE+OO 

. 713477E+02 

390 

. OOOOOOE+OO 

. 750118E+02 

391 

. OOOOOOE+OO 

. 788647E+02 

392 

. OOOOOOE+OO 

.829151E+02 

393 

. OOOOOOE+OO 

. 871728E+02 
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394 

. OOOOOOE+OO 

. 916494E+02 

395 

. OOOOOOE+OO 

.  963565E+02 

396 

. OOOOOOE+OO 

.  101305E+03 

397 

. OOOOOOE+OO 

.  106508E+03 

398 

. OOOOOOE+OO 

.  111977E+03 

399 

. OOOOOOE+OO 

.  117728E+03 

400 

. OOOOOOE+OO 

.  123775E+03 

END  OF  PROBLEM  it  1 
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